UNCLASSIFIED 


AD  NUMBER 

AD875392 

NEW  LIMITATION  CHANGE 
TO 

Approved  for  public  release,  distribution 
unlimited 


FROM 

Distribution  authorized  to  U.S.  Gov't, 
agencies  and  their  contractors; 
Administrative/Operational  Use;  Apr  1970. 
Other  requests  shall  be  referred  to  Air 
Force  Propulsion  Lab.,  Edwards  AFB,  CA. 


AUTHORITY 

AFAL  ltr,  18  Dec  1989 


THIS  PAGE  IS  UNCLASSIFIED 


Aerotherm  Report  No.  UM-70-14 


0> 

CO 

iO 

i> 

00 

a 

< 


AFRPLTR  70  92 


USER'S  MANUAL 

AEROTHERM  CHARRING  MATERIAL  THERMAL 
RESPONSE  AND  ABLATION  PROGRAM 
VERSION  3 


Volume  II-  Fortran  Variable  Names, 
Flow  Charts,  and  Listings 

April  1970 


Air  Force  Rocket  Propulsion  Laboratory 
Director  of  Laboratories 
Edwards,  California  93523 
Air  Force  Systems  Command 
United  States  Air  Force 


This  document  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 
be  made  only  with  prior  approval  of  AFRPL  ( REPORT/ST  I NFO ) , 
Edwards,  California  93523. 


AFRPL-TR-7Q-92 


Aerotherm  Report  No.  UM-70-14 
April  1970 


USER'S  MANUAL 

AEROTHERM  CHARRING  MATERIAL  THERMAL 
RESPONSE  AND  ABLATION  PROGRAM 
VERSION  3 

Volume  II  -  Fortran  Variable  Names, 
Flow  Charts,  and  Listings 


Prepared  Under 
the  Sponsorship  of 

Air  Force  Rocket  Propulsion  Laboratory 
Director  of  Laboratories 
Edwards,  California  93523 
Air  Force  Systems  Command 
United  States  Air  Force 

Project  Officer,  Robert  J.  Schoner/RPMCH 


This  document  is  subject  to  special  export  controls  and  each 
transmittal  to  foreign  governments  or  foreign  nationals  may 
be  made  only  with  prior  approval  of  AFRPL  (REPORT/STINFO) , 
Edwards,  California  93523. 


! 


FOREWORD 


This  report  is  one  of  two  computer  program  user's  manuals  prepared  by 
Aerotherm  Corporation  under  USAF  Contract  F04611-7Q-C-0012.  Included  herein 
is  Volume  II  of  the  manual  for  Version  3  of  the  Aerotherm  Charring  Material 
Ablation  code.  This  volume  presents  definitions  of  Fortran  variables,  flow 
charts,  and  program  listings.  The  report  was  first  published  as  Aerotherm 
Report  No.  UM-70-14.  The  work  was  administered  under  the  direction  of  the 
Air  Force  Rocket  Propulsion  Laboratory,  Motor  Component  Development  Branch 
with  Mr.  R.  J.  Schoner  as  project  officer. 

Mr.  M.  R.  Wool  was  program  manager  and  principal  investigator.  Sig¬ 
nificant  additional  assistance  was  also  provided  by  Dr.  C.  B.  Moyer. 

This  technical  report  has  been  reviewed  and  is  approved. 


R.  J.  Schoner 
Project  Engineer,  AFRPL 
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ABSTRACT 


This  two-volume  report  describes  a  Fortran  IV  computer  code  which  com¬ 
putes  the  transient  thermal  and  ablation  response  of  a  charring  insulation 
material  structure.  The  program  is  for  one-dimensional  bodies,  but  can  treat 
a  variety  of  shapes,  including  planes,  cylinders,  spheres,  and  more  general 
thermal  "stream  tube”  bodies.  The  program  can  treat  complex  systems  including 
a  main  ablating  material,  several  charring  back-up  materials,  and  a  multiple 
non-charring  material  back-up  structure. 

An  unusual  feature  of  the  code  is  the  very  general  heated  surface  bound¬ 
ary  conditions,  which  can  account  for 

•  Simple  specified  temperature  and  recession  rate 

•  Specified  heat  flux  with  no  recession 

•  General  thermochemical  erosion  model  incorporating  complete 
chemical  erosion  computations,  both  equilibrium  and  non¬ 
equilibrium,  for  any  material  exposed  to  any  environment. 

The  code  has  seen  extensive  use  for  thermal  performance  studies  of 
ablating  heat  shields,  rocket  nozzles,  and  spacecraft  structures. 

Volume  I  of  this  report  contains  descriptions  of  the  problem  treated, 
the  equations  solved,  the  input  information  required  of  the  program  user,  and 
the  program  output  information.  It  also  provides  a  card-by-card  user’s  input 
guide  and  a  number  of  sample  problem  input  and  output  listings.  Volume  II  of 
the  report  contains  supplemental  information  on  the  specific  Fortran  IV  codings. 
It  includes  program  listings,  flow  charts,  and  definitions  of  Fortran  variable 


names. 
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When  U.S.  Government  drawings,  specifications,  or  other  data 
are  used  for  any  purpose  other  than  a  definitely  related 
Government  procurement  operation,  the  Government  thereby  Incurs 
no  responsibility  nor  any  obligation  whatsoever,  and  the  fact 
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way  supplied  the  said  drawings,  specifications,  or  other  data. 
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or  conveying  any  rights  or  permission  to  manufacture,  use,  or 
sell  any  patented  Invention  that  may  In  any  way  be  related 
thereto. 
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SECTION  1 


INTRODUCTION 

The  computer  program  described  in  this  user's  manual  is  a  revised 
edition  of  the  Aerotherm  Charring  Material  Ablation  coaiputer  program.  The 
previous  code,  CMA  Version  2,  is  described  in  References  1,  2,  and  3.  The  cur¬ 
rent  program,  CMA  Version  3  (denoted  CMA 3  or  CMA/CBM) ,  solves  all  problems  that 
earlier  versions  could  solve  and  provides  additional  computational  capabilities, 
featuring  in  particular  an  added  charring  back-up  material  capability. 

The  purpose  of  Volume  I  of  this  user's  manual  was  to  enable  an  un¬ 
familiar  user  to  utilize  effectively  Version  3  of  the  Aerotherm  Charring  Mate¬ 
rial  Ablation  computer  program.  It  contained  a  general  description  of  the  prob¬ 
lems  that  CMA 3  solves,  an  input  data  deck  preparation  guide,  and  sample  problem 
input  and  output.  Volume  II  of  this  manual,  included  herein,  contains  the 
following  additional  program  documentation: 

uM 

e  Definitions  for  important  Fortran  variables  used 

e  Flow  charts  of  program  logic  for  each  Fortran  routine 

e  Listings  of  Fortran  IV  source  decks 

These  are  given  in  Sections  2,  3,  and  4,  respectively.  Another  document  of 
potential  value  to  a  reader  desiring  a  more  detailed  exposition  of  the  theoreti¬ 
cal  fundamentals  of  the  CMA  program  is  Reference  4. 


SECTION  2 


FORTRAN  VARIABLE  NAMES 

The  names  and  definitions  of  all  Fortran  variables  used  in  impor¬ 
tant  subroutines  in  Version  3  of  the  Aerotherm  Charring  Material  Ablation  code 
are  given  in  this  section.  Since  the  CMA3  code  is  divided  into  several  routines, 
the  majority  of  all  variables  are  included  in  COMMON  statements.  These  variables 
are  defined  first.  Fortran  variables  which  are  used  only  locally  in  the  CBM, 
INP0UT,  and  THERMS  subroutines  are  then  defined.  Variables  used  locally  in 
minor  subroutines  are  not  given. 


LIST  UF  FORTRAN  VARIABLES 

APPEARING  in  common  statements 


VARIABLE 

NAME 

DESCRIPTION 

units 

AREA  (  ) 

ORIGINAL  (INPUT)  CRosS  SECTION  AREA  OF  NOOE 

FT2 

ASTER 

1H* 

6a 

ARRHENIUS  PKE-tKPONtNTIAL  FACTOR  IN  DECOMPOSITION 
KINETIC  EuUATIUN  FOR  COMPONFNT  A 

SEC-1 

8b 

ARRHENIUS  PRt -EXPONENTIAL  FaCTuR  IN  DECOMPOSITION 
KINETIC  EuUATION  FOR  COMPONEN I  d 

SEC-1 

HuBII.J) 

ARMEN I US  PRE-EaPOnEnTIAL  FACTOR  IN  DECOMPOSITION 
KINETIC  EuUATION  OF  J  Trt  component  OF  I  IN 

DECOMPOSING  HACK-Um  (ANALOGOUS  TO  BA.  Bb*  BC) 

SEC-1 

flC 

ARRHENIUS  PNE -EXPONENTIAL  FACTOR  IN  DECOMPOSITION 
KINETIC  EUUATION  FOk  COMPONENT  C 

SEC-1 

BLANK 

1H 

BKEA 

INPUT  BURNING  MATE  EXPONENT  USED  IN  MODIFYING 

INPUT  TRAnSFEP  COEFFICIENT 

BMP 

BLOWING  REDUCTION  PARAMETER.  LAMBDA 

CnC«1 

\ 

Cmm 

CHAW  INTERFACE  CRlTtRlAl  CHAN  LINE  DEPTH  IS  DEPTH  OF 
LINE  OF  OtNSITf 

RHU(2)  *lhCRI*  (RmU  (  1 )  -RHO(<?)  ) 

RATIO  OF  CM/CH 

OtL  (  1 

NODE  THICKNESS 

FT 

OtLMO 

ENTHALPY  uF  FORMAT Iuh  OF  PYROLYSIS  GASES 

BTU/Ltf 

OtLM 

minimum  allukeo  Thickness  of  last  (shrinking) 

MODE  IN  ActLATlNG  MAItRlAL 

FT 

OtN 

INTERPOLATION  FRACTION  FUR  FUNCT IONS-OF-1 IME 

TABLES 

U«C ( l ) 

ENTHALPY  uf  1 uRMAl Iuh  OF  char  OF  1  TH  DECOMPOSING 
HACK-uP  (ANALOGOUS  lu  UH?) 

BTU/LB 

OMV ( I ) 

enthalpy  uf  formation  of  virgin  plastic  of  i  th 

DECOMPOSING  BACK-UP  (ANALOGOUS  TO  DH1) 

BTU/LB 

On  1 

ENTHALPY  OF  FORMATION  UF  VIRGIN  PLASTIC*  MAIN 

material 

BTU/LB 

« 
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LIST  Of  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 


VARIABLE 

name 

DESCRIPTION 

UNITS 

OH^ 

ENTHALPY  OF  FORMATION  OF  CHaR*  MAIN  MATERIAL  , 

BTU/LB 

OTMB 

LIMIT  VALUE  OF  TIME  ^NCREM^NT 

SEC 

DTHIN 

INITIAL  TIME  STEP  USED  AT  BEGINNING  OF  PROBLEM 

ANO  AFTER  A  CHANGE  lH  OPTION  , 

SEC 

OTPRT 

FIRST  OUTPUT  TJMt  INTERVAL 

\ 

SECOND  OUTPUT  TIME  INTERVAL  ^ 

SEC  * 

0TPR2 

SEC 

DTPR3 

THIRD  OUTPUT  TIME  INTERVAL  ,  •  t 

✓  j  ‘  .  ' 

SEC 

EA 

activation  energy  for 'component  a*  EA/R 

■0 

activation  energy  f- or  component  b.  EB/R  ■  \ 

T  DEG  R 

Eu 

.  '  DEG  R 

EC  - 

ACTIVATION  energy. for  component  c*  ec/R 

^  DEG  R 

EE(I*J) 

activation  energy  E(I*J)/R  for  component  j  of  1  th 

DECOMPOSING  BACK-UP  (ANALOGOUS  TO  FA*  E0.EO 

DEG  R 

EMA<  ) 

slope  of  area  vs  absolute  Radius  (or  distance) 

CURVE  OBTAINED  FROM  ORIGINAL  INPUT  OF  NODAL  DATA 

'  Ft 

EPSV 

EMISSIVITV  OF  BACK  HALL 

ETA 

NOT  CURRENTLY  USED  / 

A 

Ff (I.J) 

PRECALCULATED  COLLELUON  OF  TERMS  IN  DECOMP¬ 
OSITION  KINETIC  EQUATION  FOR  J  TH  COMPONENT  OF 

1  TH  DECOMPOSING  BACK-UP  (ANALOGOUS  TO  FA.  FB.  FC) 

(L8/FT3)** 

(t.-PSl(l.J) 

/SEC 

FJFM 

HALF  THE  STANOARO  NUMBER  OF  NOOELETS  (SEE  JF  ANO 
FJFS) 

* 

FJFS 

STANOARO  NUMBER  OF  NOOELETS  PER  NOOC  FOR  ALL  NOOES 
EXCEPT  THE  FIRST*  HHICH  HAS  HALF  THIS  MANY 

F1II.LI 

L  TH  TABULAR  F-FUNCUON  FOR  VIRGIN  CONDUCTIVITY 

IN  COMPUTING  CONDUCTIVITY  In  PARTIALLY  „ 

PYROLY/ED  /ONE 

FZ(I*L) 

L  TH  TABULAR  F-  FUNCTION  FOR  CHAN  CONDUCTIVITY* 

USED  IN  COMPUTING  CONDUCTIVITY  IN  P ART J ALLY 

pvroly/eu  zone 

6A(I) 

RESIN  VOLUME  FRACTION  FUR  I  TH  DECOMPOSING  BACKUP 
(ANALOGOUS  TO  GAMA)  C 

FT 3  RESIN/ 
FT3  MATERIAL 
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LIS?  UF  FURTRlN  VARIABLES 
APPEAR 1N6  IN  COMMON  STATEMENTS 


VARIABLE  DESCRIPTION  UNITS 

NAME 


GAMA  VOLUME  FRACTION  OF  RESIN  (COMPONENTS  A  AND  •)  IN.  FT 3-RES IN/ 

UNOECOMPOsED  PLASTIC*  MAIN  MATERIAL  FT3-MATERIAL 

M (  )  '  ENTHALPY  AT  a  NODE  BTU/LB 

HCONV  CONVECTIVe.  HEAT  TRANSFER  COEFFICIENT  AT  BTU/FT2-SEC- 

BACK  HALL  OE6  F 

DUMMY  EMPTY  NONOS  In  COMMON 

/  _  _  • 

NOT  USEO  IN  PROGRAM  (RETURNED  FROM  LOOK  SUBROUTINE 

VIA  COMMON;  SET  GRtATER  THAN  UNITY  IF  EXTRAP¬ 
OLATION  REQUIRED. ) 

IHI (  )  INOEX  of  last  entry  in  a  table 

II  OPTION  INDEX 

ILO<  I  INOEX  OF  FIRST  ENTRV  IN  A  TABLE 

INCH  INPUT  TAPE  UNIT  NUMBER  FOR  SURFACE  EQUILIBRIUM 

DATA 

input  *  Input  tapl  unit  number  variable 

]K(  I  ■REMEMBEHEO*  INOEX  IN  a  GIVEN  TABLE  ADJACENT  TO 

\  PREVIOUS  VALUE  FOR  •rilCH  A  LOOK-UP  HAS  PERFORMED 

IsfcNIIPP)  NUMBER  OF  ENTRIES  In  EDGE  TABLE  FOR  INDICATED 
\  PRESSURE 

*  s.  •  ‘ 

JF  STANDARD  NUMBER  of  nooelets  per  nooe 

JFh  HALF  THE  sTANOARO  number  of  nooelets  per  nooe 

JFH  ♦  1 

V  _  ># 

MARKS  LAST  ENTRY  OF  nO-ABLATION  (TEMPERATURE 
1NOEPENOENT)  PART  OF  A  CHAR  RATE  TABLE  FOR  I  TH 
GAS  NATE  A NO  j  TH  PkcSSURE  (REGARDLESS  OF  USERS 
ORIGINAL  INTENTION*  THIS  IS- TAKEN  AS  LAST  TEMPER¬ 
ATURE  BEFORE  TEHPERAlURE  ENTRIES  BE6IN  TO  OESCENO* 

IF  EVER*  THIS  EXTENSION  SMOOTHS  INTERPOLATION.) 

A.  , 

OUTPUT  TAPE  UNIT  NUMBER  VARIABLE 
LET  NUMBER  OF  LINES  REMAINING  ON  CURRENT  OUTPUT  PAGE 

MA7L  (  )  MATERIAL  IDENTIFICATION  NUMBER  FOR  A  NODE 


JFHP 

KHIU.J) 


KUUT 


IDIOT (  ) 
IEX 


VARIABLE 

NAME 

NUN 

NUN? 

MBPF 

NHUFT(J) 

NCON 

NOBU 

NFI(I) 

NF  IS 

MiUl.j) 

N1 

NL 

NLAUI 

NL0C1.J) 

NHG 

NN 

NO 


LIST  (If  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 

UC SCRIPT ION 


ALVArS  EQUAL  TO  1  PLUS  THE  ORIGINAL  NUMBER  OF 
MOOES  OF  ABLATING  MATERIAL 

EQUAL  TO  MM  IF  THEME  ARE  NO  DECOMPOSING  BACK-UPS* 
OTHERWISE  EQUAL  TO  ONE  GREATER  THAN  THE  ORIGINAL 
NUMBER  OF  NOOES  IN  THE  MAIN  ABLATING  MATERIAL  ANO 
THE  DECOMPOSING  BACK-UPS  / 

INPUT  FLAG  INTEGER  CALLING  FOR  READING  OF 

TBPF  (  •  •  I  VALUES  FROM  SURFACE  THERMOCHEMISTRY 

TABLES 

INOEX  IDENTIFYING  THE  NUMBER  OF  THE  F -FUNCTION 
TABLE  (NUMBERED  IN  UNDER  INPUT)  ASSIGNED  TO 
J  TH  CHARRING  BACK-UP 

INPUT  FLAG  INTEGER  CALLING  FOR  OUTPUT  OF  NOOAL 

thermal  conductivity  in  place  of  nooal  enthalpy  in 

STANDARD  OUTPUT  BLOCK 

NUMBER  OF  DECOMPOSING  BfCX-uP  MATERIALS 
NUMBER  OF  THE  FIRST  m60E  IN  THE  I  TH  DECOMPOSING 

back-up  material 

INPUT  FLAG  INTEGER  galling  FOR  USE  OF  FISSURE 
MOOEL  IN  BLOVldG-  REDUCTION  OF  CONVECTIVE 
TRANSFER  COEFFICIENTS 

MARKS  TOP  ENTRY  IN  ABLATION  PART  OF  A  CHAR  RATE 
TABLE  FOR  I  TH  GAS  RATE  ANO  J  TH  PRESSURE 

NUMBER  OF  ISOTHERM  POINTS  CALLED  FOR 

NUMBER  OF  LAST  NOOE  OF  ABLATING  MATERIAL 

NUMBER  OF  THE  LAST  NODE  IN  THE  I  TH  DECOMPOSING 
BACKUP  MATERIAL* 

MARKS  BOTTOM  ENTRY  IN  ABLATION  PART  OF  A  CHAR 
RATE  TABLE  FOR  I  TH  GAS  RATE  ANO  J  TH  PRESSURE 

TOTAL  NUMBER  OF  GAS  KATE  ENTRIES  IN  EACH  PRESSURE 
GROUP  OF  SURFACE  EQUILIBRIUM  UATA 

FLAG  CALLING  FOR  PUNCH  OUTPUT  OF  THERMOCOUPLE 
ANO  ISOTHERM  OATA  / 

NUMBER  OF  THERMOCOUPLE  TEMPFRATURES  CALLED  FOR 
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LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 


variable 

NAMt 

INSCRIPTION 

UNITS 

NO  I 

NO  ♦  Nl 

- ... 

NPG 

CURRENT  PAGE  NUMBER  FOR  OUTPUT  LISTING*  USED 

IN  COMMUNICATION  WITH  SUBROUTINE  LCOUNT 

NPR 

NUMBER  OF  PRESSURES  IN  SURFACE  THERMOCHEMISTRY 

OATA  OECK 

NR 

FLAG  VARIABLE  CALLING  FOR  BURNING  RATE  OR  RAOIUS 
RATIO  CORRECTION  TO  INPUT  CONVECTIVE  TRANSFER 
COEFFICIENTS 

/ 

ALWAYS  EUUAL  TO  THE  ORIGINAL  TOTAL  NUMBER 

OF  NODES 

NUMN 

OMG 

1.  -  GAMA 

FT3-REIN- 

F0RCEMENT/FT3 

MATERIAL 

OMGA(I) 

1.  — ' GA ( l )  - 

FT3-REIN- 

F0RCEMENT/FT3 

MATERIAL 

P  ( I ) 

f 

RH0V(I)/(kH0V(I)-RM0CUM  *  ANALOGOUS  TO  PETE 

PtT 

DENSITY  FACTOR.  (CHAh  DENSITY  •  PETE) 

LB/FT3 

PETE 

DENSITY  TERM,  VIRGIN  OENSITy/  (VIRGIN  DENSITY  - 
CHAR  OENSITY) 

PP(I) 

RmOC ( I ) *P l I ) *  ANALOGOUS  TO  PET 

LB/FT3 

PSI <l,J) 

DECOMPOSITION  REACTION  ORDER  FOR  THE  J  TH  COMPON¬ 
ENT  OF  THE  I  TH  DECOMPOSING  BACA-UP  MATERIAL 
(ANALOGOUS  TO  PSIA,  PSIB*  PSlC) 

PSIA 

DECOMPOSITION  REACTION  ORQER  FOR  A  COMPONENT 

PSIB 

DECOMPOSITION  REACTION  ORDER  FOR  B  COMPONENT 

PSIC 

DECOMPOSITION  REACTION  ORDER  FOR  C  COMPONENT 

PyCRI 

PYROLYSIS  INTERFACE  CRITERION!  PYROLYSIS  LINE  OEPTH 
IS  DEPTH  OF  LINE  OF  OENSITY 

RHO (?) ♦ PYCRI* (RHO ( I ) —RHO (2) ) 

RA(  ) 

CURRENT  LOCATION  Ur  NODE.  MEASURED  FROM  ORIGINAL 
LOCATION  OF  HEATEO  SURFACE 

IN 
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variable 

name 


LIST  OF  FORTRAN  VARIABLES 
APPEAR 1 NO  IN  COMMON  STATEMENTS 

DESCRIPTION 


UNITS 


RAV  (  ) 

RCI  ) 

RECORD (  > 
RFT 
RHO(  ) 

RHOC ( I ) 


RHOO(I.J) « 
J*  1*2 


RHOQ(I,J) 

RHOOA 

RHOOB 

RHOOC 


RmOKB 

RhORC 


SAVED  ORIGINAL  VALUES  OF  NODAL  LOCATIONS  RA(  »  IN 

CONTACT  RESISTANCE  BETWEEN  INDICATED  NODE  AND  (FT2-SEC- 

AND  NEXT  NODE  DOWN  DEG  Ft/ 8TU 

ALPHAMERIC  TITLING  INFORMATION 

NOT  CURRENTLY  USED 

OENSITY  OF  VIRGIN  PLASTIC  ())•  PURE  CHAR  (2) »  LB/FT3 

OR  BACK-UP  MATERIALS  13-10) 

CHAR  OENSITY  OF  I  TH  DECOMPOSING  BACK-UP  MATERIAL  LB/FT3 

(ANALOGOUS  TO  RH0(2)) 

ORIGINAL  DENSITY  OF  J  TH  COMPONENT  OF  1  TH  DECOM^  L8/FT3 

POSING  BACK-UP  PER  UNIT  VOLUME  OF  RESIN  (ANALOGOUS 
TO  RHOOA f  RHOOB) 

ORIGINAL  DENSITY  OF  THIRD  COMPONENT  OF  I  TH  OECOM-  LB/FT3 
POSING  BACK-UP  FOR  UNIT  VOLUME  OF  REINFORCEMENT 
(ANALOGOUS  TO  RHOOC) 

ORIGINAL  OENSITY  OF  COMPONENT  A  PER  UNIT  VOLUME  LB/FT3 

OF  RESIN 

ORIGINAL  OENSITY  OF  COMPONENT  B  PER  UNIT  VOLUME  LB/FT3 

OF  RESIN 

ORIGINAL  OENSITY  OF  COMPONENT  C  PER  UNIT  VOLUME  LB/FT3 

OF  REINFORCEMENT 

RESIDUAL  OENSITY  OF  J  TH  COMPONENT  OF  I  TH  OECOM-  LB/FT 3 

POSING  BACK-UP  PER  UNIT  VOLUME  OF  RESIN  (ANALOGOUS 
TO  RHORA*  RHORB) 

RESIDUAL  OENSITY  Of  TH1R0  COMPONENT  OF  I  TH  OECOM-  LB/FT3 
POSING  BACK-UP  PER  UNIT  VOLUME  OF  REINFORCEMENT 
(ANALOGOUS  TO  RHORC) 

RESIDUAL  DENSITY  OF  COMPONENT  A  PER  UNIT  VOLUME  LB/FT3 

OF  RESIN 

RESIDUAL  DENSITY  OF  COMPONENT  B  PER  UNIT  VOLUME  LB/FT3 

OF  RESIN 

RESIDUAL  OENSITY  OF  COMPONENT  C  PER  UNIT  VOLUME  LB/FT3 

OF  REINFORCEMENT 
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VARIABLE 

name 

LIST  OF  FORTRAN  VARIABLES 

APPEAR I No  IN  COMMON  STATEMENTS 

INSCRIPTION 

UNITS 

RMOV ( I ) 

VIROlN  DENSITY  OF  1  TH  DECOMPOSING  BACK-UP 

MATERIAL  (ANALOGOUS  TO  RHO(t)) 

LB/FT 3 

RuA(  > 

CURRENT  DENSITY  Of  COMPONENT  A  PER  UNIT  VOLUME 

OF  RESINt  AT  A  NOOELET 

L8/FT3 

RUB  (  ) 

CURRENT  DENSITY  OF  COMPONENT  B  PER  UNIT  VOLUME 

OF  RESINt  AT  A  NOOELET 

L8/FT3 

ROC  (  ) 

CURRENT  DENSITY  OF  COMPONENT  C  PER  UNIT  VOLUME 

OF  REINFORCEMENT 

LB/FT 3 

ROCOMt I » J) • 

J  «  1*2 

DENSITY  OF  THE  J  TH  COMPONENT  IN  THE  I  TH  NODE 
(BACK-UPS  ONLY)  PER  UNIT  VOLUME  OF  RESIN  (ANAL¬ 
OGOUS  TO  MOA(I)t  ROB ( I ) ) 

LB/FT3 

ROCOMt I, 3) 

DENSITY  OF  THE  THIRD  COMPONENT  IN  THE  1  TH 

NODE  (BACK-UPS  ONLY)  PER  UNIT  VOLUME  OF  RESIN 
(ANALOGOUS  TO  ROC (I) ) 

LB/FT3 

RSV 

ORIGINAL  SURFACE  RADIUS  (NEGATIVE  IF  EXTERNAL) 

IN 

SU(  ) 

THERMOCOUPLE  DEPTHS  ANO/OR  ISOTHERM  TEMPERATURES 

IN  OR  OEG  1 

SMELL 

PROPORTIONALITY  FACTOR  IN  ChAH  SMELL  CORRECTION 

of  heat  transfer  coefficient 

TA(  1 

NODAL  TEMPERATURE 

DEG  K 

TBPF ( 1 •  J*K) 

TABULAR  VALUE  OF  B  PRIME  FAIL  (THAT  ISt  MOOT 
FAIL/RHOE-Ut-CM)  IN  SURFACE  THERMOCHEMISTRY  TABLE 

FOR  1  TH  ENTRY t  J  Tn  PYROLYSIS  GAS  RATE*  AND  K  TH 
PRESSURE 

TBRP(I) 

1  TH  ENTRY  IN  TINE  TABLE  VALUES  OF  BLOVING 

REDUCTION  PARAMETER  LANOA 

TCBUtl.J) 

TABULAR  VALUE  OF  SPECIFIC  HEAT  FOM  I  TM  TEMPER¬ 
ATURE  ENTRY  AND  J  Th  MATERIAL  INOCX  NUMBER  FOR 
DECOMPOSING  BACK-UPS  (ANALOGOUS  TO  THZlItJ)) 

BTU/LB 

OEG 

TCMEM 

(ltJ.K) 

QUANTITY  (CM/CM  •  CMEM  PROO  -  MB)  FOR  I  TH  CHAR 

RATE*  J  Th  GAS  RATE*  AND  K  TH  PRESSURE 

BTU/LB 

TCMC  ) 

TABULAR  VALUE  OF  HEAT  TRANSFER  COEFFICIENT 

RHUE-UE-CM  In  FUNCTIUNS-UF-TIME  TABLES 

LB/FT2-SEC 

TCP t l » J) 


TABULAR  VALUE  Of  SPECIFIC  HEAT  FOR  I  TM  TEMPER 
AYURE  ENT*Y  AND  J  TH  MATERIAL 


BTU/LB-OCS  R 


.  ■uup.'i 


VARIABLE 

NAME 

TCPSEN(I.J> 
TtNT (I *J) 

TtP(I*J> 

THE  (  I 

THFIN 
THG<  > 

TMSENUtJI 

THZtl.JI 

TMZRO 

TKBU(I*J) 

TRPdtJ) 
TLMCd  t  J»K> 

TMOd.J) 

THBT 
TPK  I 

TPRI  ) 


LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 

OESCR1PTION 


SLOPE  OF  THSEN  VERSUS  TV  AT  I  TH  TEMPERATURE 
AND  J  TH  PRESSURE  IN  EDGE  TABLE 

TABULAR  VALUE  Of  SPECIFIC  HEAT  FOR  I  TH  TEMPERATURE 
ENTRY  AND  J  TH  MATERIAL  INDEX  NUMBER  FOR  DECOMP¬ 
OSING  BACK-UPS  (ANALOGOUS  TO  THZ(1*J>> 

TABULAR  VALUE  Of  EMisSIVlTY  OF  I  TH  TEMPERATURE 
ENTRY  ANO  J  TH  MATERIAL 

TABULAR  VALUE  Of  BOUNDARY  RECOVERY  ENTHALPY  IN 
f UNCTION-OP-TIME  TABLES  (EOuAlS  SURE ACE  TEMPERATURE 
IN  OPTION  2  ANO  VIE*  FACTOR  IN  OPTION  3) 

TIME  AT  ENO  OF  PROBLEM 

TABULAR  VALUE  OF  PYROLYSIS  GAS  ENTHALPY  AS  A 
FUNCTION  Of  TEMPERATURE 

VALUE  OF  HER  AT  1  TH  TEMPERATURE  ANO  J  TH  PRESSURE 
IN  EOGE  TABLES 

TABULAR  VALUE  Of  SENSIBLE  ENTHALPY  Of  I  TH  TEMP¬ 
ERATURE  entry  ano  j  th  material 

INITIAL  PROBLEM  TIME 

TABULAR  VALUE  Of  THERMAL  CONDUCTIVITY  FOR  I  TH 
TEMPERATURE  ENTRY  ANO  J  TH  HATEMIAL  INDEX  NUMBER 
f OR  DECOMPOSING  BACK-UPS  (ANALOGOUS  TO  TKP(ItJ)) 

TABULAR  VALUE  Of  THERMAL  CONDUCTIVITY  FOR  I  TH 
TEMPERATURE  ENTRY  ANO  J  TH  TEMPERATURE 

TABULAR  VALUE  OF  LN(HOOTC/RMO£-UE-CM>  IN  SURFACE 
THERMOCHEMISTRY  TABU  FOR  1  TH  CHAR  RATE*  J  TH 
PYROLYSIS  GAS  RATE*  ANO  K  Th  PRESSURE 

TABULAR  VALUES  OF  GAS  RATES  APPEARING  IN  SURFACE 
THERMOCHEMISTRY  TABLES*  I  GIVES  GAS  RATE  ENTRY* 

J  GIVES  PRESSURE  INDEX 

NOT  CURRENTLY  USEO 

TABULAR  VALUE  OF  LN  OF  PRESSURE  IN  FUNCTIONS- 
Of-TIME  TABLE 

TABULAR  VALUE  OF  LN  OF  PRESSURE  IN  SURFACE  THERMO¬ 
CHEMISTRY  tables 


UNITS 

BTU/LB-OEG  R 
BTU/LB 


BTU/LB  OR 
OEG  R  OR  — 

SEC 

BTU/LB 

BTU/LB 

BTU/LB 

SEC 

BTU/FT-SEC- 
DEG  F 

BTU/FT-SEC- 
DEG  F 


LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 


VARIABLE 

NAMt 

DESCRIPTION 

UNITS 

TPR2 

TIME  OF  TRANSITION  FROM  INITIAL  TO  SECOND 

OUTPUT  TIME  INTERVAL 

Sec 

TPR3 

TIME  OF  TRANSITION  FROM  SECOND  TO  THIRD 

OUTPUT  TIME  INTERVAL 

SEC 

TwlR  (  ) 

TABULAR  VALUE  OF  RAUIANT  HEaT  FLUX  IN  FUNCTIONS- 
-OF-TIME  TABLE,  FOR  OPTIONS  1  AND  3  (OPTION  2* 
SURFACE  RECESSION  RATE) 

BTU/FT2-SEC 
OR  HILS/SEC 

TKAC(I.J) 

LIMIT  TEMPERATURE  BELOB  WHICH  COMPONENT  J  OF 
DECOMPOSING  BACK-UP  I  IS  NOT  ALLOWED  TO  DECOMPOSE 
(EFFICIENCY  MEASURE)  (ANALOGOUS  TO  TRACA*TRACB* 
TRACC) 

DEG  R 

TWACA 

LIMIT  TEMPERATURE  BLLOW  WHICH  COMPONENT  A  IS 

NOT  ALLOWEO  TO  OECOMPOSE  (EFFICIENCY  MEASURE) 

OEG  R 

TRACB 

LIMIT  TEMPERATURE  BELOW  WHICH  COMPONENT  B  IS 

NOT  ALLOWED  TO  DECOMPOSE  (EFFICIENCY  MEASURE) 

DEG  R 

tkacc 

LIMIT  TEMPERATURE  BELOW  WHICH  COMPONENT  C  IS  NOT 
ALLOWED  TO  DECOMPOSE  (EFFICIENCY  MEASURE) 

DEG  R 

TRACm 

MINIMUM  OF  TRACA*  TNACB*  TRACC 

TMEFd) 

REFERENCE  TEMPERATURE  FOR  HEATS  OF  FORMATION 
ANALOGOUS  TO  TZ>  FOR  I  TH  DECOMPOSING  BACK-UP 
MATERIAL* 

DEG  R 

TR£S 

RESERVOIR  TEMPERATURE  TO  WHICH  BACK  WALL  IS 

EXPOSED 

OEG  R 

TTH(  > 

TABULAR  VALUE  OF  TIME  (lNOEPENOENT  VARIABLE)  IN 

FUNCT IONS-OF-T INC  TABLE 

TTS(1*J*K) 

TABULAR  VALUE  OF  TEMPERATURE  IN  SURFACE  THERMO¬ 
CHEMISTRY  TABLES  FOR  1  TH  TEMPERATURE*  J  TH 
PYROLYSIS  GAS  RATE*  ANO  K  TH  PRESSURE 

OEG  R 

ttsenii • J) 

TABULAR  VALUE  OF  TEMPERATURE  (INDEPENDENT  VARI¬ 
ABLE)  IN  EOGE  TABLES  FOR  1  TH  TEMPERATURE  ENTRY 

IN  J  TH  PRESSURE  GROUP 

OEG  R 

TIH  ) 

TABULAR  VALUE  OF  TEMPERATURE  (lNOEPENOENT 

VARIABLE)  IN  PYROLYSIS  GAS  ENTHALPY  TABLE 

OEG  R 

TT2(l*J) 

TABULAR  VALUE  OF  TEMPERATURE  (lNOEPENOENT 

VARIABLE)  IN  MATERIAL  PROPERTIES  TABLES  FOR 

1  TH  TEMPERATURE  ENTRY  IN  TABLE  FOR  J  TH  MATERIAL 

OEG  R 
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LIST  OK  FORTRAN  VARIABLES 
APPEARING  IN  COMMON  STATEMENTS 

VARIABLE  DESCRIPTION  UNITS 

NAMt 


TIS(I.J)  TABULAR.  VALUE  OF  TEMPERATURE  (INDEPENDENT  VAR-  DEG  R 

I ABLE)  IN  MATERIAL  PROPERTIES  TABLES  FOR  1  TH 
TEMPERATURE  ENTRY  IN  TABLE  FOR  J  TH  INDEX  NUMBER 
FOR  DECOMPOSING  BACK-UPS  (ANALOGOUS  TO  TT2) 

-  TA ( I  * J)  TABULAR  VALUE  OF  PLASTIC  MASS  FRACTION  X»  LB  VIRGIN/ 

INDEPENDENT  VARIABLE  IN  TABLES  OF  Fl(I,J),  F2(I»J)  LB  MATERIAL 

VFZ  INPUT  OPT1UN  1  VIEW  FACTOR 

VR  TABULAR  INTERPOLATION  RATIO  RETURNED  BY  LOOK 

SUBROUTINE 

WT  NUMBER  OF  THE  F-F UNCTION  TArLE  (NUMBERED  IN  ORDER 

INPUT)  ASSIGNED  TO  THE  MAIN  CHARRING  MATERIAL 

X(  )  MASS  FRACTION  OF  VIRGIN  PLASTIC  PER  UNIT  MASS 

OF  MATERIAL 
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variable 

name 


Ad) 

ASU 

B(I) 

B) 

BF 

RPRM 

BPHHG 

BK 

C(I) 

C« 

CHZ 

CHO 

CMOL 

CMFl 

CMMA 

CHHI 

CNC(  ) 


LIST  OF  FOR JR AN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 

DESCRIPTION 


COEFFICIENT  ON  NEW  TEMPERATURE  OF  NODE  1-1 
IN  EQUATION  OF  NODE  I  IN  IMPLICIT  MATRIX 

CURRENT  CROSS  SECTION  AREA  OF  SURFACE  NODE 

COEFFICIENT  ON  NEW  TEMPERATURE  OF  NOOE  1  IN 
EOUATION  OF  NOOE  I  IN  IMPLICIT  MATRIX 

AN  UNCORRECTED  B  PRIME.  EOUaL  TO  (MOOTS  ♦  MDOTC)/ 
RHOE-UE-CM.  ALSO  SEE  NEXT  ENTRY 

REPLACES  BBB(I.J)  FOR  CONVENIENCE.  ALSO  SEE 
ABOVE  ENTRY 

BPRIME  «  (HOOTS  ♦  MUUTC) /WHoE-UE-CM 
BPR1MEG  *  MDOTG/RHOE-UE-CM 
CH/CHO 

COEFFICIENT  ON  NEW  TEMPERATURE  OF  NOOE  1*1  IN 
IN  EQUATION  OF  NOOE  I  IN  IMPLICIT  MATRIX 

HEAT  TRANSFER  COEFFICIENT  -  RMOE-UE-CH 

HEAT  TRANSFER  COEFFICIENT  NOT  CORRECTED  FOR 
BLOWING  ' 

CHAR  RATE*  MOOTC 

LN  (MDOTt/RHOE-UE-CHI 

VALUE  OF  THAT  PART  OF  EROSION  RATE  DUE  TO  FAILING 

MINIMUM  (FROM  VARIOUS  TABLES)  VALUE  OF  AVERAGES 
OF  NEXT  TWO  SUCCESSIVE  TABULAR  VALUES  OF  TLMC 
(  .  *  )  ABOVE  CURRENT  CMOL.  USED  TO  LIMIT 
CORRECTION  SIZE  ON  CMOL 

MAXIMUM  VALUE  (FROM  VARIOUS  TABLES)  OF  AVERAGES 
OF  NEXT  TWO  SUCCESSIVE  TABULAR  VALUES  OF  TLMC 
(  *  .  )  BELOW  CMOL.  USED  TO  LIMIT  CORRECTION 
SIZE  ON  CmOL 

VALUE  OF  CHAR  THERMAL  CONDUCTIVITY  AT  A  NODE 
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•  .  .'i 


UNITS. 


BTU/FT3  - 
OEG  F 

FT2 

BTU/FT3  - 
DEG  F 

SEC- 1 


BTU/FT3  - 
DEG  F 

LB/FT 2-SEC 
LB/FT2-SEC 

LB/FT2-SEC 

*  ) 

LB/FT2  SEC 


BTU/FT-SEC 
OEG  F 


r>  * 

LIST  OF  FORTRAN  VARIABLES 

APPEARING  IN  SUBROUTINE  CBN 

V AMI ABLE 
name 

U£ SCRIPT  ION 

UNITS 

cnoii) 

THERMAL  CONDUCTIVITY  OF  NOOF  I  GENERATED  FOR 

OUTPUT  PUhHOSES  IF  NCON  DIFFERS  FROM  ZERO 

BTU/FT  SEC 
OEG  R 

CHT 

TOTAL  AMOUNT  OF  CHAR  ABLATED  ABAY.  PER  UNIT 

AREA  OF  ORIGINAL  SURFACE 

LB/FT2 

CN(  ) 

VALUE  OF  IHERHAL  CONDUCTIVITY  AT  A  NODE 

BTU/FT-SEC- 
OEG  F 

COLO 

SAVED  PREVIOUS  VALUE  OF  CHAR  INTERFACE  DEPTH 

IN 

CH(  ) 

SPECIFIC  hEAT  AT  A  NUDE 

BTU/LB-OEG  F 

CPC(Ii 

SAVED  VALUE  OF  CHAR  SPECIFIC  HEAT  FOR  NODE  I 

BTU/LB  CHAR* 
OEG  R 

CPE(  ) 

char  or  pyrolysis  interface  depth 

IN 

CPGAS 

SPECIFIC  heat  of  PYROLYSIS  GAS 

BTU/LB-OEG  F 

CPNL 

specific  heat  at  last  node  of  ablating  material 

btu/lb-deg  f 

CHV  (  > 

SPECIFIC  hEAI  OF  VIRGIN  PLASTIC  COMPONENT  AT  A 

NODE 

btu/lb-oeg  F 

CHI 

SPECIFIC  hEAT  OF  TOP  NOuELET  IN  CURRENT  NODE 

Bit  'LB-DEG  F 

C4 

RELATIVE  THICKNESS  PARAMETER  USED  IN  DRDPPINb 

LAST  NODE 

0(1) 

RIGHT  HAND  TERM  IN  IMPLICIT  MATRIX  EQUATION 

FOR  NODE  1 

BIU/FT3 

OLOT 

ABSOLUTE  VELOCITY  OF  CHAR  INTERFACE 

IN/SEC 

Ok  COM 

CURRENT  ACCUMULATED  (NODE  BY  NODE)  AMOUNT  OF 

ENERGY  ABSORPTION  IN  DECOMPOSITION 

BTU/FT2-SEC 

OECOHT 

TIME  INTEGRATED  VALUE  Of  OECOM*  AOJUSTEO  TO 

UNIT  AREA  OF  ORIGINAL  SURFACE 

j 

BTU/FT2 

OEOT 

CURRENT  ACCUMULATED  (NODE  BY  NODE)  RATE  OF  ENERGY 
ABSORPTION  IN  SOLlD^t  (BACK-UPS  EXCLUDED) 

BTU/FT2-SEC 

OLOTT 

TIME  INTEGRATED  VALUE  OF  OEOT*  AOJUSTEO  TO  UNIT 

area  of  original  surface 

BTU/FT2 

DtLCrf 

MINIMUM  OF  DtL(l)/FJFH  ANO  DELM/FJFS 

FT 
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LIST  Of  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 


VARIABLE 

NAME 

DESCRIPTION 

0 

UNITS 

OfcLR 

thickness  ratio  used  in  dropping  last  node 

OfcNOLD 

SAVED  VALUE  OF  DENSITY  OF  PREVIOUS  NOOELET  USED, 

IN  INTERPOLATING  CHAR  AND  PYROLYSIS  INTERFACES 

LB/FT3 

OfcP(I.J) 

J  TH  DEPTH  OF  THE  I  TH  ISOTHERM  IN  THE  BACKUP 
MATERIAL 

IN 

OERR 

(11  IN  ABLATING  SECTION  OF  SURFACE  ENERGY  BALANCE 

BTU/FT2-SEC 

OERR 

RATE  OF  CHANGE  OF  ERR  KITH  lN  8  PRIME  CHAR  UN 
ABLATING  SECTION  OF  SURFACE  ENERGY  BALANCE  CALCU¬ 
LATIONS)  OR  WITH  SUkFACE  TEMPERATURE  (IN  NON¬ 
ABLATING  SECTION) 

BTU/FT2-SEC 

OR 

BTU/FT2-SEC- 
DEG  F 

OIOT 

RATE  OF  SURFACE  RtCtSSION 

IN/SEC 

OM06 (  ) 

ACCUMULATED  AMOUNT  Of  PYROLYSIS  GAS  GENERATION 

IN  A  NODE  PER  UNIT  AnEA  AND  TINES  NUMBER  OF 

NCDELETS  PER  NOOfct  FINALLY  AOJUSTEO  TO  AMOUNT  , 

OF  PYROLYSIS  GAS  GENERATION  IN  A  NODE 

LB/SEC-TT2 

SURFACE 

On  IV 

DERIVATIVE  OF  EMISSlVlTV  BlTH  RESPECT  TO 

TEMPERATUkE 

OEG  F— 1 

DNCPl  ) 

CHAR  ANO  PYROLYSIS  LINES  CRfTERiAL  DENSITIES* 

SEE  CHCRI  ANO  PYCRI 

L8/FT3 

ONS 

CURRENT  NOOELET  DENSITY 

LB/FT3 

OPOT 

ABSOLUTE  VELOCITY  OF  PYROLYSIS  INTERFACE 

IN/SEC 

ORL 

OENSITY-THICKNESS  PROOUCT  FOR  NOOE  ABOVE  NOOE  TO 

BE  DROPPED 

LB/FT2 

OMLC 

OENSITY-THICKNESS— SPECIFIC  HEAT  PROOUCT  FOR  NOOE 
ABOVE  NOOE  TO  BE  DROPPED. 

8TU/FT2- 
OEG  F 

ORL  CP 

I 

OENSI TV-THlCKNESS-SPtCIF IC  HEAT  PROOUCT  FOR  NODE 

TO  BE  DROPPED 

BTU/FT2  - 
OEG  F 

ORLP 

OENSITY-THICKNESS  PROOUCT  FOR  NOOE  10  BE  DROPPED 

LB/FT2 

OHOAC 

RATE  OF  CHANGE  Of  OtNSITf  Of  COMPONENT  A  AT 

CONSTANT  A  DOE  TO  CONVECTION 

LB/FT3  - 

RESIN-SEC 

• 

DhCAT 

RATE  OF  CHANGE  OF  OENSI TY  Of  COMPONENT  A  AT 

CQNSTANT  Y 

LS/FT3  - 
RESIN-SEC 

a 


■/ 


■ 


/ 


VARIABLE 

NAME 


LIST  Uf  FORTRAN  variables 
APPEARING  IN  SUBROUTINE  CBM 

INSCRIPTION 


OKOBC 

OROBT 

DROCC 


RATE  OF  CHANGE  OF  DENSITY  OF  COMPONENT  B  AT 
CONSTANT  X  DUE  TO  CONVECTION 

RATE  OF  CHANGE  OF  OENSITY  OF  COMPONENT  B  AT  / 
CONSTANT  Y  .0 

RATE  OF  CHANGE  OF  OENSITY  OF  COMPONENT  C  AT 
CONSTANT  X  DUE  TO  CONVECTION 


DNOCT 


RATE  OF  CHANGE  of  density  of  component  c  at 
constant  V 


DKOOTD  TOTAL  RATE  OF  CHANGE  OF  DENSITY  AT  CONSTANT  Y  FOR 

CURRENT  NU0E.  -OMOG(l)/(RR(T)*OEL(l> > 

DROT(J)  LOCAL  HATE  OF  CHANGE  1*  DENSITY  OF  COMPONENT  J*  „ 

WHERE  J  IS  THE  COMPONENT  INDEX  (J*l»2.3>  IN  THE 
DECOMPOSING  BACK-UPS  (ANALOGOUS  TO  DRUAT.  OROBT. 
DROCT) 

DS  SURFACE  RECESSION  OOH1NG  TImE  STEP 

OSDT  SURFACE  RECESSION  RATE 

DSDTB  NEW  VALUE  OF  SURFACE  RECESSION  RATE 

OSI  SURFACE  RECESSION  DURING  TlM6  STEP 

DsS  PARAMETER  FJF/OEL(I)*OS/RR(!> 

DIA  TEMPERATURE  INCREMENT  PER  HALF  NOOELET  FOR 

LINEAR  INTERPOLATION  IN  SPACE 

DTH  TIME  INCREMENT 

OTHC  INPUT  maximum  allowable  tihf  STEP 

•I 

DIHS  SAVEO  PREVIOUS  VALUt*OFvBTM 

DIS  CMAf*6£  .IN  SURFACE  UMPERATUrE  DURING  PREVIOUS 

COMPUTATION 

OVB  CONDUCTIVITY  PARAMETER 


OVBS 


SAVEO  VALUE  OF  DVB  FOR  OOMP 


UNITS 


LB/FT3  - 
RESIN-SEC 

L.B/FT3  •- 
RESIN-SEC  . 

LB/FT 3 
REINFORCE¬ 
MENT-SEC 

LB/FT 3 
REINFORCE¬ 
MENT-SEC 

.  ,  LB/FT  3-SEC 


LB/FT 3-SEC 


FT 

FT/SEC 

*> 

Ff/SEC 

IN 


OEG  F 


SEC 
SEC 
SEC 
OEG  F 


BTU/SEC- 
DE6  F-F 12- 
SURFACE 

BTU/SCC  OCR  r 
FT 2- SURF ACC 
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LIST  Ufr  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBN 


VARIABLE 

NANt 

DESCRIPTION 

UNITS 

01 

VALUE  USCO  IN  0R0PPIN6  LAST  NOOC 

on  i 

OUHHY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 
DERIVATIVES  OF  OCPENOCNT  VARIABLES 

VARIOUS 

D2I  ) 

DUMMY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 
DERIVATIVES  OF  OEPEMUENT  VARIABLES 

VARIOUS 

OJ<  ) 

DUMMY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 
DERIVATIVES  OF  DEPENDENT  VARIABLES 

VARIOUS 

e 

REPLACES  ££(I * J)  FOR  CONVENIENCE 

OEG  R 

EoO 

ENERGY  LEAVING  SURFACE  WITH  PYROLYSIS  GASES* 

GS*MGAS  AT  SURFACE  TEMPERATURE 

BTU/FTZ-SEC 

EITER4  1 

VALUE  OF  ERROR  TERM  IN  SURFACE  ENERGY  BALANCE* 

SAVtD  FOR  POSSIBLE  DUMP  FOR  ALL  ITERATIONS 

BTU/FT 2-SEC 

EH|V 

ABLATING  SURFACE  EHlSSlVlTV 

V 

EHO(  t 

GENERAL  VARIABLE  NAME  FOR  CuRVE  SLOPE 

VARIOUS 

EHFX 

COMBlN^rrWn*  TERMS  IN  SURFACE  ENERGY  BALANCE 
EQUATIONS  WHICH  00  NOT  CHANGE  DURING  ITERATIONS 

BTU/FTZ-SEC 

E«R 

ERROR  (DEPARTURE  FROM  ZERO)  IN  SURFACE  ENERGY 
BALANCE 

BTU/FTZ-SEC 

EKRC 

CORRECTION  APPLlEO  TO  LN(MOOTC/KHO£-UE-CM>  OR  TO 
SURFACE  ENERGY  BALANCE  OPERATION 

-  OR  OEG  R 

ERRS 

SAVED  VALUE  OF  ERROR*  ERR 

BTU/FTZ-SEC 

'It 

FACTOR  USED  IN  DROPPING  LAST  NOOE 

F 

REPLACES  FF(I.J)  FUR  CONVENIENCE 

ILB/FT3)** 
(I.-PSI (1*3)) 

FA 

DEFINED  As  (1.-PSIA)«BA*(RH00A**|1.-PSIAI) 

ILB/FT3I** 

(1.-PSIAI/SCC 

FACT1 

FACTOR  DTh/(DELII)*RH(I>) 

SCC/FT 

fact? 

FACTOR  GSM/U£L(1>*RH(1) 

LB/FT3-SEC 

Fb 

OEFINEO  AS  4i.-PSIB)«B8»tRH00B«*ll.-PSIB)> 

ILB/FT3I** 

(1»-PSIB)/SEC 
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LIST  Of  FORTRAN  VARIABLES 
appearing  in  subroutine  cbm 


variable 

name 


description 


FC  DEFINED  AS  (I«-PSIC)*BC*(RHrtOC**(l*-PSIC)  ) 


FJf  NUMBER  OF  NOOELETS  IN  CURRENT  NODE 

Fa  FACTOR  USED  IN  COMPUTING  CONVECTIVE  CONTRIBUTIONS 

TO  OENSITt  CHANGES  IN  NOOELETS  OF  LAST  NOOC  OF 
ABLATING  MATERIAL 

FA  VALUE  USED  IN  DROPPING  LAST  NOOE 

n  vXluc  USED  IN  DROPPING  LAST  NOOE 

GAN  REPLACES  GA(I)  FOR  CONVENIENCE 

GS  TOTAL  PYROLYSIS  GAS  FLO*  OUT /4 ME  ABLATING  SURFACE 

6SEGR  ENERGY  TERM*  SUN  OF  HGAS*ON06(  I  OVER  NOOCS 


6SM  ACCO^JLATtD  PYROLYSIS  GAS  FLOW  RATE  ENTERING  A 

NOOE*  ON  A  UNI T-AREA-OF -SURF ACE  8ASIS 

GSNS  PYROLYSIS  GAS  FLOW  OUT  THE  ABLATING  SURFACE 


GSMT  TOTAL  (TIME  INTEGRATED)  PYROLYSIS  GAS  FLOM  OUT 

THE  ABLATING  SURFACE  ADJUSTED  TO  UNIT  AREA  OF 
ORIGINAL  SURFACE 

GSM2  ACCUMULATED  PYROLYSIS  GAS  FlOV  RATE  ENTERING  A 

NOOE*  ON  A  UNIT-ARCA-OF-SURFACE  BASIS*  WHICH 
DERIVES  FROM  CHARRING  BACK-UPS 

GSM2S  FLOW  RATE  OF  PYROLYSIS  GAS  ORIGINATING  IN  DECOM¬ 

POSING  BACK-UPS  OUT  THE  ABLATING  SURFACE 

G5H2T  TOTAL  (TIME  INTEGRATED)  FLOW  OF  PYROLYSIS  GAS 

ORIGINATING  IN  DECOMPOSING  BACK-UPS  AOJUSTED  TO 
UNIT  AREA  OF  ORIGINAL  SURFACE 

Gt  FACTOR  USEO  IN  DROPPING  LAST  NODE 

HAPM6  COMBINED  ENTHALPY  PER  UNIT  SURFACE  AREA  OF  NOOE 

TO  BE  DROPPED  AND  ADJACENT  NODE 

HnAR  TEMPERATURE  DEPENDENT  REACTION  ENTHALPY  DEFINED 

AS  (VIRGIm  DENSITY  •  VIRGIN  ENTHALPY  -  CHAR 
DENSITY  •  CHAR  ENTHALPY) /(VIRGIN  OENSITY-  CHAR 
DENSITY) 


UNITS 


(LB/FT3)»« 

(l.-PSIO/SEC 


FT 3  RESiN 

LB /FT 2  SEC 

/ 

BTU/SEC-FT2 

SURFACE 

LB/SEC-FT2 

SURFACE 

LB/SEC-FT2 

SURFACE 

LB/FT2 


LB/SEC-FT2 

SURFACE 


LB/SEC-FT2 

SURFACE 

LB/FT2 


BTU/FT2 


BTU/LB 


VARIABLE 

NAME 


LIST  Of  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 

UESCR1PTION 


UNITS 


HBARS 


HC(  ) 


SAVED  VALUE  OF  HBAR  FOR  DUMP 
ENTHALPY  OF  CHAR  AT  A  NOOE 


BTU/LB 


BTU/LB 


TOTAL  ENTHALPY  OF  EuGE  GASES  FOR  OPTION  I  (EQUALS 
SURFACE  TLHPERATURE  FOR  OPTION  2  OR  VIEW  FACTOR 
FOR  OPTION  3) 

ENTHALPY  OF  PYROLYSIS  GAS 

enthalpy  OF  PLASTIC  AT  A  NOOE 

TOTAL  HEAT  TRANSFER  COEFFICIENT  AT  BACK  WALL 
(INCLUDING  RADIATION) 

ENTHALPY  OF  EDtoE  GASES  AT  WALL  TEMPERATURE 
(OPTION  I)  OR  ENTHALPY  OF  PYROLYSIS  GASES  AT 
WALL  TEMPERATURE  (OPTION  3) 

ENTHALPY  OF  TOP  NOOLLET  IN  A  NOOE 

COUNTING  INOEX.  USUALLY  FOR  NOOES 

FLAG  USED  TO  OETECT  FIRST  PASS  THROUGH  ABLATING 
SURFACE  ENERGY  balance  package 

COUNTING  INDEX 

INOEX  FOR  GAS  RATES  (MOOTG/RHOE-UE-CH) 

INTEGER  ZERO 

INDEX  FOR  PRESSURE  tNTRIES 
SOMETIMES  REPLACES  1R(  I  FOB  EFFICIENCY 
SOMETIMES  REPLACES  IR(  )  FOR  EFFICIENCY 
SOMETIMES  REPLACES  IR(  )  FOR  EFFICIENCY 
SOMETIMES  REPLACES  1R(  )  FOR  EFFICIENCY 
NOT  CURRENTLY  USED 

SAVLP  MATERIAL  IDENTIFICATION  NUMBER  FOR  FIRST 
NOOE  OF  BACK-UP  MATtNlAL 

NUMBER  OF  PASSES  THROUGH  MAIN  COMPUTATION  LOOP 


BTU/LB* 

DEG  R»  OR  -• 


BTU/LB 

BTU/LB 

BTU/FT2-SEC- 
DEG  F 

BTU/LB 


BTU/LB 
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LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBN 

VARIABLE  DESCRIPTION  UNITS 

NAME 


ITL  LIMIT  NUMBERS  (ASSUMES  2  VALUES)  ON  ITS  TO  CON¬ 

TROL  ITERATION  EVENTS  IN  ABLATING  SURFACE  ENERGY 
BALANCE 

ITS  COUNTER  FOR  NUMBER  OF  ITERATIONS  ON  SURFACE 

ENERGY  BALANCE 

11  REPLACES  SUBSCRIPTED  INOEX  FOR  EFFICIENCY 

12  REPLACES  SUBSCRIPTED  INDEX  FOR  EFFICIENCY 

U  REPLACES  SUBSCRIPTED  INOEX  FOR  EFFICIENCY 

I<»  REPLACES  SUBSCRIPTED  INOEX  FOR  EFFICIENCY 

J  GENERAL  INDEX*  NOOELET  COUNTER  IN  A  GIVEN  NODE 

IN  DECOMPOSITION  BLOCK 

J1  INDICATES  NUMBER  OF  FIRST  NODELET  IN  A  NODE* 

COUNTED  NOOE  BY  NODE  (EQUALS  JFH  ♦  1  FOR  FIRST 
NODE*  1  FOR  OTHERS) 

K  GENERAL  UTILITY  INOEX 

K1  TEMPORARY  INOEX  IDENTIFYING  AN  ISOTHERM  WITHIN 

A  BACKUP  MATERIAL 

KK  COUNTER  ON  THERMOCOUPLE  ANO  ISOTHERM  OUTPUT 

KKSV  RETURNED  ARGUMENT  IN  CALL  OF  SStftCH  (N*KKN) 

KSCT  SCRATCH  TAPE  LOGICAL  UNIT  NUMBER  VARIABLE 

KT  INDEX  REPLACING  SUBSCRIPTED  MATERIAL  NUMBERS 

MATL (  )  FOR  EFFICIENCY 

L  UTILITY  INDEX 

LL  UTILITY  INDEX*  USUALLY  LOWER  LIMIT  OF  DO  LOOP 

INDEX 

Lu  UTILITY  UPPER  LIMIT  OF  UO  LOOP  INOEX 

M  UTILITY  INOEX 

N  GENERAL  UTILITY  INOEX*  USEO  FOR  COUNTING  TOTAL 

NUMBER  OF  NODELETS  IN  DECOMPOSITION  BLOCK 

NISO(KI)  NUMBER  OF  DEPTHS  AT  WHICH  THE  KI  TH  ISOTHERM 
OCCURS  AT  CURRENT  OUTPUT  TIME 
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variable 

name 

NOR 

NLl 

NLH 

NZ 

0 

OMGAM 

OVR 

PGPU 

PGPUf 

PHI 

POLO 

POW 

PHES 

olmem 

OCHEMT 

OCONO 

QCOMOT 

QCONV 

QCONV T 


LIST  Of  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 

DESCRIPTION 

NUMBER  Of  NODES  WHICH  HAVE  BEEN  DROPPED 
NUMBER  Of  LINES  Of  NOOAL  DATA  OUTPUT 
NL-1 

CURRENT  TOTAL  NUMBER  OF  NOOELETS 

FILLER  FOR  UNNEEOEO  VARIABLE  NAMES  IN  LOOK  CALL 

REPLACES  OMGAII)  FOR  CONVENIENCE 

MEASURES  APPROACH  Of  DENSITY  TO  RESIDUAL  DENSITY 
EGO  -  GSEGR 

TINE  INTEGRATED  VALUE  OF  PGPU*  ADJUSTED  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

PARAMETER  2.*8RP*DF 

SAVEO  PREVIOUS  VALUE  OF  PYROLYSIS  INTERFACE  DEPTH 

EXPONENT  IN  PYROLYSIS  DECOMPOSITION  RATE 
CALCULATIONS 

CURRENT  VALUE  OF  LOG  OF  PRESSURE 

FROM  TABLE  LOOK-UP.  IS  INITIALLY  TERM 
(CM/CH)*<tHEM  PROO)  •  HR), 

LATER  ADJUSTED  TO  RHOE-UC-Cm* ICHEM  PROO!  -OCMEM 

TIME  INTEGRATED  VALUE  OF  OCHCM  AOJUSTEO  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

OCONO.  ENERGY  CONDUCTED  INTO  CHARRING  HA TER l At  AT 
HEATED  SURFACE 

TIME  INTEGRATED  VALUE  OF  OCONU  AOJUSTEO  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

INITIALLY  OUNMY  NAME  FOR  HER.  AOJUSTEO  TO 
OCONV  ■  RM0E-UE-CM«|HE-HW» 

TIME  INTEGRATED  VALUE  OF  QCOMV  AOJUSTEO  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 
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*  * 


UNITS 


FT3  REIN¬ 
FORCEMENT/ 
FT3 

LB/FT3 

BTU/SCC-FT2- 

SURFACE 

8TU/FT2 


8TU/L0  OR 
BTU/FT2-SCC 

BTU/FT2 

BTU/FT2-SCC 

BTU/FT2 

BTU/LB* 

BTU/FT2-SEC 

BTU/FT2 


variable 

name 

QLOSS 

QLOSST 

00 

QHA 

QRP 

QRPT 

RAO 

RAOT 

RAT  (  ) 

RO 

RkX 

RU(  ) 

RON  (  ) 

ROOZ 

RUTC  » 

ROI 
RM(  ) 


LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBH 

DESCRIPTION 


ENERGY  LOSS  QLOSS  AT  REAR  FACE  OF  ABLATING  MATER¬ 
IAL*  PER  UNIT  AREA  OF  FRONT  SURFACE 

TIME  INTEGRATED  VALUE  OF  QLOSS  AOJUSTEO  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

OUMMY  NAME*  INITIAL  QCONV  FOR  INTERPOLATION 

RAO I ANT  HEAT  FLUX  TO  SURFACE  ORAD  IN  OPTIONS  1 
ANO  3  (IN  OPTION  2*  SURFACE  RECESSION  RATE) 

RAOIANT  HEAT  ABSORBED  AT  SURFACE*  EQUAL  TO 
EMISSIVITY  •  ORA 

TIME  INTEGRATED  VALUE  OF  ORP*  ADJUSTED  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

RATE  OF  ENERGY  RA01ATI0N  AWAY  FROM  SURFACE  < 

SIGMA  •  EMISSIVITY  •  VIEW  FACTOR  •  TS**4 

TIME  INTEGRATED  VALUE  OF  RAD*  ADJUSTED  TO  UNIT 
AREA  OF  ORIGINAL  SURFACE 

CONDUCTION  RESISTANCE  OEL(  )/(CN(  )*RR(  )) 


DENSITY  DIFFERENCE  USED  IN  PYROLYSIS  CALCULATIONS 

BURNING  RATE  EXPONENT  PARAMETER  ( 1.6-.2*BREX)/ 
(l.-BREX) 

DENSITY  OF  A  NODE  (USUALLY  OLD  DENSITY) 

OENSITY  OF  A  NODE  (USUALLY  NEW  DENSITY) 

DENSITY  PARAMETER*  EQUAL  TO  ZERO  OR  TO  FJFH  TIMES 
UNIFORM  NOOCLET  DENSITIES  IN  AREAS  WITHOUT 
DECOMPOSITION  REACTIONS  (PURE  PLASTIC  OR  PURE 
CHAR)*  USED  TO  EXPEU1TE  DENSITY  CALCULATIONS 

DENSITY  OF  TOP  NOOELET  IN  A  GIVEN  NODE 

REPLACES  ROT (  )  FOR  EFFICIENCY 

OBTAINED  FROM  OGLE  AS  CROSS-SECTION  AREA  OF  A 
NOOC*  LATER  NORMALIZED  ON  CROSS-SECTION  AREA  OF 
SURFACE  NODE 
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UNITS 

8TU/SEC-FT2- 

SURFACE 

BtU/FT2 

BTU/LB 

BTU/FT2-SEC 
OR  MILS/SEC 

BTU/FT2-SEC 

BTU/FT2 

BTU/FT2-SEC 

BTU/FT2 

(FT2  SURF- 
SEC-OEG  F)/ 
BTU 

LB/FT3 

LB/FT 3 
LB/FT3 
LB/FT3 

LB/FT3 
L8/FT3 
FT2  OR  — 


*-  , 


LIST  OF  FORTRAN  VARIABLES 

APPEARING  IN  SUBROUTINE  CBM 

variable 

NANt 

DESCRIPTION 

UNITS 

RSU 

CURRENT  SURFACE  RA01US  (TOTAL  RECESSION"  FOR 

PLANE  GEOMETRY 

IN 

SA 

CURRENT  VALUE  OF  TOTAL  SURFACE  RECESSION 

IN 

SONET 

RATE  OF  RECESSION  AFTER  RATF  OF  CHAR  SWELL  IS 
ACCOUNTED  FOR 

IN/SEC 

SIG 

STEFAN-BOLTZMANN  CONSTANT 

BTU/FT2-SEC- 
OEG  R**4 

SNET 

NEW  SURFACE  RECESSION  AFTER  SWELL  IS  ACCOUNTED  FOR 

IN 

SOEGR 

TERM  HGAS«DM0G(  )»  SUMMED  OVER  ALL  NODES 

BTU/FT2  - 
SURFACE-SEC 

T 

TERM  USED  IN  EVALUATING  SOLID  CONVECTION  ENERGY 

BTU/FT3 

TABC 

INTERPOLATED  TEMPERATURE  APPROXIMATING  BOTTOM  OF 
ABLATION-SECTION  OF  SURFACE  THERMOCHEMISTRY 

TABLES*  COMPARED  TO  OLD  VALUE  OF  SURFACE  TEMPER¬ 
ATURE  TO  JUDGE  IF  ABLATION-SECTION  SHOULD  BE  USED 

OEG  R 

TAS 

NOOELET  TEMPERATURE 

DEG  R 

TB 

SUMMATION  OF  TN*DSDT*RR ( I )  OVER  NODES  OF  ABLATING 
MATERIAL 

BTU/SEC-FT2- 

SURFACE 

TEMP 

TEMPORARY  NOOAL  TEMPERATURE 

DEG  R 

TERH1 

FACT  I  *  OVB 

# 

BTU/FT-DEC  F 
FT2-SURFACE 

TERM2 

LUMPING  OF  TERMS  USEFUL  IN  COMPUTING  IMPLICIT 
TEMPERATURE  COEFFICIENTS*  SFE  STATEMENT  14  ♦  6 

BTU/FT2  - 
OEG  F 

TERM3 

LUMPING  OF  TERMS  USEFUL  IN  COMPUTINO  IMPLICIT 
TEMPERATURE  COEFFICIENTS*  SEE  STATEMENT  25 

BTU/FT3- 
DEG  F-SEC 

TM 

TIME 

SEC 

TMOS 

VALUE  OF  TIME  AT  A  TIME-TABLE  DOUBLE  ENTRY 
(SHIFT  OF  OPTION) 

SEC 

THPRT 

OUTPUT  TIME 

SEC 

TN 

TERM  USED  IN  EVALUATING  SOLID  CONVECTION  ENERGY 

BTU/FT3 

tim  > 

THERMOCOUPLE  TEMPERATURE  OR  ISOTHERM  DEPTH 

DEG  R  OR  IN 
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LIST  UF  FOR TRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 


variable 

name 

DESCRIPTION 

UNITS 

TOPI 

COMBINED  HASS  OF  NODE  TO  BE  OROPPED  AND  ADJACENT 

NODE 

LB/FT2 

TOP2 

COMBINED  THERMAL  CAPACITY  OF  NOOE  TO  BE  DROPPED 

AND  ADJACENT  NODE 

BTU/FT2 
DE6  R 

TUP3 

COMBINED  ENTHALPY  OF  NODE  To  BE  OROPPED  AND 

ADJACENT  NODE 

BTU/FT2 

TKACMI 

MINIMUM  VALUE  OF  A  SET  OF  TRACd.J)*  J*  1*3 

DEG  R 

TS 

SURFACE  TEMPERATURE*  EOUlVALENCED  TO  TA(1) 

OEG  R 

TSAVE 

SAVED  VALUE  OF  SURFACE  TEMPERATURE  FROM  PREVIOUS 

TIME  STEP 

DEG  R 

TSMA 

MINIMUM  (FROM  VARIOUS  TABLES)  VALUE  OF  AVERAGES  OF 
NEXT  TWO  TABULAR  VALUES  OF  TTSC  *  *  > ABOVE  CURRENT 
TS.  USED  TO  LIMIT  CORRECTION  SIZE  ON  TS 

DEG  R 

TSMI 

MAXIMUM  (FROM  VARIOUS  TABLES)  VALUE  OF  AVERAGES  OF 
NEXT  TWO  TABULAR  VALUES  OF  TTS  (  *  *  >  BELOW 

CURRENT  TS*  USED  TO  LIMIT  CORRECTION  SIZE  ON  TS 

DEG  R 

TSSO 

SQUARE  OF  SURFACE  TEMPERATURE 

DEG  R«»2 

TT 

TIME  INTEGRATED  VALUE  OF  T6  ADJUSTED  TO  UNIT 

AREA  OF  OHIG1NAL  SURFACE 

BTU/FT2 

Vf 

VIEW  FACTOR 

VITER (  > 

VALUE  OF  LN(MDOTC/RHOE-UE-Cm)  OR  SURFACE  TEMP¬ 
ERATURE  (FOR  A  NON-ABLATING  CASE)  SAVED  FOR  EACH 
ITERATION  OF  SURFACE  ENERGY  BALANCE*  FOR 

POSSIBLE  OUMP 

—  OR 
DEG  R 

VOL 

VOLUME  OF  NEW  AMALGAMATED  LAST  NODE  IN 

MAIN  MATERIAL 

FT3/FT2- 

SURFACE 

VRM 

REPLACES  VR  FOR  PYRuLYSIS  GAS  RATE  INTERPOLATION 

VRP 

REPLACES  VR  FOR  PRESSURE  INTERPOLATION 

XP1 

REPLACES  Ad)  FOR  EFFICIENCY 

XI 

X<  )  FOR  TOP  NOOELET  IN  A  NODE 

YU  ) 

DUMMY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 

DEPENDENT  VARIABLES 

VARIOUS 
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VARIABLE 

NAME 

Y2<  ) 

V3(  1 


LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  CBM 

OESCRIPTION 


DUMMY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 
DEPENDENT  VARIABLES 

OUMMY  VARIABLE  NAME  USED  IN  LOOK  CALLS  FOR 
DEPENDENT  VARIABLES 


UNITS 

VARIOUS 

VARIOUS 
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LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  INPOUT 

VARIABLE  DESCRIPTION  UNITS 

NAHE 

At  EXPONENT  ON  RAOIUS  GIVING  NoDAL  CROSS-SECTIONAL 

AREA  VARIATION  WITH  RAOIUS 

A V  ALPHAMERIC  IDENTIFICATION  VARIABLE  FOR 

COMPONENT  A 

B  REPLACES  BLANK 

BP  8  PRIME  *  (MDOTC*MDOTG)/RMOf-UE-CM 

BPG  BPRIMEG  >  MOOTG/RHOE-UE-CM 

B9  ALPHAMERIC  IDENTIFICATION  VARIABLE  FOR 

COMPONENT  B 

INPUT  VALUE  OF  CMH 

UNUSED  PLACE  HOLOING  VARlABl  E  IN  LOOK  CALL 
UNUSEO  PLACE  HOLOING  VARlABl E  IN  LOOK  CALL 
ALPHAMERIC  IDENTIFICATION  VARIABLE  FOR  COMPONENT  C 
TEST  VALUE  OF  PYROLYSIS  GAS  RATE*  MOOTG/RHOE-UE-CM 
UNUSEO  PLACE  HOLDING  VARlABl  E  IN  LOOK  CALL 

MASS  FRACTION  OF  RESIN  (COMPONENTS  A  AND  B>  IN  LB  RESIN/ 


UNDECOMPOSED  PLASTIC 

LB  VIRGIN 

HCH 

ENTHALPY  OF  CHAR 

BTU/LB 

HE 

TOTAL  ENTHALPY  OF  EOGE  GASES  AT  A  GIVEN 
TEMPERATURE 

BTU/LB 

HGA 

ENTHALPY  OF  PYROLYSIS  GAS 

r 

BTU/LB 

HSH 

DATUM  ENTHALPY  AT  REFERENCE  TEMPERATURE  TZ 

BTU/LB 

HZ 

Z  ENTHALPY  TERM  FOR  EOGE  GAS  AT  A  GIVEN 
TEMPERATURE 

8TU/3B 

I 

UTILITY  INOEX 

IIN 

INDEX  USED  IN  OUTPUTTING  PYROLYSIS  GAS 

enthalpy  table 

\ 

CMHS 

CTI 

CT2 

C9 

QMS 

QUM 

GAMAM 
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I 


VARIABLE 

NAME 


list  or  fortran  variables 

APPEARING  IN  SUBROUTINE  INPOUT 
DESCRIPTION 


UNITS 


IG  FLAG  USED  TO  MARK  A  MAXIMUM  TEMPERATURE  ENTRY 

IN  A  CHAR  RATE  TABLE  , 

IN  UTILITY  INOEX 

I OPT (I)  HEATING  OPTION  (1*2*  OR  3)  FOR  TINE  TABLE  TABULAR 

ENTRY  I 

IP  INOEX  ON  PRESSURE  IN  SURFACE  THERMOCHEMISTRY 

TABLE  INPUT 

IPN  INOEX  ON  PRESSURE  IN  SURFACF  THERMOCHEMISTRY 

TABLE  INPUT 

IS  SAVED  VALUE  OF  II 

JT  UTILITY  INOEX 

IA  FLAG  ON  TYPE  OF  ERROR  STOP  IN  SURFACE  THERMO¬ 

CHEMISTRY  INPUT 

IH  l  OUTPUT  INOEX  ARRAY  FROM  SUBROUTINE  OROERD 

J  UTILITY  INOEX 

JbU  COUNTER  FOR  NUMBER  OF  CHARRfNG  BACK-UPS  ASSIGNED 

TO  AN  F— FUNCTION  TABLE 

JNG  FLAG  TO  IDENTIFY  TYPE  OF  SURFACE  THERMOCHEMISTRY 

TABLE*  -I  FOR  E06C  TABLE*  0  FOR  ZERO  CHAR  RATE 
TABLE  (INDEPENDENT  SURFACE  TEMPERATURE)*  GREATER 
THAN  0  FOR  SURFACE  EQUILIBRIUM  (SURFACE  TEMPER- 
/  ATURE  #OEP£NO£NT) 

J)  INOEX  LIMIT  EQUAL  TO  NUMBER  OF  NOOELETS  IN 

CURRENT  NOOE 

K  UTILITY  INOEX 

KH  ORDINAL  RANK  OF  DECOMPOSING  BACK-UP  MATERIAL* 

COUNTING  FROM  BACK  MALL  OF  MAIN  ABLATING  MATERIAL 

KtiS  SAVEO  VALUE  OF  KH 

si 

KKSW  RETURNED  ARGUMENT  IN  CALL  Of  SSMTCH  (NtKKSW) 

KMTL(I)  INPUT  ARRAY  OF  MATERIAL  SEQUENCE  NUMBERS  (ZERO  FOR 

MAIN  MATERIAL*  1  TO  5  FOR  CHARRING  BACK-UPS) 
ASSIGNED  TO  CURRENT  f -FUNCTION  TABLES 


• ;  \ 

,  t  V  4 


LIST  OF  FORTRAN  VARIABLES 
APPEARING  IN  SUBROUTINE  INPOUT 


j 


VARIABLE 

NAME 

DESCRIPTION 

A  * 

UNITS 

KNST 

* 

CHECK  INTEGER  TO  FLAG  FIRST  SURFACE  TABLE  INPUT 

KSV(J) 

SAVED  NON-ZERO  VALUE  OF  KMT|  (1) 

KT 

INOEX  REPLACING  SUBSCRIPTED  MATERIAL  NUMBERS 

MATLI  >  FOR  EFFICIENCY 

L 

UTILITY  INOEX  '  *  * 

LL 

UTILITY  INOEX 

* 

LLL 

NUMBER  OF  LINES  OF  OUTPUT  Or  A  GIVEN  CHAR  RATE 
TABLE 

LU 

UTILITY  UPPER  LIMIT  ON  DO  LodPS 

N 

INDEX  FOR  NODELETS  « 

NC _ 

FLAG  MARKING  PYROLYSIS  GAS  FNTHALPY  CARD  PAIRS 

AND  tNO"OF  FUNCTIONShOF-TIMF  TABLES 

NOBUCH 

CHECK  VALUE  FOR  NOSU  SUMMED' DURING  'NODAL  DATA 

input  r 

AlAAT 

NUMBER  OF  ENTRIES  IN  CURRENT  MOOTC  TABLE 

NOP 

EQUAL  TO  NO  •  1 

NOPT 

ACCUMULATED  NUMBER  OF  OPTIONS  APPEARING  IN  A 

time-table 

NSEN 

NUMBER  OF  ENTRIES  IN  CURRENT  EDGE  TABLE 

NST 

INPUT  FLAG*  NON-ZERO  VALUE  CALLING  FOR  RE-USE  OF 
PREVIOUSLY  INPUT  SURFACE  THFRMOCHEMISTRV  TABLES 

NTH 

NUMBER  OF  TIME  POINTS  IN  FUNCTIONS-OF-TIME 

TABLES  1 ACCUMULATED)  * 

NT1 

. INDEX  USED  TO  COUNT  ENTRIES  IN  PYROLYSIS  GAS 
ENTHALPY  TABLE 

PSV 

.  SAVEO  TEST  VALUE  OF  PRESSURF  DURING  INPUT  OF 

SURFACE  THERMOCHEMISTRY  TABLES  * 

ATM 

RKR 

TEMPORARILY  REPLACES  RHOC(K)  OR  RHOV(K)  FOR 
CONVENIENCE 

LB/FT 3 

RSVN 

POSITIVE  VALUE  FOR  RSV 

INCHES 
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LIST  or  FOOT* AM  VARIABLES 


UNITS 


MKMW  IN  SUBROUTINE  INPOUT 

VARIABLE  DESCRIPTION 

NAME 


TCZSCNll.J)  DERIVATIVE  Of  FROZEN  EDGE  Ii3  Z-CNTHALPY  MITh  BTU/ 

RESPECT  TO  TEMPERATURE  AT  X  TH  E06E  TAKE  POINT  LB-OEG  R 

IN  TAKE  FOR  J  TN  PRESSURE 

TSEN(l)  RE  AO  IN  AS  ENTHALPY  IM  FROZEN  CORE  TAKES  AND  BALL  SlU/UT 

ENTHALPY  IN  SURFACE  THERMOCHEMISTRY  TABUS  BUT 
LATER  CONVERTED  TO  FROZEN  El)SE  ENTHALPY 

“  >  R 

T SURF (I)  NAME  OF  SURFACE  SPECIES  FOR  I  TM  TAKE  ENTRY 

,TZ  REFERENCE  TEMPERATURE  FOR  KAY  OF  FORMATION  0C6  R 

tzseniivJ)  z -Enthalpy  tern  at  i  th  entry  in  frozen  edge  tabu  btu/lb  o£e  m 

FOR  J  TH  PRESSURE 


,  VKH  floating  KH 

MLO  UNEQUAL  OIFFUSION  EAPQNENT 

N. 

MLS  SAVED  TEST  VALUE  CF  MLO 


NMl 

CM«I» 

OCP(ltJ) 

OIS 

0152 

1 

INNS 

IPlS 

ISO 

ISOT 

J»KiLtN 

N 

N1S0IKI) 

01 

RAISO(l) 

RAT1I I 

RK(I) 

1 

Ti 


list  or  fortran  tmiaous 

APPEARING  IN  SUBROUTINE  THERM 
OE SCRIPT! ON 


VALUE  Of  THERMAL  CONDUCTIVITY  OF  THE  I  TH  NODE 


J  TH  OEPTh  or  THE  I  ih  isotherm  in  the 
BACKUP  MATERIAL 

DISTANCE  BCTmCEN  A  THERMOCOUPLE  LOCATION  AND  AN 
ADJACENT  MORAL  LOCATION  FOR  INTERP&LAT ION  PURPOSES 

INTERPOLATION  QUANTITY  WHICH  ACCOUNTS  FOR*  CONTACT 
RESISTANCE  IN  THE  EVALUATION  OF  THE  TEMPERATURE 
OF  A  THERMOCOUPLE  AT  AN  INTERFACE 

LOOP  INOEA 

INTEGER  USED  IN  SEARCHING  FqR  ISOTHERM  LOCATION 

INTEGER  USED  IN  SEARCHING  FqR  ISOTHERM  LOCATION 

NUMBER  OF  BACK**  NODE  OR  HUMBER  OF  BACKUP 
HOOFS  •  1 

CURRENT  INOEX  IDENTIFYING  THE  ISOTHERM  TEMPERATURE 

TEMPORARY  NOOE  IDENTIFIER  OR  LOOP  INOEX 

NUMBER  OF  DEPTHS  OF  CURRENT  ISOTHERM 

NUMBER  OF  OEPTHS  OF  THE  KI  TH  ISOTHERM 

HEAT  FLUX  BETWEEN  ADJACENT  MOOES  OF  DISSIMILAR 
MATERIALS*  USED  FOR  INTERPOLATION  IN 
ISOTHERM  LOCATING  LOGIC 

DEPTH  OF  IME  I  TH  NUDE  BELOW  THE  CURRENT  SURFACE. 
I»1  IS  THE  LAST  NOOE  OF  MAIN  MATERIAL 

CONDUCTION  RESISTANCE  OETVEFN  THE  I  TM  NOOE  ANO 
1*1  NOOE 

RELATIVE  AREA  OF  THE  I  TH  NoOE  (AS  OEFINEO  IN  CBM) 

TCiPERATURE  Of  CJRRtMT  ISOTHERM 

TEMPERATURE  OF  NOOE  ABOVE  OR  BELOW  CURRENT 
sotherm 


UNITS 

BTU/FT  SEC 
OEG  F 

IN 

FT 

FT 


BTU/FT2-SEC 

IN 

FT2-SEC- 

OEG-F/BTU 

OEG  R 
OEG  R 


LIST  Of  F OUTRAN  VARIABLES 
APPEARING  IN  SJRHOUTINE  THERM 


VARIABLE  description  UNITS 

NAME 

T*  TEMPERATURE  Of  NODE  BELOV  OR  ABOVE  CURRENT  OEG  R 

ISOTHERM 

TATHNUt  TEMPERATURE  Of  THE  1  TM  NOOf*  1*1  IS  THE  LAST  OF.O  R 

NOOE  Of  MAIN  MATERIAL 

TIL  THE  LOVER  TEMPERATURE  AT  AN  iNTERfACE  HAVING  OEG  R 

CONTACT  RESISTANCE 

TlU  THE  UPPER  TEMPERATURE  AT  AN  INTERFACE  HAVING  OEG  R 

CONTACT  RESISTANCE 

■tv  f  , f 

TO(I)  TEMPERATURE  Of  THE  1  Trt  THERMOCOUPLE  OR  DEPTH  OF  OEG  R 

ISOTHERM  IF  ONLY  ONfc  LOCATION  OCCURS 
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it 


/ 


-  ,i . 


OR 


SECTION  3 


6 


FLOW  CHARTS 

i  *' 

i  Computer  generated  flow  charts  were  produced  and  are  given  in  this 
section.  The  flow  charts  show  transfers  as  lines  on  the  right  edge  of  the  fig¬ 
ures  and  DO- loop  blocks  as  lines  on  the  left  edge  of  the  figures. .  The  order  of 
presentation  corresponds  to  that  used  in  Section  2  for  Fortran  variables  name 
definitions. 


I 


t 


A 


3-1 


y- 


/ 


uuuououu 


CBM 


PflBE  NO  1 


.L 


•gTPffifJTTiT  MWl - 

CHARRING  MATERIAL  THERMAL  RE8P0N8E  AM)  ABLATION  PROGRAM  ALLOWING 
FOR  UP  TO  FIVE  DECOMPOSING  BACK-UP  MATERIALS 

SEE  CBM  LISTING  FOR  COMMON  STATEMENTS  ♦•♦♦♦ 

DIMENSION  STATEMENTS  ***** 

EQUIVALENCE  STATEMENTS 

FORMAT  STATEMENTS 


1390 


0NCP<3>-99933. 

KSCT-17 
SI6-. 41  IE-12 

CHAN  AND  PYROLYSIS  ZO*E  CRITERIA-  DENSITIES 
0NCP<1)*RHO<2>*CHCRI*<RHO<1)-RHO<2  >) 

DNCP<2  > -RHO  <  2  ) +PYCRI*  <  RH0<1  > -RHO«  » 

INITIAL  VALUES  FOR  TIME  LOOP 

ITER— l 

DTHC-OTHB 

IAB4) 

ITS-0 

SA-0.0 

SNET-O. 

GSMTMJ.O 

BSMfT-O. 

6SMS=0. 

CMD-O.O 

RSU-ABS<RSV> 

BS-0.0 

CMFL-0.0 

CMT-O. 

OSOTB-O.O 

DSOT-O. 

0SI«O. 

DIOT-O. 

CPE<1)“0. 

CPE<2>=0. 

COLD-0. 

POUVO. 

OCOT-O. 

DPOT-O. 

QCOND-O. 

acoNV-o. 

QCHEM-O. 

RAOm. 

ORfMJ. 

QCONVT-C. 

QRPT-0. 

RAOT-O. 

QCHEMT-O. 

QCONDT-O. 

PGPUT-0. 

OECOMT-O. 

oeorr-o. 

QLOSST-O. 

TT-O. 

TB»0 . 

SCESR-O. 

GSEGR-0. 

PBPU-O. 

OECOM-O. 

E80-C. 

HW-O, 

BR-0. 

KK-0 

TH*»THZRO 

OTA-O.O 

THOS-THZRO-OTHIN 
THPRT-TH 
REWIND  KSCT 

CALL  L COUNT  <-2,LCT.NPB«RECORD<95>> 

WRITE  (KOUT.S42> 

IS*0 

OTH-OTHIN 


3-3 


non 


1  It  HRE9*S I8*ePSW*aR<NUMN>*TRES>*<TFKNUMN>»*2+TRES««2 >+HCOHV 


W  1 


Pfwfc  no  r 


do  TO  1873 


1871  IP  (NtMf-NBM)  1870,1874,1874 


1874  KI*I-NO 


IP  {NISO(KM-l)  1870,1870,1872 


*• ... - - — « _ L_ — . _ 

1872  f>UNCH  581,I,TH,SO(I) ,  (DEP(KI,J)',J-1,N) 

"  r  “ 


1878  CONTINUE 


188  NRITE<K8CT>  TH.TS.T0,0B».NIS0 
WC»MC«1 


ISO  CONTINUE 
N-HT 


IF  <NCON>  SOU. .30 18. 3008 


f* 


8003  1*1. NL 


IF  <HT>  SOOT. 800T. 3008/ 

;  /• 


3008  CflLL  LO0K<31*M,X<I>.TX<l.N>.Fl<l.H>.F2Cl.M7.0.0. Y1.D1.2) 
CN0< I>-Y1 « )«CN< I >»Y1 <8 )*CNC< I ) 


60  TO  3009 


i«< 


3007  CNO<I>-X<t>«CN<I>t<l.OHUI>>«CNC<I> 


3003  CONTI 


IFCKKU  3083.30X3.308t 


SOIt  DO  8030  1-l.HDBU 


V 

1 

* 

r  ■ 

* 

' 

• 

t 

. 

♦ 

*  , 

.  - 

\ 

i 

. 

4  « 

1 

’  ■' 

" i 

„  11 

* 

, 

lot* 


can 


pme  no  t« 


CIO  IF<01H-.000001>  ict.urcoo 


Itt  MUTE  (K0UT.502>  TH. OTH.OTHS.DTHB. THOS.  DTS,  OELCR.  DSOl 
’  TM-THFIN 


60  TO  3000 


INTERNAL  OECOWWITION  —  0EN6ITY  CALCULATION 


«00  N--JFM 

SPECIFY  SURFACE  CHANGES  DURING  THIS  TINE  INTERVAL 

DSDT«050TB 

JJS^OSOT*DTh 

OSI-1I.O«OS 

SA*SB*OSI 

RSU*ABS<RSV«SA) 

DTMB-OTHC 

OBL<NL>*OEL<NL>-OS 

FK»0 .0 

FJF-FJFN 

JI-JFHP 

DEN0LD=RH0<2> 

cold«cpe<d 

P0LD-CPE<2> 

CPE<  i>*RA<NL>««  .«OEL<N.) 

CPE<2>-CPE<1> 

IE*  I 

ROOZ-O.O 

ISV*WTL<NL*1> 

HATL<NL*I)»nATL<NL) 

TB<M.-»I>-TA<NBH> 

OEl<NL*l>*OEL<NBH> 


00  252  I-l.NL 


onoeci>*o.o 

RON<I>*«OOZ 

'■*  OSS*F‘JF/'OEL<I)*OSOT 


00  255  J-Jl.JF 


N*N+Jl-l 


\  . 


3-16 


1  ,,'  H“  ir  *•*».* 


icrNNn*  -ta>  •  ’mi-#**.  f*i  ■  **  ' 


*  '  v  -  /  ,  ’  •*  7  V •”tiv  '# 

/  *  .  .  .  .y  .  r ..  e  S',  -i  ,  *  . 

•  v  -  .  >,vr;  .  :'V  - 


can 


mw  no  is 


IF  <J-JFFP>  953.955  253 


|  <59  IF 


\ 


U- 


<NRTUIVHm<I*l>-0  M9.2I0.MI 


Ml  IF<flMRXt(Tfi<I>, Tft<J+l»-TRflCf1)?IUMl.MI 


Ml  IF  <N8TL«>*HRTl<I)*»«n.<I*t>-9>  2M.2C0.2S2 


200  J1«JFHF 

-  —  ■  —  "  - ; - i 

ROOZ-f  JFM*«XI*1> 

RONCD-fOK  D+fJFl+NKK  I) 

IWrWFH  s 

V. 

80  ID  M* 


U 


Kt  8002*0.0 

DTO*<Tft<I*l>-TO<I»/<FJF8^JF/OEL<I>*OOt<I^l>/WKI)«lwa*l» 
Jl«l  •  . 

TW8*Tft<I> 


MS  N*H*1 ' 


IF(flB8(D98>- .000001)  1011.1010.1011 


1010  DRonc-o. 

.  DROBOO. 
OROCq-O. 


80  TO  1021 


10  If  lF<I-NL)  i0  IS.IOIS.  10 15 


io is  rK-FK-i.o 


CM 


MO  II 


'  t 


11 


till  0NOBT-<H8>«(<NDmNH>- 

i«T7poh»/oth 


FB*€XP<  -€■/  TH8>*OlM>** 


SO  TO  <19* 


11 


tut  oM^MoSB^SiofHtSrra/iHPwT^Tonr 


l 


ri33  ROB<N>-«OB<N>+C 


:oKvr«aR08C>«oTH' 


IFCR08 ( H > -RHOR9?  2 133.223. 223 


I 


2 133  RO8(N>-MB<N>~<0*OT4aROBC)^OTO 
DROBT“(RMDRB-ROB(N> J/DTM-OROBC 
«08<:N)=4tt«R8 


11- 


2(3  IFCTflS  -TRHCC)  20S.MS.20l 


205  DHOCT-O.b 


11- 


eo  TO  2193, 


201  1FCRH0RC-*0CCN»219. 205.209 


213  RD°ROC(N)-«HOBC 
POM-l.-PSlC 


H- 


IF  <P0U)  2191.2192.2191 


2151  OROCT»(-SO*<<RD**FOM>- 
1<1./P0H»/0TH 


FC*UP«C/  TflS>3PTH>~ 


60  TO  2193 


* 


ft( 


3-20 


IF  <NDBU>  71fl.71tl.71ZZ 


710  1*1. 


7150  L*NBUFT(I) 

CALL  L00K(J1'*'L.X<J>.  TX<l,L),Fl<l,L).F2<l.L>,0,0,  Y1,D1,2> 
CN(d>=Yl<l>«CN<J>+YlC2>*CNC<J> 


mo  ssh2S=bsw 


f 


PflBE  NO  2S 


25  TERW-<-FFCT2*CPefiS-0S0T»R0l4CPl/0Q.<O>»0TH 
C<D=€<I>+TERK3 

D(I>«O(I>+TR<l+l>*TERH3+<FRCTf*H6RS*0SDT*ROl»Hl/DEl<D>*0TH 
TB=TB-TN*OSOT#RR< I> 


20  CONTINUE 


TT-TT+TWJTH/WEfK  1  >*ASU 

NON  THE  LAST  RBLHTING  NODE  REQUIRES  DIFFERENT  TREATMENT 
OROOTD=-OMD6<M_>/  <RR(»L>«OEL<  ft.)  > 

FHCTIOTH/  <DEUNL>«RR  <NL  » 

R<NL>=-FACT1*DVB 

DVB=1.0/<0.5*<RAT<N_>+RHT<NBH>)+RC(NL>/RR<NL» 

C<NL>=-FACT1«0VB 

TERM2 -RON  <  NL  >  *CP  <  NL  >  -  C  CPSHSMKOCmH  <  DSOT/ UEL  <  NL  »♦ 

<RO  (N_>«CPNL  -R01«CP1»«0TH 
B<HL)=TERH2-C<NL)~A<NL> 

D<NLJ=TR<NL)*TER7t2+OTH*CH6RS*ORODTD-mAR*<RON<lO- 

RD<NL>>/DTH*DSDT*CRO  <NL>*HttO-fi01«Hl>/OEL(M.>> 

K  =  NL 

6Sn=eSM-OHDBiNL> 


1FOC8U)  ri70.rm.TiTo 


TITO  FflCT2=eSM/<CCL<NO«RR<ND> 

TERttt=FRCT2*CPGfl6«DTH 
B<NL>=8(M_)  ♦FERT12 

0<NL)=0<K.)*rfKNL)*TEW1Z-FfCT2*HB«S*OTH 

CALL  LOOKC2. TR<NBH>, TTL, TH6. 0. O.O.HBAS.CPBAS. 1> 

H6RS-H6HS+0EI.H6 

C<NL>=C<«_>  FflCT2«CP6AS*DTH 

0  (  NL  ) * 0<M->  -  TH<NBM>«FACT2  *CPBHSN3TVHPRC  T2«H6flS«OTH 


TUI  CONTINUE 

: - NON  FOR  DECOMPOSING  BACK-UPS  IF  ANT 


IF  <NOBU>  7250. 7250. 7251 


7251  00  720  L-l.NOBU 


LL-NFHO-1 

LU=N-R<D 


00  725 1  t-LL.UJ 


3-28 


OROOrD»-{M)6<I  >/  <HR<  I  >  *OEL  <  I » 

FICTl^JTH/<UQ.a>*fiR<I» 

FflCT2«GSM/C0Q_Cn*«RCr>> 

AOC)=-fHCTl*OVB 

dvb=i  .0/  <0 .5»tiwrc  o+wrrc  i*  i>  j+rccd/rrc  i  >  > 

TERM*  I  >«CP<  I  )-0TH»<CP8flS»<DRGDTI>-fflCT*  )  ) 

TEnm*^ncTi«ovB 

BOO-TERM* -flOO-*TERMl 
COO-TEHMl 

DOC)=TFKI)»TERW*<HERS*DRDOTD-WHR*<RDNCI>HW<I?)/OTM 

-FflCT*«H0flS>»OTH 


ira-NUKMDeu>>  itt.iti.iti 


<I*l>*<i.0-X<I*l» 


WHR-P<L)«»»,CI*l>-l>PCL>/RH0VCL)*HCcr*O 
CflLL  LOOK  O.IWItl).  TTl.THB.O.O.O.HGflS.CPGRS.  1) 
HBRS=HGRS*0ELH6 


IFO-NFlO.)*!)  ««-T«.?*r 


COO=COO*TERM3 

D*IO-0OO*TH<Itl>«TERMJ*ERCr?«HGflS*0TM 


T^«  CflLL  SSUTCM<S.KKSU> 


GO  TO  (mO.rJll).KKSU 


<-{10  URtTEOCDUT.WQ)  I.K.L.OEL<I>.RR<I>.CP<I>.CPV<I>.Cf*C<n.W»<t>. 

lHC<I).CN<O.fWrci>.fl0a>.K)N<I>.X(O.0fl0Ga>.0fl0DTD.P<U.Pfa>. 

SHBRS.CPGflS.OVB.OVBS.MHt.HBnRS.eSn.FflCTt.FflCrt.  TEflrtl.  TERM*. TOGO. 
S8fl<L>.  OM0FKL>.  fl<K> .  BOO ,  COO .  DOO 
???0  FCRMRT<5X.3I5. 10E10.3/«0X.  10E10.D) 


i-29 


DO  90  l»*.K 


nno 


(HtT.  UlD.fTIS 


n 


tt«IL0<15> 


m  tfa>«mi>«vRN*m<i4C>-rfa)> 


mt  iF(EnMms)  wmnmiu 


CBM 

1 

PRGE  HO  40 

H 

60 

1  7 

OCOW-OCDNV  *VRP«  (  OO-OCONV  > 


1435  OCHEM“(OCHEM*OCONV>'*<» 
HH=OCONV 

QCONV-CM-CHE-CCONV> 

DS0TB=CHD/RH0<2> 


irtHBPf)  M5r,ior,*ST5 


4375  IFCIRB)  1437.1437.4376 


4376  CMFL-eXPCTJ  (3J)«CM^MH 


OO  55  1-8. H_ 


ROO  -  SOHO 

TEMP  -<0<I  )-fXI  >«TWI-1»/B<I  > 

OEDT-OEDT+RON<I>'*CP<I>*<TE>f*-TR<I>)*OEL<I>-*fiR<I>/DTH 


55  TRO-TE7P 


IF  <NUMN-WM>  57.S6.SS 


36  K=HL 

TRChBM-O  =  TR<HL> 


DO  31  I  -  HBM.NUMN 


i. 


ItPOUT 


PHBE  NO  1 


L 


SUBROUTINE 

INPOUT 

c 

— —  fiFF 

INPOUT  LISTING  FOR 

COMMON 

•MM 

c 

DIMENSION 

c 

»HM 

EQUIVALENCE  STATEMENTS 

Wf< 

c 

<»MB 

FORMRT 

■'ll.  IPi  "  ;l  ■ 

c 

MMB 

OHTR 

INPUT  =5 

KOUT-S 

INCH-5 

UT=0 .0 

_ _ _ 3 

\ _ 

_j  00  3500  1-1.5  | 

/ 

_ _ _ / _ 1 

f 

J  3500  NBUFT<I>=«  | 

_ 3 

l _ _ _ 

C  INPUT/OUTPUT 

C  / 

1  NP6=1 

WRITE  (K0UT.551NP8 

_ 3 

> 

REHD  < INPUT.  431.  END=333>  <RECQRD<I>.  1=1. 3S1 

■  '  _ 3 

. _ ....  -- 

BO  TO  330  |j 

_ _ _  •  '  J 

„ . 

_ 1 

333  CONTINUE  j 

_ _ _ L_3 

3 _ ! _ 

STOP  | 

»•  CONTINUE 

WRITE  <K0UT<S02XREC0RD<  I >.  1=1.30 
WRITE  <KQUT.503> 

WRITE  (KOUT.504) 

WRITE  <KQUT.505> 

WRITE  <KOUT.50«> 

REHD  (INPUT. 5S1> 

1B.PSIB.EB.TRRC8 


.RHORB.H 


R9 .  RHOOR.  RHORH.BR.  PSIR.  ER.  TRHCH.B9.RHOOB 

_  _ C9 .  RHOOC,  RHORC.BC.  P8IC.EC.TRHCC.NDBU 

WRITE  <KQUT .  507  >R3 .  RHOOR .  RHORR.  BR.  PSIR. EH. TRHCfl.  B9 .  RHOOB.  RHORB.  BB. 
IPSIB. EB.  TRHCB  . C3 .  RHOCC. RHCRC. BC.  PSIC.EC.TRRCC 
RERO  (INPUT. 5I3>  JF.NUHN.NN.NO.NI.THZRO.TfFIN.DTPRT.DTPRi.DTPRB. 
1DTHB.BRP.TPR2.TPR3.0EU1.DH1.DH2.DELHB.BRMR.TZ 


_ 3 

l _ 

IF(JF- l>  in.J70.lT2 

_ _ • _ I _ 3 

t _ 

170  JF=2 

3-46 


iw»out 


PflBE  HQ  S 


imvr 


PRBE  HO  • 


TOW  CONTINUE 
<09  OELOtMN>*OEL<MUHN>/19 . 

CMJL  ajKXMtt4.RR.HREn.EHR> 
mm  nl»i 


ifocbu)  root,  row.  root 


j  root  lg«H0JK«BU><T' 


h 


00  TO  TM1 


row 


root  CONTINUE 


-6- 


IFCRSV)  1091. lOSt. 1033 


«0X1  RSVMrHWV 

MUTCCKOUr .  991 MSVN.  RE 


00  TO  901 


1099  WITEIKOUT.ISI) 


ISO  IF<KT-t>  9SO.SM.riO 


■ENxueur.sri)  mc.tt5<it.i>.tc«u<it.o.tk8u<it.i> 


Ir<NC>  ri3.TW.ru 


rtf  IHKIWO-IT 

ILOUWO-l 

IR(I+ti>»l 

TEMT<1.T)=0. 

oo  ru  j^.it 


TU  1ENT<J.  O-TEWKJ-l.O+OCSUCJ.  D+TCOUCJ-t. 1»/t.«(TT5(J.  I>-TTS<J 
i.I» 


CALL  LOOK<£i»I.7REF(K).TTi(l.  n.TEWTCl.  D.0.0. O.HSH.OUM.  O 


US  TEHT<J.  0=TENT(J.  O-HSH 


if rir.no.rir 


TIC  RWWWOCOO 


BO  TO  Til 


nr  mommovoo 


n«  CHJ.  LC0LKT<0tIT.LCT.M>6.REC0R0O5» 

L-r*n 

HRITE(KOUT.52IO)  L.RRR.  <TT5<J.  O,  TC8U<  J.  O.  TKBUCJ.  I>.  TENUJ.  O. 

IJ-l.  IT) 
i*ru 
IT-0 


IF(NC)  r 19.353. 39) 


INPOUT 


PAS E  NO  1* 


PYROLYSIS  an  ENTHHLPY 


353  NT1-0 


311  IN-l*MTl 
NT144NT1 

READ  < INPUT. «r«  >NC. <TT1 ( I> . 1-IN. NT1> . <THBCI >. I-IN.  NT1 ) 


3M  NTl^lTl-1 


SOT  IUK33-1 
IR<»-1 
IHK»-NT1 

CALL  LGMJNT(9*«NTl4t)/5).LCT.NP8.REC0RD<35» 
3 IS  WRITE  <K0Ur.5tt> 

IFIM 


INB<NTl.IFNOT> 

TE  (KOUT.ftl>(TTt(I>.I-IN.IFN  > 
IE  <X0UT.S3S>ITHB<I>.1-1N.  IFN> 


INPOUT 


PMC  NO  IS 


3-60 


3-61 


INPOUT 


PR8E  NO  !• 


nmuT 


PH8E  NO  to 


INPOUT 


PRBE  NO  21 


<•13  CULL  L COUNT (11+2 ♦NR, LCT  .NP8.RECORD<35>> 

HRITE  (KOUT,53f! 
write  (Kour.srar)  cmh.hlo. vfz 


-FI.3.1XUH*  CHRR  THICKNESS! 


iiraur 


mo  n 


I 


Itr  IF<NST>  ttl«O.2ft«0.13*O 


I 


U 


Milt  CALL  LCOMT<LU.«.IXT.m.fCCORD<ff» 
WRITE<ICOUT.5fl2>PtV 


-C 


00  til*  LL-I.LLL 


I 


Mlf  MtriE<KOUr.SrM>  CTTSDKL.  »»>.  THSOKL. IPJ.L-U..MSDI-LU.) 


•o  id  mt 


mo  NSBt-0 
IX-f 


if <0*4-1  .>  tKi.tm.tm 


mi  i  x-* 


IF<MLO)  2t2«.2t2C.2t2« 


mi  MtITE  OCOUT.STJJ)  IX 


STOP 


»»  IF<IP-I>  IKMKl.mt 


■JJ 


tin  if<tts<j.mp»  wn  iHMiH 


2129  IFCPSY-TPR<IP»  tm.2t30.2tl2 


I 


3-66 


IMKIPW 


[tml  Hen-T9EN<K> 


INPOUT 


pnec  no  z« 


INPOUT 


pnee  no  25 


3-69 


tin  HBITE<K0UT.8fM>  BREX 


IMPOST 


PB0E  NO  >• 


|  EM) 


3-72 


6 

r — — — — — — 

MM  XFOC-MMO 

I _ : _ i _ ! _ 5 _ 

ten  oem-tyi<i.*-tk.  - 

■  .  «rtiw>r  .  ’■ 

atr<t.n>-*AlLT-*.*€B.lL>rM/OOt 

_ : _ ; _ i _ t _ : _ a _ _ _ ! _ _ _ .... 

V  r 

*• 

*  .  ' 

t.  j  \  .  ' 

[  (IS  1WNDTM,  _  ■  “ 

_ _ _ _ _J 

MHh.MiliifirVt')]-' 
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SECTION  4  . 


LISTINGS  OP  FORTRAN  IV  SOURCE  DECKS 


Listings  of  Fortran  IV  source  code  decks  are  presented  in  this  sec 
tion.  A  driver  routine  which  allows  all  computations  to  be  performed  by  sub¬ 
routines  is  presented  first.  Then  the  three  routines  unique  to  this  program 
CBM,  INP0UT,  and  THERMS  are  given.  Lastly,  utility  subroutines  used  by  the 
three  main  subroutines  are  listed  in  alphabetic  order. 


C  MAINLINE  DOM Mf  ROUTINE  WHICH  DRIVES  CBM  AND  INPOUT  SUBROUTINES  MAIN  1 

1  CALL  INPOUT  MAIN  2 

CALL  CBM  MAIN  3 

GO  TO  1  MAIN  4 

END  MAIN  S 


4-2 


1 

2 

3 


f  SUMMtxJl  lME  CbM  CSM 

C  Ouu-RING  <UIt«lM.  THERMAL  RESPWSE  AM)  ABLATION  PROGRAM  tUIMIMC.  COM 

C  F  O'  JP  TO  Five  JECONPOSIMG  BACK-UP  XAJLa  IJU_%  CSX 

C  aeRuTmEnm  CORPORATION  km  almiall  C  mote*  CSm 

common  aGuT. ICa.uEn.vR  CSM 

COMMON  IhI  OS  I  .iLO(3dl .|R(36>  .IT2C30.10I.ICPI 30. lOl.TKPiJ*. 10)  .TmZCBm 
lC30.18>.TEPC30.14I.TlMC7#I.ItCI30I.TORC30).TCN(3#>.Tll 130) .IMGC30ICSM 
2.DHI2I2)  .  RECORD  (jo  I  ,SO(?*l  CSM 

COmON  RiOllol  CSM 

.Common  haTl<  SOl.DELI  S|I|(M  SI) .hcSO) .RCI50) .RAC  SOI  *  CSM 

-  1APEA I  SOI -EMA(  SO) .RAVI  SOI  CSM 

common  R00(  SOOl.RDBI  So 01  .MOC(  S00.  CBM 

C-v*kXi  TPM  1 31  «)RM..  IDIOT  (21  .TNGC15.3I  .NL0  ( 1S.3)  «IOil  (  IS*  31  *  CBM 

1*hI  us, 31  .TTSCNCZS.3)  .ThS£N(25.3)  ;TCPS£N(2S.3I  .TLMCC2S.1S.3I  •'  CSM 
21  SEN  (3) .TPI (301 .1  IS :2S. 15.31 .TCHENC2S.1S.3) .VFZ.CMH  CSM 

COMMON  TSFFI2S.lS.3l  ..  CSM 

<C0*9tCN  •#>*>  CSM 

COMMON  LCT  .MPG.1I  .NBM.Nip4N.M-.0ELHG.DELH.Rri  .RMORA.RrtORB.RMORC.TRACSM 
lCA.TeAC8.TMACC.Mrt00A.RHo0e.HH00C.EA.EB.EC.B4.88.BC.PSlA.PSIB.PSIC.C8M 
2IRACM.PET  .PETE.MSV.ETA.DTPR3.DTPR2.0TPRT  •TWC3.TPTt2,TK2RO.TK'lM,sf  .CSM 
3TMMT .GAMA.OMG.NGi.F  JFH.F-^S.  JF.  JFHP.  JFH.  iiput .  DTHIN.BRP.HCONV.CBM 

♦EPS*.  IRES  CSM 

COMMON  INCH.OTHb  CSM 

COMMON  tM.Hl .mol  '  CSM' 

Common  ChcrI.PyLRI  Com 

COMMON  TBRPC30)  f.  CSM 

.COMMON  Ml  CSM 

3 COMMON  TX30.-6I  .FI  130.61  .F2C30.0)  CSM 

V Common  mcon  com 

(common  NBPF.HF1S  ‘  •  .  CSM 

common  BREA,  Smell  CSM 

common  BBS  IS. 31  «E£ (S.3)  ,FF  IS.3I.PSI  IS .3)  .RM00CS.3I  .RMOR  (5.3l  .  '  CSM 
1ROCOM(SO«31  .DmCcSI  .Om4is1.KHOCISI.RMOV (51  .pcsi.ppisi.trefisi.saisicsm 
2.0MOAC51  .MF1IS1  «NLA  (S7.TTS I 3S.101  .TEMT  |3S.  1SP.TKSOI3S.  101  .TC8UI30.C8M 
31*1  .X(S«I  .MCieo. MM2  CSM 

common  TRACIS. 31  CSM 

COMMON  NSuFT  (SI  CSM 

COMMOn/OATA/ASUM.  BLANK  CSM 

DIMENSION  uttOTlJl  CSM 

:  DIXNSION  CMOISvl  CSM 

i  DIMENSION  0MCPI3I  .CPEI3I  CSM 

OUCNSION  TO (201  CSM 

OltCNSlON  DEPI10.SI.  NISOIIOI  CSM 

jOItCNSION  VITERIS1I  .EITFUIS1I  ,  CBM 

OltCNSlON  T2I2NI  .02(241  f 73(81 .03 (SI  .T1I23.01I2I  CSM 

DIMENSION  CPC  I  SOI.CPVC  SSI.CPI  SOI  .HP  (  SOI.HCC  SOI  .CM  (  SOI.  CSM 

IRAT (Sll  .ROT (SOI  .ROISOI  .PON(SO)  .00)61501  .RRtSOI  •  'N  CSM 

2CNCI  SOI.* I  SOI. 01  S0I.C1  501.01  50) .CMO(  SO)  CBM 

E oul VALENCE  (Orti,oHl2(li )  .  (0M2.0M12I2I)  .(TS.TA)  CBM 

S2V  FOmMAI  I|)h  ITEKATIUN-STOR  1  '  CBM 

5A2  FORMA1  (/33A1AM - OUTPUT - 1  CSK 

SA3  FORMAT  I//BA2BM - - - F9.A.37M  SECONDS  -  -  -  CSM 

l-  -  -  - - -  -  -  -  I  CSM 

SNA  FORMAT  (6X.AMTIME. 2*. AMSURF.2A.AMPR0S.2X.TMSURF ACC. SX.brtH  BALL.  CSM 

lAX.krtH  E0CE.6A.  10.NRLAT  COEFF.kX.0MCM/CMO2BA.AMSTEP.2X.AMlTER.2X.AMCSM 
20PTN.2X.BMRAD  ( IN) .3X.BNIBTU/LB) .2X.SH (8TU/LB) . 3X. 1 AM ILB/SO  F T-SECCBM 

3)  I  CSM 

SaS  FORMA! (6*. IA.2Ik.Fl*. A.F11.2.F10. 2. F12.A.TX.F0.5/1M  I  CBM 

SAN  FORMAT  I33X.20M - ABLATION  RATES - 1  CBN 

SAT  F0RMA1  (6X.7MB  PNlME.3X.9tAV  PRIIC  G.3A.14HN  DOT  CMAR.3X.9tM  DOT  GASC8M 
,  16X.6HM  CMAH.TX.SHN  GAS/TAX .  1AMILS/S#  FT-SEC I  .  1 1 X .  ISM  (LB/OR  16  SO  FTCBM 
-  2)1  CSM 

SAS  FORMAT  ISA .F6.S.2X.FS.5.A I3X.F 10.6I/1H  I  CSM 

SA«*  FORMAT  I27X.32M - RECESSIONS/RECESSION  RATES: - /  CRM 

133X.19HIIN)  /  (TN/SeC)/  CSM 

f  16A.7HSOMFACE.lkX.6MCMAM  (.FA.2.1H) .11X.I IMPTROLTSIS  (.FA.2.1H) I CSM 

saoT  format i5x.3(ax.fio.t. in/. F9.T)/ih  1  csm 

Sa02  F0RMA1  (27I.3TH - SURFACE  ENERGY  FlUA  TERMS - /25A. J7MC0RRENT  RXTESCSM 

1  (BTO/SQ  FT  SURF  ACE-SEC)  /2AX . 3BMAM0  INTEGRATED  VALUES  (BTU/CRIG  SOCSM 

2  FT)/  CRM 

J3  13X.10MC0NVLCTED  .Ax.  ION  RADIATED  .AX. TOM  RADIATED-  .AX. TOM  CHE  CSM 
.MICXL  «AX  .  10MC0NDUCT  ION/1 7X.2HIM.  12X  *2MI  V.  1 1  A.  3M0UT  .  BX  *  1 OMGENENAT I  CSM 
50N.7X.AMXMAT)  CBM 

3  FORMAT  l6X.AMftA(£.3A.S(El«.3.AX>/64.SrtTOTAL.2A.SIE10.3.AX>/lM  )  CSM 
SABA  FORMAT (30X.Z7H - INTERIOR  ENERGY  TERNS - /2SX . 37HCURRENT  RATES  ISTCBM 

(1U/SD  FT  SURFACE-SECI/2Ax.3Bt(Af<)  INTEGRATED  VALUES  (BTU/0R1G  SO  FTICSM 
2/  CBM 

3  13X.9t*»YK0L  GAS.7X.6*C£COHP,6X.14HCONVECTION.6X.7><STORAGE.  CSM 

A7X.THL0SS  AT/1AA.7NP1CK  UR.kX.  1 OMABSORPT  I0N.3X.J INHITN  SOLIOS.SX.  CSM 
SOMlN  SOLID. 0X.9MREAR  FACE)  CSM 

SABS  FoRMAI  (6X .AMRATE.3A.S (E I0.3.AX) /kA.5MT07AL.2X*S  (E  10.3. AA) /1M  I  CSM 
5A9  FORMA)  CkXBHNOOC  MAT3XAMTEMP3XTteCNSITV3XSHENTMALPY2iatWOOE  MAT 3AACBM 
1MTEMPjA7mOENS1TT3ASHENTkALPT/1SATMIDES  R)IIM  CLB/CU  FTI9M  (BTU/LBICSM 
21 lATHlOEG  RI11H  (LB/CU  FTI9M  IBTU/LBI)  CRM 

SA90  FORMAT  IkX.CRMOOC  MAT  .  3X.AMTttO>.  3X.TH0CNSrT  Y.2X.9HCONO  (OTU/.  21  .OMNOCBM 
IDE  MAI.3*.AMTEMP.3A.7OCNSnT.2A.9MC0ND(8TU//lSA,7M(D€&  RI.I1M  (LBCBM 
2/CU  FII.9H  FT  SC  F ) .11A.7HC0EG  RI.11H  (LB/CU  FTI.9H  FT  SC  F)l  CBM 

550  FORMAT  I 5X2 l A «F V « 2. F 10. 3. 2XF0 .2. IA2IA.F9.2.F10.3.FI0.2)  CBM 

55*0  FORMAT (5A.2IA.F9.2.F10.3.2X,FB.6.1X.21A.F9.2aFI0.3.F10.6)  CBM 

551  FORMA)  (lHllOXGSMALROTtCRM  ChARRING  MATERIAL  TtCRMAL  RESPONSE  AM)  ACBM 
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32?*  MITC  IXOul . V>*il  OCOMV. •NP.tAO.tOCP. DCOM). DCONVT.MPT. MOT .OCMCMCaM 

29* 

it.oconot 

CM 

295 

MITC  IK0UI.S*t*l 

CM 

296 

MITC  IKOuT.SMsIPtPU.DrCOM.  IS.OEDT.tLOSS.PSPUT.OCCMt. 

TT'OEOTToCMM 

287 

i gloss i 

CM 

296 

ir  IMCONI  Jtn.36M.362l 

CM 

299 

3*71  Mitt  ((OUT. 5*3*1 

CM 

366 

CO  TO  3*22 

CM 

361 

3626  •altclK0UT.S*9l 

CM 

382 

3622  CONTIMUC 

CM 

363 

IFIM01I  1H.IH.III 

CM 

36* 

163  CM3.  SLOPO  (NL.HAIll  .1*1  II  .CM0I11 1 

CM 

36S 

iriMOl  162. 162.166 

cm 

384 

166  C*Lt  06tC  INO.SU. T0UI.M..UI1I.T6III.CMI1II 

CM 

361 

162  IF  INI1  166.166.ltS  v 

cm 

366 

16S  00  166  1*1. ML 

CM 

369 

166  CH0lll*l./CN0lll 

CM 

316 

CALL  06t£  INl.SOIMO.il. TOlMO.ll.NL.tA.aA.CNOI  ' 

CM 

311 

166  lFINUMMHHMIttt.6tl.ttl 

CM 

312 

661  CALC  T MCIIMS INISO.T0.0CP.NM.CN.t6fl 

CM 

313 

■66  IFINNI  16t.ltt.ltf 

CM 

31* 

161  00  1613  1.I.N01 

CM 

315 

IFIl-MOl  16l6.lMft.10fl 

CM 

316 

1616  PUNCH  SMI  .1 .  IN,  SOtll.  TOIII 

can 

311 

00  10  1873 

CM 

316 

1611  IF  INUM-NtM)  Iel6.ltf6.1tf6' 

CM 

319 

1874  KI*I-ftO 

CM 

326 

IFINlSOIKII-ll  16ft.ltft.ltf2  \ 

CM 

321 

1812  PUNCH  Sai.l.Tn.SOtli  •  lOFPIKI.JI  .Jm1.ni 

CM 

322 

ItlllCONIlHUC 

CM 

323 

188  TH.fS.IOtOtP.NlSO 

CM 

324 

**■**♦! 

CM 

325 

188  CONTINUE 

CM 

324 

NsMT 

CM 

321 

IF  INCOM)  3612.3612.3666 

CM 

321 

3884  00  3808  UloNL 

ca* 

329 

IF  («U  3807o30«7.J088 

CM 

336 

3666  CALL  LOOK  I  J1HI.A41I *TX 1 1 .Nl .FI  1 1.NI .F2II.MI.6.6.T1.01.2I 

CM 

331 

.  CN0<l)*Tl(ll*£Mtl)»Tll2l«CNC(ll 

cam 

332 

00  10  3889 

CM 

333 

3661  CNOI ll*Alll *CN(II • 1 1 .*~R  Cl  1 1 *CNCI 1 1 

CM 

334 

3889  CONTINUE 

can 

335 

iriNDoui  3623.3623.3626 

CM 

336 

3824  00  3838  UloNDBu 

CM 

331 

LL-MFItll 

can 

338 

LUMNLAlII 

CM 

339 

LMOUFTIII 

can 

348 

00  382S  uNUtLU 

CM 

341 

1FIL1  3626.3626.3621 

CM 

3*2 

3621  CALL  L00K131.L.AI Jl « Till .LI tf l I 1 »LI *F2( l.LI .6.6.T1.D1.2I 

CM 

343 

CNOt  JI*T1 1  |)*€N(J)  *Tl  (2|*CNCC J)  ’jf 

00  TO  3825  * 

CM 

344 

CM 

3*5 

36M  CNOUl.llflMCNUl.ll.t-lUirKNCIJI 

can 

344 

362S  CONTInUC 

CM 

347 

3636  COM3 1NUC 

CM 

348 

3623  IFINUMNHWM2)  3612.3626.3626 

CM 

349 

36M  DO  3629  l*MDN2.MUMN 

CM 

3S6 

3829  CNO(ll«CNlKl 

CM 

351 

3812  CONTINUE 

CM 

352 

DO  3611  3*1. NLI 

CM 

353 

L»J 

CM 

35* 

IF  ILH6.I  3662.3662.3661 

CM 

3SS 

3661  L.L.NON 

CM 

356 

K»NLt 

CM 

351 

GO  TO  3883 

CM 

356 

3662  IF  1L.NLI-ML)  3663.3663. 3665 

cbm 

359 

300S  X.NLl.NOa 

CM 

348 
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M3  Ks«IM«<fftJNN»ft»L> 

CON 

361 

IF  INCONI  306*  o  3004.  JO  1 6 

CON 

302 

Ml*  M(K  (AOUT.SiOl  ll.HATllll.TAII>. Mill. Mill,  1-L.N.XI 

CON 

303 

60  10  MU 

CON 

304 

Ml*  ■PIICikOUT.SSOOMI.MTLiII.TOIII.ROIII.CMIII.IfL.H.XI 

CON 

305 

3011  COMTlMUE 

CON 

300 

IF  I St>t  1X1  JIM.jlll.lIM 

CON 

367 

310  *  SONET- (l..SOtlXI»OIOI-SotU.»OCOI 

CON 

360 

■RITEIXOUT  .59*1  KCI.SMI 

CON 

369 

wi  roMui(/i«i,4iMSdNfAa;  accession  after  shell,  umcksi 

•  roc on 

370 

I.4/!U«4«HSMrAa  ACCESSION  DATE  ilTH  StfCLL  (INCmCS/SCO 

■ro.41  con 

371 

31*1  CONTI Hu( 

CON 

372 

Call  SSxTCM( J.aasa) 

CON 

373 

SO  TO  1745.7461  .|CK  SO 

CON 

374 

7*5  IFINOOUI  M.IWdlT 

CON 

37S 

7*7  LU-M6H2-1 

CON 

376 

LfLU-nON-1 

CBN 

377 

U.H/t*2 

CON 

370 

COLL  L COUNT (L.LCT .MRO. RECORD OSI 1 

CO  i 

379 

DO  7*4  1-HUH.LU 

COI 

300 

740  ■RITE|«0UT.7*4I  (HOCONII.JI.X-1,3) 

CmN 

301 

MUTE  IX0UT.7491  (0ND6IJI  .J  MH.LUI 

CON 

302 

■RITEIXOUT. 7**1  (.SMS.GSM2S.GS.6SM2T.6SM 

CON 

303 

LU-HUHM  NOB 

CON 

304 

MUTE (0007.74*07  (OHO  .01X1  .CIKI  .010)  .A-1.LUI 

CON 

305 

749  FORM!  ll0A.6E10.il 

CON 

306 

7490  FORMTUOX-4E10.il 

CON 

307 

740  CONTINUE 

CON 

300 

mts  ir(TH-TMri»-o.o*«*miSi.2  .2 

CON 

309 

2  IF  INUII  1-1-3 

CON 

391 

3  RENIMI  XSCT 

CON 

391 

NM 

CON 

392 

CALL  LCOUNT  (-10.lCT.HPG. RE  CORO.  751) 

CON 

393 

■RITEIXOUT. 5521  N0-NI 

CON 

394 

■RITEIXOUT. SS21I 

CON 

395 

■RITEIXOUT. SS*>  TH.TS.  ISOlII  .Ixl.MI  1 

CON 

396 

■RITE  IKOUT  -SS3I 

CON 

397 

00  *  A-l.AA 

CON 

390 

Be AO (ASCII  TH-TS.IO-OEF.NISO 

CON 

399 

00  *05  1-1.N1 

CON 

400 

IF  (M-MISOtll  1  004.00S.04S 

CON 

401 

006  MfNISOIII 

CON 

*02 

00S  CONTINUE 

CON 

403 

COU.  LCOUNT  (  1  -LCT-NR6. RECORD  I3SI 1 

CON 

404 

4  (MITE  (K0UT.SS4I  TH.TS.  UOUl.  If1.MII 

CON 

405 

lFIMUHN  NONI  1 -*02-002 

CON 

406 

002  IFINKNI  1.1.003 

CON 

407 

003  REBIND  ASCI 

CON 

400 

COLL  LCOUNT  (-7.LCT.MS-0EC0R0 13511 

CON 

409 

■RITEIXOUT .5601 

CON 

410 

JfM-1 

CON 

411 

■RITEIXOUT.  70411  ISOIII.IfJ.MII 

CON 

412 

■RITE (AOUT. 7042* 

CON 

413 

S60  FORMI  (/0X06MTHL  FOLLOMIM  BLOCK  BITES  THE  OUTPUT 

TINE  AMO  1#  TO  fCOM 

414 

1IVE  LOCATIONS  /  OX-SSMOF  THE  IMICOTEO  ISOTHERMS  MITNIN 

TH E  OACKUOCON 

415 

2  MTERIOLS.I 

CON 

416 

7041  FORMT  I//I2X-9I20.F  10.4)  1 

CON 

417 

7042  FORMT  I/I 

CON 

410 

00  004  XfI.XK 

CON 

419 

READ  (ASCII  TH-TS.TO-OEP-Nl  SO 

CON 

420 

00  004  JFI.N 

CON 

421 

COLL  LCOUNT  1  1-LCI-NPO-REC0MI3SII 

CON 

422 

004  URITCIX0UT.S54I  TH- IOEFI I • Jl -IfI.NII 

CON 

423 

GO  TO  1 

CON 

424 

1151  1FITH-TPR2.0. 000011  154.150. ISO 

CON 

425 

ISO  DTPRT-0TPR2 

CON 

426 

OTPR2-OTFRJ 

CON 

427 

TPR2NTPR3 

CON 

420 

TPR3FTHFIH 

CON 

429 

154  CONTINUE 

CON 

430 

IFITH-THPRT*. 000011  743.7*2.742 

CON 

431 

742  THPRTaOMINKTHPNTHITPRT.TPREI 

CON 

432 

743  CONTINUE 

CON 

433 

4410  OTHfANINKOTHR.UELCR/IOWTO-. 0000001 1.TH-THOS.  50. 

||/(AIS( TSA Vt-T A  ICON 

434 

llll-.ll*OTHI 

CBM 

435 

TSOJFE-TOllI 

OTM 1 TMRT-THI  /  I01NT 1 ITMRT-TMI  /DIN*  1 .01  * 

COM 

436 

COM 

437 

144  THFTH.0TH 

COM 

430 

C 

COM 

439 

c  functions  of  tire 

COM 

440 

604  IfIRIII 

COM 

441 

TF»TF2 

COM 

442 

601  lFtTTHII.ll-TH.0.000011  614.604.6*4 

COM 

443 

614  IFILI-IHIUII  602.604.404 

COM 

444 

602  I«l-I 

COM 

445 

IFITTHI 1.1 l-TTHIl 1 1  6*1.603.601 

COM 

446 

603  ThfTH-OTH 

COM 

447 

ThOSfTTHIII-OTHIN 

COM 

440 

OTHbONOXKOTHIN.TTHIII-THI 

COM 

449 

ThfTH-OTH 

COM 

450 

\ 

\ 

\ 
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*1 

. — — — 

GO  TO  Ml 

CRN 

651 

M*  (fcM<l IM-I7HI i> 1/ IIIMII«| 1-llMlll 1 

CRN 

652 

■ 

IF  If \X 

CRN 

653 

ocn*o. 

CRN 

656 

MS  Crt*TCM||>*OCN«<ltN<I*|>.rCM<lM 

can 

455 

O0**7Jtf|ll»OCM«lloatl.li-I*n(lll 

CRN 

656 

P*CS«IP1UI.DEMMIPI«I.|I-TPHII» 

can 

657 

w(»IHt  (||  »OCM#  (TMC 1 I* 1 )*VMC (II  1 

can 

659 

SOPMlMP  ( 1 1  •OUf  1 T*«P  1  l.ll -TBHP 1 1 1 1 

CRN 

659 

11*1 

CRN 

6M 

IF  (CHI  MOl  aM91*M9 

can 

661 

6001  11*/ 

CRN 

662 

Crt**.0 

can 

663 

CRN 

666 

6092  11*3 

CRN 

665 

! 

*F-*C 

CRN 

666  ' 

j 

MF-f. 

CRN 

667 

699 

CRN 

66R 

IFdUfcl  619*  1 13*619 

CRN 

669- 

j 

||3  OtOTM* 

CRN 

679 

IT£«*1 

can 

611 

■ 

GO  1C  3009 

CRN 

672 

610  If  I0!«-.0#0»01I  162.162.694 

CRN 

673 

162  PSltl  IKOUf.Sa/l  fH.OTM.D7MS.DlNB.Tl«S.07S.OELC»,DSOI* 

CRN 

676 

. 

Trt*Tr*  IN 

CRN 

615 

GO  TO  3999 

CRN 

676 

C 

CRN 

677 

c  imtemmal  oecoMPubiuoN  —  ocmsity  cm.cu.ai  ion 

CRN 

679 

699  N*-JFrt 

can 

679 

C 

CRN 

6*9 

• 

C  SPtCIf  1  SUPfACC  CHAMCCS  WMIN6  INIS  TIME  INTERVAL 

CRN 

6R1 

OSOI-OSOIB 

CRN 

692 

OS*OSoT*OTH 

can 

6*3 

OSl-12.9«OS 

CRN 

6R6 

? 

SA«SA»0SI 

CRN 

6*5 

HSO*MS|RSv*SA) 

CRN 

696 

i 

OTH8-OTHC 

CRN 

697 

OELINU-OCtlNLI-OS 

CRN 

699 

: 

Fr*0.o 

CRN 

699 

F  JF*F  JFH 

CRN 

699 

ji-jtmp 

CRN 

691 

0CH0U)*RM0I2I 

CRN 

692 

; 

COLO-CPE ( 1 I 

can 

691 

POLO*CP€ (/) 

can 

696 

| 

CP£II>«IA<NL>.6.*0EL<M.) 

can 

695 

•  i 

CKI2IKKI1I 

can 

696 

1 

lt»l 

can 

697 

j  • 

9002*0*9 

CRN 

699 

( 

ISV«*ATLINL«1> 

can 

699 

6ATL<NL»1I«MAIL<I*.> 

CRN 

599 

i 

T  A INL ♦ 11 «T4 (NbMI 

CRN 

5*1 

[ 

0CLINL<I1<0CL<NM<I 

CRN 

592 

t  i 

fiRUC'll-MHNRNI 

CRN 

593 

00  Z5C  1*1 *NL 

CRN 

596 

OMOGI1»*0.9 

can 

5*5 

RONI  11*900/ 

can 

596 

*  OSS*T JF/DCLd>*U*OT 

CRN 

5*7 

*61  Oo  »i>  KJl.JF 

can 

599 

N*N*31-| 

CRN 

599 

IF  IJ-JfM>l  fS]i2Mi»3 

CRN 

51* 

f59  If  <NA1L<I1.NA?LI1<11-6|  Mi<2M(lU 

CRN 

511 

261  IF  I AMAX1  <TA<I).fA<t<l)  |  *I*ACM)  261 .261  <2*2 

CRN 

512 

261  IF  (NATL  d)*HATl<l)«lUTL <1*1 1-3)  2M*2M./62 

CRN 

513 

269  Jl*JF|iP 

CRN 

516 

P002*f  JFH*RO(I*l) 

can 

515 

RON  (If  *R0N  1 1 1  *F  jFH*RO  ( 1 1 

CRN 

516 

N*N*  JTH 

can 

517 

GO  TO  264 

can 

sia 

262  R 002*9*9 

CRN 

519 

OTAMiAii.ii>iA<in/<rjrs>rjr/OELiii*OELU*ii/miii*M*ii*iii 

can 

529 

JI*1 

can 

521 

TaS*TA III 

can 

522 

253  N*N* 1 

CRN 

521 

IF  i*asiossi-.#»®»»n  1*16. 1116. 1917 

CRN 

526 

1016  O0OAC*O. 

can 

525 

OROOCaO* 

CRN 

526 

OROCC-9. 

can 

527 

00  TO  1921 

can 

521 

1917  IFIl-MLI  1915.1919.1915 

CRN 

529 

191*  FieFK-1.9 

CRN 

5M 

1F1FK1  1919.1916.1929 

CRN 

531 

1919  FF-F  jf -1.9 

can 

532 

1929  tSS*OSOT/OEL(II»F« 

can 

511 

1915  IF  IM-ll  11I5.U1S.I11S 

CRN 

536 

■ 

1115  OPOAC<<MAIN<I  I -ftOAINH  >055 

CRN 

535 

M 

OPO*C<  1 00*  1 M*  1 1  MO*  1 N 1 1  *055 

CRN 

516 

B 

OPOCC.IPOCIN.II-MOCINMaOSS 

can 

517 

M  ' 

00  10  1921 

can 

S3* 

m  ■ 

1215  IF  1 1 1-1 1  1115.1115.1115  '  "s 

can 

539 

m 

Ills  DP0AC<IR0AIN.ll-NN0*AI<1)$S 

CRN 

5M 

1 
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1 

«r  ■<«.  . 

i 

• 

I 

cm  3*1 

CM  3*2 
can  3*3 


UMXM.  *  ( *08  <  N*  1  > -mows  >  *r>ss 
MOCC*  (HOC  (*•  1  >  -tmMC  I *033 

aoiili*Hnoi2i 

1021  TASMAS.QTA 

1FIT*$-TH*CM|  221.227.229 
22V  IF  11*3  -IM*C*I  201.201.202 
2*1  000*1*0.0 
60  TO  2113 

2*2  IFlMdOA-eoAInll  211.2*1.2*1 
211  00*00*  INI  «O*10OA 
000*1.-031* 

If  1000)2111.2112.2111 

2111  ueO*T*(-«O.|<O0«*OOO>-  FA-CAPI-CA/TASl-OTHI*- 

1 11.200*1  l/OTH 

60  10  2113 

2112  0*OAT<ttl*(tAP|-aA-OTH*tlPI-eA/TASD-l.l/OTH 

2113  00*1*1*00* (HI .  IiMOAT*OOOAC)  *0TM 
iFIMOAtNI-HMmAI  211**221.221 

2(1*  OOA<N>«OOA<HI-l0O0Af»0OO*CI*OIH 
000*1*  I0H00A-00*  (HI  I  /OTrt-OOOAC 
oo*ihi«ohoo* 

221  1FITA3  -10AC6) 2*3. 2*3.2** 

203  O0OHT-*.* 

60  10  2133 

20*  IF  (MMUOO-HO0  (HI  1213.2*3.2*3 
213  8D*NO* IHI-OHOnH 
00**1.-051* 

IF (00*12131 .2132.2131 

2131  00081.<-a0.(  (HO**tM»l-  Fa*£JU>(-£a/TASl-OTMI-« 

1  ( 1./00HI  l/OTM 

60  TO  2133 

2132  URO*1*aO*lClPl-aa*OTH*CxPI-£a/TASII-l.l/OTM 

2133  OOB<HI*«OBINI.|DMaT*MoeCI*OTN 
IFONtolHI-OHOlMI  213*. 223*223 

213*  KOOlN)*HW(N)-(DOO*TKttO*CI*t>TH 

ooo*t-  iomommmm  (HI  I  /DfM-OMOC 
(MMIHIaONOM 

221  IF ITAS  -TO ACC I  205.245.20* 

2*5  OflOCT**.* 

60  10  2153 

2*6  IF  langHC-MClMI  1215.2*5,2*5 
215  00*HOC  (HI -UMOnC 
P0--I.-PS1C 

IF  100*1  2151.2152.2151 

2151  OOOCT-I-*0.(  (00*— 00*1-  FC*£*P(-CC/T*5I*0TH>  — 

in. /com)  i/stm 
CO  TO  2153 

2152  OOOCTaOO-ICXP  1-*;*0TM*E|P l-EC/T ASI  I  -1 . 1 /DTH 

2153  «OCIHI>«OC(HI.UMOCT<MOCCI*OTM 
IFIOOCIHI-HHMCI  2154.225.225 

215*  M0C(N)*O0CIII>-lllHOCT<oaoCCI«OTH 
OOOCT*(OnOOC-HOC(H)  l/DTw-OHOCC 
OOC(NI*OHOMC 

225  01*161 11  ■OW6UI-mLUl*HOHO*THMO>Tl*)MH*«OH**OnOCTl 
2?r  OhS*<HO*|HI*RO*IH1I*MH*HM6«<IOCIH> 

MOM  ( 1 1  **0H  1 1 1  .«a 
1FIH-1I  225I.225J.225I 

2251  IFIOHV-OHCOIICII  2253..'  -52,2252 

2252  COCllCI*OU.IIl*liri.O*Tt  1-*.5l/f’JF.|.-FlMTlIll*l2..ll*(ll 
lF(**5(0«$-OCMacill-1.0*-2*l  225*. 2254.225* 

2254  CPCIItl*  COCIIE)-OC1.I1I/IIMS-OCHOUII*NM$-OHCOIICD/FJF«12.* 
225*  CPCIU.fMAXltCFtllEi.Sal 
IE*IC-I 

2253  OTMUMW 
255  IA5*T*S-0T« 

26*  0*06111*01*16111/1  JO  *0*1 1 1 
«OH(l>*«ONII>/Tjr 
IFII-H  25T.25T.J5* 

25?  FJT-FJFS 

0?A*01A/DC(.  11 1  »ua.  121/2. -WII2I 
60  TO  252 

25*  OTA*OfA/OCL(l)*UtLf  I-1)/H*III*MH1-|) 

252  CONTINUE 

IF  IKUOU)  2522.2522,2523 
2523  00  ?•(  l.l.WOu 
6*w  6*H) 

0H6AM*|.  -**N 

U.MVIUI 

LIMM.AIII 

i**CHI.*N1n1  ITOACII.il,  T**C  11.21.  TUCH.3II 
00  T*2  JHL.LJ 
OHOOIJI**. 

IFI1A1JI-I0ACH1I  T*2.T*4,T*3 
T«3  00  1*9  K*| ,3 

IF (TAI JI-TBACl 1.61 1  F *4.7*5.?** 

T*5  0*0? (*!■*• 

GO  TO  T*4 

TM  HO-HOCONI  J.AI-dHOOll.K) 

IF  IH0-*.*1I  T *5.7*5. T*4* 

1*6*  F-FF  ( I  ,KI 


can  5** 
Can  5*5 
can  5*6 
CM  5*7 
CM  5*« 


CBN 

549 

cm 

$5# 

cm 

Ml 

cm 

M2 

cm 

553 

cm 

M4 

cm 

55b 

cm 

556 

cm 

55 7 

cm 

556 

cm 

559 

cm 

566 

cm 

561 

cm 

562 

cm 

563 

cm 

564 

cm 

565 

cm 

566 

cm 

567 

cm 

566 

cm 

569 

cm 

576 

cm 

571 

cm 

572 

cm 

S73 

cm 

574 

cm 

575 

cm 

S76 

cm 

577 

cm 

576 

cm 

579 

cm 

569 

cm 

561 

cm 

962 

cm 

563 

cm 

564 

cm 

565 

cm 

566 

can 

567 

cm 

566 

cm 

569 

cm 

596 

cm 

591 

cm 

592 

cm 

593 

cm 

594 

cm 

595 

cm 

596 

cm 

597 

cm 

596 

CBM 

599 

CBN 

666 

CBN 

661 

CBN 

662 

CBN 

663 

CBN 

664 

CBN 

665 

CBN 

666 

CBN 

667 

CBN 

666 

CBN 

669 

CBN 

616 

CBN 

611 

CBN 

612 

CBM 

613 

CBN 

614 

CBN 

615 

CBM 

616 

CBN 

617 

CBN 

616 

CBM 

619 

CBN 

626 

CBN 

621 

CBN 

622 

CBN 

623 

CBN 

624 

CBN 

625 

CBN 

626 

CBN 

627 

CBN 

626 

CBN 

629 

CBM 

636 

£•€£■!.«> 

Wi>i.>niiui 
iriWal  707.709.707 

7*7  Mgii»<i-«*ii«MnMHauri-{/itiJii>iiiMMii./nini/tni 
60  70  7*9 
7««  ofmmhi.ii 

MOI  (1  1  MO*  (UP  i -WQtmtJP  i  -c/t*  u»  » -i .  i  /dth 

709  MCbM<J.«laMC4N(J.»l.M0TISI*0TM 

791  OOM(Jl«U»OCOM(J.11.0»lCOM(J.2H.*lM.OMg*M*BOCOWIJ.31 
702  CONTINUE 
700  CONTINUE 

c  aoo  wturr  necessmv  nth  post -decomposition  pmopebties 

2S22  J«-PHP 

N1TL(M.<1)«1S* 

oo  os  I'l.w. 

IF  Hlte(ltONII>-«M><lll-.(l)  01.01.02 
01  MATL(I**1 
*111.1.0 
60  10  OS 

02  IF  I10S(NONI1|-«MO(2II'.01I  03.03.0* 

03  N1TLI1W 

1(11*0. 

CN(ll*CNC(lt 
60  70  OS 
0*  MlTLIIIM 

■  IIIMfTl'Oll/MlIIl 
IP  (M7I  OSOO.OSOO.OM1 
OSOI  L-#7 

ClIX  L00MJ1»L.1II>. 7*11 0.1  .FK1H.) .roil. Ll.0.0.71 .01.2* 
CMIU-Tl  (|)*CN(il«  71 (2I*CNC  til 
60  70  OS 

osoo  CMiii«xiii«anii.ii.o>iimocHeiii 

OS  NlT(l|.OCL(ll/ICM(II*M(lll 
IF  (MUl  7121.7121.n22 
7122  00  710  U1.MMU 
U-MFtlll 
LIMUdl 
00  7||  J4X.LV 

1f(*os(*on(J)  MMH)(i;i-.01l  T12.7I2.713 

712  HITLUI*  20.2*1 

101*1.0 

60  70  716 

713  IFIKtStOOMOl  *>aC(lll-.01l  716.716.719 

71*  Ml  IL  01621.2*1 

101*0. 

BOMlJlXMOCIll 
CMIJI*CMCOI 
60  70  716 

715  <M7L(JI*0 
1IJI*P(II*PPI1I/OOMOI 

IFIIOUF7I1II  7ISI.71SI.71M 

7150  OOWIIII 

CALL  iM(3K.lUI.7|tlilOI)U)rftlHt.O.I.714ld) 
CNOI6T1 1 11*0(01  .Tl(|l  *CK  IJI 
60  70  716 

7151  C0M7IN UE 

CM(  Jl*l  (  JI*CM(  J1  *  ( 1 .0*1 1  Jl  1  *OK  (  J1 

716  NlT(JI.OCL(Jl/(CMIJ|lOMIJI> 
in  Coni  imuE 

710  CONTINUE 
7121  OSH*  0.0 

00  122  1*1 .ML 
122  *5N*0»M'Cn00U1 
WW4W 

69N7*w9N7  «*SNS»1M«/1MC6  ( 1 1  *OTM 
6SM2*0. 

IF  (MUOUI  7IM.7100.7101 
7101  00  717  UI.MOOU 
UMOlltl 
UMUIII 
00  710  J4X.LM 
710  •SM2^SM2«0M00tJI 

717  COMTtNUC 
7106  UlOUiMt 

IWT>OOWT4Ml.tl(l/MUI|l*OM 

6MWMM1 

HoOl 

UC07«  (CPC  1 1 1  -COL0I  /OTM 
OPOTKCPCIll-POLOI/OTN 

C 

C  CHOU7IM  or  HTL1CIT  TCMPEBITVBE  COEFFICIENTS 

C  MO  1MTEOM1L  110107  un  TOM 

C  Ml  IN  0L0C1 

070*0. 

sac6o*o. 

6SEM*0. 

70*0. 

■17(11.2. *097(11 

CPM.*CP(NLI 


COM 

Ml 

CON 

632 

cm 

633 

cm 

63* 

cm 

635 

cm 

636 

cm 

637 

cm 

630 

COM 

639 

COM 

6*0 

COM 

6*1 

COM 

6*2 

CON 

6*3 

COM 

*** 

COM 

6*5 

COM 

6*6 

cm 

6*7 

cm 

6*0 

cm 

6*9 

cm 

650 

cm 

651 

cm 

652 

cm 

653 

cm 

65* 

cm 

655 

cm 

65* 

cm 

657 

cm 

650 

cm 

659 

COM 

6*0 

CON 

6*1 

COM 

6*2 

CON 

6*3 

COM 

6t* 

COM 

6*5 

CON 

6*6 

COM 

6*7 

CON 

6*0 

CON 

6*9 

CON 

*70 

CON 

*71 

CON 

*72 

cm 

*73 

COM 

•7* 

COM 

*75 

CON 

*7* 

CON 

•77 

CON 

•70 

CON 

•79 

COM 

606 

COM 

601 

COM 

M2 

CON 

603 

COM 

6** 

COM 

60S 

CON 

606 

CON 

MT 

CON 

600 

CBM 

609 

CON 

690 

CON 

•91 

CON 

692 

CON 

*93 

COM 

•96 

CON 

695 

CON 

696 

con 

697 

CON- 

690 

CON 

6*9 

CON 

7M 

CON 

701 

CON 

702 

CON 

703 

CON 

70* 

COM 

70S 

CON 

70* 

CON 

707 

CON 

7M 

CON 

709 

CM 

710 

CON 

711 

CON 

712 

CON 

713 

CON 

71* 

CBN 

719 

CON 

71* 

CON 

717 

CON 

710 

CON 

719 

CON 

720 

.1 

T 
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rl 

ININ** 

00  3*  1«UUN.NLM 

i r  <11  is, is, i» 

14  GSMafeSH-OtOOdl 

COOC  1 0>-<ll(IC  1 1 1  / 1  AD  1 1 1  »OCt  1 1 1 1 
FAC11aOTh/tOCLlll*NN<ll) 

FACT2^SN/|0£Llll*Mfll » 

aiii— c«cti*ovo 

dv**i.*/(».s.(Ba1(ii.«ati1.ii).«C(Ii/iw;iii 

T£AM2UMNtll*CPfll-OTH*fCI>«*$A|OHOOTO-CACT2l 

i  -osot.noi.cim/deliiii 

temi.facti»dv8 

8  ( 1 1  «TEIM2>*  1 1 1  .]£M1 

cm— tcmi 

0 1 1 1  »1A  111  .TENNc.  (nNAS.OttOOTO-H9AB.  IN0N1 1 1 
1  -NOCII  l/OT— F*CTZ«H6AS'OS0T.*01.m/0EL<II  l«Or»( 

is  Roi«ttT<i«ii 

M-OCIE-PtT/BOl 

CP  1  «CF  T 1 1 . 1 1  •*  1  .CPC  1 1 « 1 1  •  1 1 .  *-A  1 » 

CPI1>1I«CP«I1<1I«1II»1UCPCII«1I*I|.»'III«1II 
IN— A01*f*M*«01*f<l>T 

hbm>pe  te-no  <  l  •  1 1  -oc  1/hho  1 1  ■  «hc  1 1 » 1 1 

T-H01«H1«M1NWII 

28  CIU  LOOK  l2lTtlIi|)il!|ill«,l,M,M«S,CNU<ll 
HGAS*«AS»OCU« 

($EHMX«KNMMW«I1.II 

sotawofainowigwii.il 

IF  111  24.24.2S 

24  EbO-tSMSOItAS 
HM-MG4S 

60  TO  3* 

25  TEMl>KICtIKFWMaT<MI>(F|/liall)l>nil 
C|ll«Cfll.TCM3 

OlH^)IH.TAIl.||»TEOW3.IFACT21WA»(OSOT«OOI»W|/BCHll>»OTH 
TB*!8-W*DS0T*80|!I 
30  CONTINUE 

•fII*OTH/OCL«ll 

II>II.IMTN/MU|l|a«S|| 

C  NON  TK  LAST  AOLATINO  NODE  0C0U10ES  DIFTEOENT  TOCATKNT 
OBOUTO— OMKll«.l/l|Nl|NLl*UCLIM.II 
F4CT1^)TN/ IOCLI«A.I«NNm.ll 
AIM.  I— CACTI  •UNO 

0V*.l.*/t*.5*f*AII«.l**ATl*MI|.*CfaLI/**U*.ll 
ClM.1— r*C71«0*0 

TERN2WUNIW.I  «CP  1(0.1  •ICPOAS^INOOTO.  fOSOT20CLUA.il* 

1  (NO  UO.I*CNIO.  -N01*CPIII*0TN 

OlNLI*TCNIUK<NLI-AIIO.I 

PUQJ*TA  11^1  NTCHNg*0 TIC  |U081(QNOOTP  CIO*  INONIALI- 
1  NOUA.il/OTNHMOTAIOO  UO.IANUO.I-NOIANII/IICLUA.II 

K  »  NL 

OSWOW*  OMOGINU 
ifinonui  ri7N.mi.iir* 

TIT*  FACT2*«SK/IOCLtM.|AMHa.ll 
TENN2*FACT2*CNOAS*01N 
BINLIa*  INII  *TCM2 

D INL I  *0 INL1  •  T  A INL I  ATUMP-r  ACTfANMOAtTN 

CALL  LOOOU.TAINONI.TT1.TI«.*.*.*.HOAS.CMAS.1I 

KAteAALtflMS 

C|NLI*CUA.I-FALI2*CP«AS*0TM 

OlWU*OWO.»-TAINONTAroCMHWAI»OniifOCT»A|M*t*OtN 
T1T1  CONTINUE 

C - NON  FUN  OBOOWIA*  Mtt-Ol  ir  AMT 

IF  UMUI  TXM.T290.T2Sl 
TXSl  00  TX*  L-l.NOOu 
UNflAIH 
LOU*.  I IH 
00  7241  1*U.LU 
IFII-MIlLI'll  T2S.TM.TX9 
TXS  ISM  Mm  ONOEI1I 
UNOSauM 
K*A*1 

ONOOtO.  0*0*111/ 1001  H*OEL  111  I 
r  ACT  l«TN/ IOCL 1 1 1  •••  1 1 1 1 
FACTI-OON/ IOEL II I  AMI  1 1 1 
A(«l— CACT1«P*9 

0*0*1. */l*.».|*ATIIIUUt  |l*ll  lUKIII/MHII 
TEAM*— llll«CNl«l«HeiCO0*O»H— TO  TACTXII 
TEOHI^ACUAONO 
OtKI*TEONX**IOlATEON| 

CIKI— TCMI 

QIK|aT4HIaT*0NX*IN»AS*0N«0T0  iMHHOMm.NOllll/OTN 
C  -CACTXH*WOTM 
IriltCAIHMUII  TM.TXT.TtT 
TX*  CO(|.||ACOVil*ll*AII.|IKXCI|.|)A|I.OHlll*lll 

MOAINO  a  I  •(#  1 1  *  1 1  *00  |L  I /WM IU  «NC  1 1  *  1 1 
CALL  LOOK  IX.TAII.il.  TTI.TNP.*.*.t.MPM.<PPM.I» 
HAA9.IMAT.QCLN* 


CM 

TX1 

CM 

722 

CM 

723 

CM 

T24 

CM 

72S 

CM 

72* 

CM 

727 

CM 

72* 

CM 

729 

CM 

73* 

CM 

731 

CM 
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CM 

733 

CM 

73* 

CM 

T3S 

CM 

73* 

CM 

737 

CM 
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CM 

73* 

CM 

74* 

CM 

741 

CM 

742 

CM 

743 

CM 

7*4 

CM 
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CM 

74* 

CM 

747 

CM 

74* 

CM 

7*9 

CM 

75* 

CM 

7S1 

CM 

7S2 

CM 

753 

CM 

754 

CM 

755 

CM 

75* 

CM 

757 

CM 

750 

CM 

75* 

CM 

7*0 

CM 

7*1 

CM 

7*2 

CM 

7*3 

CM 

7*4 

CM 

7*5 

CM 

7** 

CM 

7*7 

CM 

7*0 

CM 

7*9 

CM 

77* 

CM 

771 

CM 

772 

CM 

773 

CM 

77* 

CM 

775 

CM 

77* 

CM 

777 

CM 

77* 

CM 

rr* 

CM 

/** 

CM 

7*1 

CM 

702 

CM 

703 

CM 

70* 

CM 

705 

CM 

70* 

CM 

7*7 

CM 

70* 

CM 

70* 

CM 

T9* 

CM 

7*1 

CM 

792 

CM 

7*3 

CM 

79* 

CM 

79* 

CM 

79* 

CM 

7*7 

CM 

79* 

CM 

7*9 

CM 

•M 

CM 

Ml 

CM 

M2 

CM 

Ml 

CM 

M* 

CM 

MS 

CM 

M* 

CM 

MT 

CM 

•M 

CM 

M9 

CM 

•I* 

»27 


if  n-*r  I  <Li>it  li*. ll*. ill 
It**}*  l•f<CIMIlWil•lllll 

cmmcmi.icmu 

OMI<OWI<IIIHIFKMM<CIMMMIII 
72*  COLL  SSKTCfMSfSKSkI 
6e  10  1 11*0. 11*1 1 <KKS* 

72*0  MIUW0OT.T2III  l.K.L.OCLU  1  <MII|I  .CPUl  <CPVU1  .CPCill  -Mill  . 

<OMIi  <SAT  <11 -S0<  1 1  .sohUI  <0111  <0M*<  11  <0MOTO<P<L1<PP1LI< 
2HMSi»US2l«l<OiM<MM<i«M$<tM<Heil<r<CT2lT(MI<l(M2ilUnl 
MilU  >IM(<UI<<WI<OMI  ,(MI  <0111 
111 S  f0OMIIi>i]»tl«il<.VIMItllCIO.lll 
72*1  COWflMUC 
72«  Cox  Tim* 

C  MOW  FM  SAC*  UP*  17  MIT 
TIM  IFIMMHM2I  »*<•*<** 

*S  00  SS  1  <<MM2<M<FM 
««1 

7 *C 1 1 «07h/ IOC*  1 1 « «MI  1 1 1 
0<K1  «  <f<CIl<0W 

ov*»  i .  •/ 1  •  .  s<  <ao  n  i  •  <•*  r  <  i  <  1 1 1  <«c  1 1 1  /mu  1 1  • 

C|l|F-rKTl<OM 
K<KIM<OIII<CPI  I  I<C<K1<0I0I 
SO  0<K1<70<1I<<MUI<CPIII 

C  MOM  17  THE  LOST  MM  MS  1KSULOTE0  K E  MIST  SCP0IS  LOST  I  ON  C 

oo  iriHatsi  m.ts.ks 

70  K<KIMMKI<C<KI 
C<K1<0« 

C  WOK  KL70WE  SOUSS  BCOUCTlOK  WE KW1TC  LOST  LUC  OP  NOTMIO 
00  0<KI<IMKI<C<K1<ISCS 
L<* 

00  •#  I<2.K 
L<L-1 

0<LI<O<LI-C<LI/K<L<ll<0<L<ll 
*0  OlLI <K<L1-C 1L1/K <L<11<0<L< II 
K<ll<O<ll/0<ll 
0<IIMl<ll/Olll 
PGPW<LOO<OSC6<1 

PGPUT<P6PUl<t'GPu*OTH/0*70MI<OSU 

OEC<PMOSCM<SOCO<I 

orcowi  <occomi  <oc  ushotp/omo  <  i  i  <osu 

C  <r< 

C  SuKEOLC  SOUMMSt  COMITIOP  POCKOOC 

c 

Cm2 -C< 

IPI<0<I1 

ITL'IO 

1IS<I 

S<CT<<1<<SKELL  1<S0<SMLL<CPEU1 
17111-21  2S02<I*20.2S01 
1*20  0*010-0*0/ 12000.0 
TOIIIwC 
M<0- 
MX0<* 

*••0.0 

0»<0>01K«I  <KOOUI<N0S<11 1<MW0<MC<U<MSI 
•oO-O. 

S»0<*. 

OOCMO. 

CM«<0. 

OO  10  l*J7 

2501  TOOC-VOOOO.O 
OS010-S. 

CS7I<<I<I1 

oocmo. 

OC(P«<0. 

60  TO  2S01 

2502  17  <Mkl  I *2*  <1*2*<1*2S 

I02S  CmkCM/IU.sMT/sOS<KS*1i<MUI 
Cm2<Cm 

102*  GO  TO  II*27<  1*2*1  .MIS 
1*27  *7<<CKfr-CML*  OS  I  AN 
CO  10  1*2* 

1*2*  S7-ICK0— CMLI/Cm 
1*2*  COOT  HOC 

17  <*•<!•. Oil  I<22<I022<I07J 
1*72  CmoOMF I  .-,S<PKI  I 
•0  10  1*2* 

1*2)  CnoCM/lfOPIPNII-l.KPMI 
1*2*  00  III  U*2M<I*201I<MI« 

1*2*0  CMaOMOAl  ICM.OS/OWOll  ITKOUI<TKSIISa*lll 
1*201  COOT  IMK 

(If 

coll  lock  u2<*o  /iCmm-wo. TM<0<S<0<S<12<02<S1 
I<<0<IS<I2I 

IMMI 

I7<I<PS<II  2S1S.A1S.2SIi 
2SII  COLL  LOOKIU<P**S<IPS<S.*<*.*.T2<M<*1 
IPS*  1*11)1 
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ILO 1 1*)  !<*(,.  1P«»1) 

can 

9*1 

I*.  U*)**Ml  1  1*0.  1P»»  1) 

COM 

9*2 

llO  <  1 7 1  »««.o  1 1N6»  1 .  IP**  i  i 

can 

9*3 

IhI(17I*MhI  IlMG.l.lPl.U 

can 

9*4 

can 

9*5 

I7*lLb(16l 

can 

9*6 

lv»lLb(17l 

can 

9*7 

2511 

IlOIU)6«.OIl«iiWl 

can 

9*6 

lHl(l*i*«MlilnG«lpai 

cate 

9*9 

HO  1 151  •**.<!  1  IMG.  1UPSI 

Com 

91* 

1h  1  1  151  **»»1  1  INS*  1 . IP*') 

cart 

911 

11*1LUI1*» 

can 

912 

12»1L0115I 

COM 

913 

ir(t*v-l>  2512.2512.2517 

cm 

914 

»n 

IF  1 l«l 1 171-1*1  *2«.*2«.251* 

can 

915 

IF llrti 1161-13)  *29. *21.2512 

can 

916 

2512 

[Flint  11*1-11)  *21. *21.756* 

can 

917 

IF  llm  (151-121  *21. *21. 2565 

can 

91* 

TABC*TT5(Il»lN&«lPS)*tflir*(TTS<12«  INS.1 .IPSI-TTSt 1 1  .INS.  IPS)  * 

carl 

919 

IF  (NFV-1)  2513.2563.2515 

cm 

92* 

2515 

TABC»iAac.7S»*.(ns(i3.iitG.iPS.ii.*s»iTTSU*.l"&«i.iPs»l>-TTSii3.iNcaM 

921 

16«(M»U)*TA0C) 

cm 

922 

2513 

ir |TS*rt-TA»C)  42l»42l»42l 

cm 

923 

C 

-  ABLATING  SUOfACL 

cm 

924 

*31 

If  1  !*•»  *22. *22 .*23 

cm 

925 

*32 

Ct«L*TLMC  « 1 1 .  IMG.  IP!)  -*0*  <  TUC 1 1 1 . 1«*.  IPSI -ILMC 1 12. INS.  1  .IP«»  1 

cm 

926 

CM>*€*P  100.1  *Oi 

cm 

927 

1aB*1 

cm 

92* 

*23 

C1U  LOOr.a4.CMDL.TLMCII.lM6.lHO.TTSlI.lHS.IPll.TOCMIl.lni.IPSl 

•CBM 

929 

UBPFii.iMS.lP0i#i.Y2iii.r2t*».3» 

CBM 

93* 

in*«i«a*i 

CBM 

931 

call  lookus.cnul.tlmcii .img.i.ips> .ttsu.ing.i.ips> .tocmu.imc-.icsm 

932 

1 

1  .IPS! .T0P7 a .1M6.1.1P1I .1.72171 .721111 *31 

CBM 

933 

isoumusi 

CBM 

934 

IruPtell  *23J. *233. *231 

CBM 

935 

-.231 

CAUL  LOMU6.CMDL.TLMCIltlN6.lPl.lltTTSU.lM6.IPlrll.TOCMIl.lMS.ICaH 

936 

1  PU.ll.TBPFIl.  IMG,  in.  >1.1, 72(13). YZIUI.JI 

CBM 

937 

iac«iaii6t 

CBM 

93* 

call  lomu7.o«il.tumc(i.ins.i.ipi.ii.ttsu.ims.i.ipi.ii.tocm(i. 

ICBM 

939 

U.1NGM.IPM.  11. 1.7211. 1  .72(221  .3) 

CBM 

94* 

160*I«UT| 

CBM 

941 

Do  *232  1*1.12 

CBM 

942 

*233 

72Ui*T2iii.vap*(T2U.i2i-r2(in 

CBM 

943 

IF(NDrF)  *237. *237. *231 

CBM 

944 

*231 

1FIV0P-1.I  *33*. *23*. *2 75 

CBM 

945 

*335 

72131*721151 

CBM 

946 

72191 *72121* 

CBM 

947 

GO  TO  *237 

CBM 

94B 

*33* 

IF  1 VHP)  *336. *237. *237 

CBM 

949 

*336 

72(31  .(72131-721 151*7*0/11. -TIPI 

CBM 

95* 

72(.l»(r2IVI-72l2|l*V8*>t/U.-YlP) 

CBM 

951 

*237 

COMTlNkC 

CBM 

952 

*233 

1)0  *2»  1*1.6 

CBM 

953 

*26 

73111 *72111 *vaM*172II.6l-T21I»l 

CBM 

954 

IF  1721  111  *261 **361 .*261 

CBM 

955 

*26* 

ITL-11S 

CBM 

956 

GO  TO  *356 

CBM 

957 

*261 

CONTINUE 

CBM 

95* 

Call  LOOK i*. 7211)  .TT2I1.ZI  .TEP 11.21 .1.1. 1. EMIV.OMIV.l) 

CBM 

959 

IF (NAIL  1 1 1 —2)  *21. *37 .*36 

CBM 

96* 

*26 

CALL  LOOM  13. 73111. TT2tl.ll.TCPU.il  .1.1.1.73. 03.11 

CBM 

961 

CMI7*tMI».APl*173-CMI7) 

CBM 

962 

OMI7*OM17.AP1*1D3-OMI9I 

CBM 

963 

*37 

T5S1**t2I1  1*72  111 

CBM 

964 

TW2U1 

CBM 

965 

1AD-5 1G»EK 1 7*1  S50*TSS0*vF 

CBM 

966 

*16 

ERf>*Cn*Y?(  2 1  .CMI  7*01A-lAI>-S  1 1 »  *rs.  E1F* 

CBM 

967 

UCBP*CM*  r  2 1 5 1  *  1 1 OftA-lAD/EM 1 7 1  *0M  1  *-*  .  /  T  S*«A£)-B  1 1 1 1  *72 1  *  1 

CBM 

96* 

CBM 

969 

7ITCBU1SI-CM0L 

CBM 

97* 

CIT£»tITS)*C«l« 

CBM 

971 

C»«L*CMDL-tBaC 

CBM 

972 

Ct#tl  •—!.£. 31 

CBM 

973 

CmmUU.I.C.JO 

CBM 

974 

IF  1 1 Lot  1*1— I RA)  *361. *363. *363 

CBM 

97S 

*361 1 

|  IFIIL0(15>-1RS>  *363.* 3*3. *363 

CBM 

9/6 

*3621 

!lFl*S»H-l)  *271. *371.4271 

CBM 

977 

*371 

IF  1  ILul  161  — IRC  1  *272^363.4363 

CBM 

97* 

*272 

IF  IILOI171-MDI  *273.6363.6363 

CBM 

979 

4?7T  CMHl *AMAA I  ULMC  ( IAO»  !*«»«  1  •  IP0»  1 )  «TUC <  K0»*l  •  1  •  1PW«  I  >  , TU«C  ( ISC 

•CBM 

9*6 

1 IMG • iPfl* 1 » *TLNC 4 1BC~ I • 1n6# l*ft« 1 M 

CBM 

9*1 

c 

W<1*B91  fCMLT*)kPl#WU«l-A*l  (CMLT«HV»I*0M1,AAI  <CMLT*lMMCIFlXAH-lNC*M 

9*2 

427* 

CMMl*AMUX  ICMHI  «TLMC(  ISU*1N6«1P«)  *TUCf  !«A*1  «TUC  4  IAS»  IM6C8M 

9*3 

CBM 

9*4 

C*OC*AMJUl  {CHOL.CIWIJ 

CBM 

9*S 

4363 

IF  « irtl  (1M-1W-U  4366.4366*4364 

CBM 

9*6 

4364 

IFII*I <1S>~1R»~1>  4366*4366*4274 

CBM 

9*7 

4274 

lFtNPft-1)  4275*4275*4276 

CBM 

9** 

4276 

XFCXHI  (Ifei-IPC-n  4366*4366*4277 

CBM 

9*9 

4277 

IF 4 (Hl-4 ] 7| )  4366*4366*4276 

CBM 

996 
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*27«  C—l».*m1p,|  ,  IlmC  (  1*0. 1  .I-O-l  .  ISS-I  >  .TLMCdtKW.INO.  1,1P«.|  I  ,  CB« 

CBM 

C  (**-!  (CSM. T.lfirM. 6*1.2. AA1  ICML!>>0J>1.6«!.1«ARI  (CHLT  .AAMC  (F  lNlN-ANCBN 

*27S  C*MA.»<«ll,l  (C-KLA,  1LNC  I  lfi».  ;  ,  JMG,  1PVI  -TLNC  <  IMA-2.  ING.IPU)  ,  TLMC  (lBS'lCttM 

I.IMG.i.lPe).TLML(I90.2.!Ml..|.IPfl|)/2.  CBM 

CmOl*AM(N|  iClUl.CMU)  CBM 

|FIII>-ITl-1i  *3o6.»3Sl.*3S2  CBM 

*3Sl  £B«S»t.-*a  CBM 

C.OL-CNK4  CBM 

•jO  10  *3*7  CBM 

*JS2  If  l£#«»£resi  *3>-  .*36  7 .*3t»3  CBM 

*3S3  C«Ol*C«»i  CBM 

uO  70  *3*7  CBM 

*js*  itl»ss  cbm 

IF«ta»Cl  *  JS5.A36  7  .*367  CBM 

*3SS  CmOl*CNNJ  CBN 

UO  10  *3*7  CBM 

*3**  IF1II>-ITL)  *J67,*3S6«*367  CBM 

*3S6  C<OC>*Mls)  <Tlhcui.1nG.IS«|)  .TLHCd2.lNG-i.IPSi)  CBM 

*367. *367, *279  CBM 

C  I  j  J.kPl.l-GMi,*!  (C»9_T.n*«PI.&Hl.3I  <C3B.T.L0NCIF1NIN*CDCBH 

*27*  CM0LKAMlNI(CM0L*lLMC<I3.lM6.fPa*lt.TLMCfl6.1N6«l.IPfl»l)l  CBM 

*3*7  CMO«t  if*  I CMOL )  *0*-C»6.  CBM 

If  « US-S9I  **0.**0.9V«  CBM 

**•  ITS*IISM  CBM 

IF  ll»ICM  1-3.  I  *372* *372 **23  CBM 

C  - MOM*  ABC  AT  INC  SURFACE  CBM 

*29  tS*TSA«f  CBM 

lAB*«  CBM 

O*>*0.9  CBM 

*39  If  1 11-31  *392.633.6392  CBM 

*392  1L0(I6)«I  CBM 

Ih3i iai*«Hl ilMC.lPRi  CBM 

lLOll*l*l  CBM 

IhI ( 1*1 **Ml I IMC»1 . IPfll  CBM 

CALL  LOOM  U*.I>.MSIl.IM*.lFNI.TC>CMII.lKC.IP*).».».».r2ll).T2(2)C8M 
I. II  CBM 

CAL*  LOOK  (19.f».7TSd.  IMG*  1  •  1PRI  •  IOCHI I  •  INC*  |  .  IPftl  ,9,9.9*Y2(3I«  CBM 
1T2I*I*1>  CBM 

*322* *322**39 3  CBM 

*39  3  IlO  129)  *  1  CBM 

IHI2»l*fnIIIll6.1M*ll  CBM 

ILOI21IM  CBM 

1~1I2I<«*hKIhG -l.lPOdl  CBM 

CALL  look  C29.  IS.  ITS!  1  *1NC.  IPR-ll  .TOEMI1  .INC*  (PR-11  ••*•*•*  *l’SI  •  CBM 
1*216). II  CBM 

CALL  look  I2I.IS.T7S II* ING-I.IOR-U  *TC»CHI 1 .lMB-l .IPH-11.9.9.9.  CBN 
1*2(71. *2(91.11  CBM 

CO  *321  1*1.*  CBM 

*321  *2(11**2(11 -**P*(*2( l-*l— *211)1  CBM 

*322  00  *323  1*1.2  CBM 

*323  *2111**2(1) .»«M«(T2(I*2l-T2(Ill  CBM 

*33  CALL  LOOK  I*. TS.TT2II. 2). IEP(1. 2). 0.9. 9. EMlV.OMIV.il  CBM 

I.*  (MAIL  (11-21  *3*.**2.*3*  CBM 

*3-  CULL  LOOK  (3.TS.T*2<1.1| .  T£P  ( I . 1 1 .9 .9.9. *3. 03.1 >  CBM 

Em1**lm1*.aP1*  (  V3-Cm]VI  COM 

0"I»*0Ml*.*f*l*(t/3-0MlVI  CO* 

**2  3SSO*TS*IS  CSM 

If (11-31  **22.*62i.6*22  CBM 

**21  I*  (EmIV)  **?:.**23.**2?  CBM 

**23  E«I V* I .0  CBM 

**22  OA0*SI(.KfMl**TS»>KTSSB6vF  CBM  ; 

*39  E(W**Cmk*2<  1 1  *EM1  V*ORA— BaB—B  ( l)  *T5.EBFJl  CSM 

0EM)«CM*T2(2I  *  (  (<I«A-MAO/CMIV(kOMIV-A./TS*IIAO-(II  1)  >  CBM 

E!W)C*E»f(20t(«(  CBM 

V1TE9 1 1 ISl *TS  CBM  ! 

E!TEO(I*S|KCew  CBM  ; 

tskIS-c«oc  cbm  ; 

IF(II-3)  *391 »*S*7 ,*S97  CBM  ] 

*391  *SMl*-l.E-3«  CBM  ] 

T  SNA* • | . £  »  30  CSM  ] 

IBA*I«(10I  CBM  I 

ibb*ik(19)  csm  i 

19C*I'.  (201  CBM  ; 

1*0*I»(21)  CBM  I 

If (IL0II8I-1BA)  *S99.*S91.*S91  CBM  | 

*S00  lrilLO<19l-IS0)  *S92.*Sol.*S«l  CBM  | 

*S02  IFlWV-ll  *S03.*S93«*S9*  CBM  I 

*SC*  IF (1LO(20I -INC)  *S9S.*S9l.*S01  CBM  I 

*S0S  If ( 1 LO ( 21 1  — 1MO)  *S96. ASg 1 .*59 i  CBM  | 

*S96  T  LM I  *  AML  1 1  ( *  T  S  ( 1*0  .  IMG*  )  .  1 3R*  l I -TTS ( 190- 1  * IMG*  1  •  IPO  •  1 )  .  TTS  ( IOC  •  CBM  I 
IImO,Ia<-1i.TTS(IkC-i.1MG.1W*..||)  CBM  I 

*S0  3  1  SMl  *AMAl  |  (  TSKl  .  1 1S  ( IBS.  IMG.  1  .  1PR I  .775  ( IRD— I  .  IMS*  1  .IPO)  *  CO*  1 

ITTSdOA.lMu.lPfll.TTSdRA-l.lMO.IOBD/Z.  COM  | 

TS"AMA*| (TS.TSMI)  /•  CBM  I 

-SOI  IFdril  (IBI-IOA-Il  *S97.*S9 7.AS99  CBM  1 

*S9»  IFIIMI 1191-IUB-l)  *S07.*S97.AS99  CBM  1 

*S5*  JF(*fM.-l»  *SI().*S19»6Sll  CBM  1 

*S1I  "IFlIxl  (20I-IBC-1)  *S97.*59  7.*S12  CBM  1 

*S  1  2  I*  (!•*!  (2K-1U0-II  *S97.*S07.*SI3  CBM  1 


«5IJ  TSMA*»*lNl(tlSClhOM.lHr,»l.lf,!<«lt«TfS(lWWrIHGM.IP«»l>, 

msu-c-i,iH(,,iKK.i)«irs(irfc-z.  img.ipr.im 
*5|0  TSMA=AMIN| itsma.i  is<iRt>.i .  ihc«i  .  :pr>*  risi  iph«2.  ino.i  .ipri  . 

•  m$(lKA.l.lHG.lP.0..1TS(|>.A.2.IK^.lP>il)/2. 

rs-*MiNi  us.  rb««i  t 

*so/  comjmuc 

irurs-soi  **i.**;.99s 
**1  !TS*1TS.) 

lF(AUS(ERft  1-1.  )  *390. *390. *30 

998  *«ITE  (KOUT .  52vl 

»RI T£ (KOUT  .582) (<1T£R(I) .LITER) 1) .1*1.51) 

•RlT£|K0uT.582>  TM.DTM.vRk.£RfA.TABC.EHIV.D»<lv.RAO.Mfl>.CH.fMI, 
1  D(II.ME.fPl.<W».T2U).»2(2),72(3F.T2I*> 

ThsTmFIn 
60  16  3000 

C  - POS1  ITERATION 

*390  lr(ll-2l  *371. *3/1. 1*37 
*371  72121=72(1) 

*372  OCHEH*72(?) 

IE (IStNl 1)1*373. *37*. *373 

*37*  GC 06**0. 

GO  70  1*39 

*373  Call,  ogle u.is.ocoNv.iSrN) ipr> .ttsen( i .ipri .tkseni i.ipri .tcpsen 
llPflii 

IECNPk-1)  1*39.1*39. 1*3A 

1*38  CALL  UGLEU.TS.OO  •  !  SfM  ( 1PR*  1 1  . T TSE* ( l . lPfi.J )  .THSEMIl.IfM.il. 

ISENll.IPH.l))  *  - 

OCONV=OCONV.VkP*IOO-OCO*VI 

1*39  OcrCM^(OCMEM»OCUMV)*Ol  t 

HhsOCOMV  i 

OCOMV*CH*  |h£-OCW)  _  </ 

050TB*CHG/RhO(2) 

IF(NBPF)  1*37. 1*37. *375 
*375  IrilAu)  1*37.1*37. *376 
*376  OTL'EAPIrZOIl'CM^ORI 
1*37  ROI1)*BON(1) 
bR*CM/CHZ 

OPP=EMlV*0«A  ' 

bC0n0=-0 1 1 ) -a 1 11 *7 A 1 1 1 
GCOfc*  T*OCOnV  T  .OCOR7»DTH/Ak£A  1 1 »  *ASO 
VCHf  H  |  *OCME»l  I  »0<XA1*»)TH/ AREA  1 1 »  •  ASO 
OCOfcO  1  «OCOMOT  »0C0*K>*0  7H/  ARE  A  1 1  >  •  A  SO 
UPPT=bRPT.ORP*0TH/AREA<l  |*A5U 
haOT«m AOT .MAO*03n/AREA ( 1 1 *ASU 
C  '■ 

OEOT=kOM ( 1 1 *CP ( 1 ) * ( TS-TsavE > *OEL < 1 1 /07n 
DO  95  1*2. NL 
POlII  •  B  0*1  II 

-r*p  *10(1  1-All  )*TA(1— 1 ) l/Bll  I 
EOT  *oeor-RON  1 1 1  *CP  1 J  )  •  «  TE77P-TA  1 1 1  1  *0£L  1 1 1  »RR  C 1 1  /DIM 
7a 1 1 1 wTEMP 

IF  19UWM  »*ax I  97.96.96 
06  <*ML  .•* 

IAIMBM-D  •  lAiNL)'  •= 

DO  96  1  *  (KIM.NuMN 
K  *  <♦  1 

9«  T  A  I  I  )  *  IO(M*AIK)«U(Ul)r/B<Kl 
97  UtCTT*O£OTT*U€0f*0TH/A*EA(l)*ASU 
I F  C wDtMJ)  7S5'7S6*7S5 

*55 

Ll#*NLA  (NO0U) 

00  757  1«4*L«LU 
7S7  ftOUlzROMC!) 

7Sfc  continue:. 


SHfllNA  (ANO  DfrOPl  OF  LAST  ABLATING  N00€ 

d£linl>*ocl<*-m-4js  isce  int  occompi 

1 F I OfcL I NL ) -OELH) l 49 • 1 *9 . I SO 
149  0«iLP«O€l_<NL)  -kO(K.)  ) 

0*|_CP»0«LP*CP  <Nc> 

NL=NL-1 

»C<NLI *«C<NL*1> 

W*L*0£L (NL)*WO INL) ®R« {NL ) 

o»lc*o«l»cpinl» 

HAPhB=0«L»M C NL I *0«LP*H I M.  • I) 

T0FI*O«L*[)KLP 

TOP2*0ttLC~O*t.CP 

T0P3*D«LC»TA(NL)  •Of»LCP*TA(HL*l> 

V0L=0£L  (M_>  **«  (ML)  ♦0£L  (NL*  1 )  (ML*  I  > 

OEL  (NLl  =C€C  *  NL)  ♦DLL  (WL*  ]  ) 

fiA<NLI*«A(NLJ  -6.*OCL  I  **.•!! 

»OINL)*TOPl/^OL 

CP CNL»*T0PZ/T0P1 

UlNLI*T0P3/T0P2 

N  <  NL  I  *NAPmB/  T  OP  1 

OCLR-OtL (NL* 1 ; •**  (NL* IT /VOL 

C7-l*0-O£L« 

oz=o.o 


CBM  1081 
CBM  108? 
CBM  1083 
CBM  1084 
CBM  108S 
CBM  1086 
C6M  1087 
CBM  1088 
CBM"  1089 
CBM  1090 
'  CBM  1091 
CBM  1092 
CBM  1093 
CBM  1094 
CBM  109S 
CBM  1096 
CBM  1097 
CBM  199C 
CBM  1099 
\  CBM  1100 
C8M  HOI 
CBM  1102 
Kl *C8M  1103 
CBM  1104 
CBM  110S 
T CPC BN  1106 
CBM  HOT 
CBM  1108 
CBM  1109 
CBM  1110 
CBM  fill 
CBM  1112 
CBM  1113 
CBM  |1|4 
C0*  1115 
CBMv11|6 
CBM  NJ 117 
CBM  1118 
CBM  1119 
CBM  1120 
CBM  1 121 
CBM  1122 
CBM  1123 
CBM  1124 
CBM  1125 
CBM  1126 
CBM  1127 
CBM  1128 
CBM  1129 
CBM  1130 
CBM  1131 
CBM  1132 
CBM  1133 
CBM  Ip* 
CBM  11JS 
CBM  1136 
CBM  1137 
CBM  1138 
CBM  1139 
CBM  1140 
CBM  1141 
CBM  1142 
CBM  1143 
CBM  1144 
CBM  H4S 
CBM  1146 
CBM  1147 
CBM  1148 
CBM  1149 
CBM  1150 
CBM  US! 
CBM  1152 
.  CBM  1153 
CBM  1  1S4 
CBM  US5 
CBM  11S6 
CSM  U57 
CBM  1158 
CBM  US9 
CBM  1160 
CBM  1161 
CBM  1162 
CBM  1163 
CBM  1164 
CBM  1165 
CBN  1166 
CBM  1167 
CBM  1168 
CB*  1169 
CBX  L170- 
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Gf-C2  4 

n-nI-jtm 

«*N 

f  2»*J€lW 

a -w 

GO  10  1TV 

07-02*1. •  . 

173  f*«02-C? 

ITS*  1 7«» 1 7S 

17*  67-Oti.fi 
,  17S  RaKO 

CZ-CZ.62 
176  tZ^CZ-O/ 

If <rZ) 17S.1TT.1/T 
>177  fiOAIN) -ROA IN) *FZ*ROA  IKI 
ROM  INI -MOM INI *F  !•« OS IK > 
ROC  INI  -ROt  INI  •r^NIOtlKl 
^irW-fcZJ  171.1S6.171 
|7I  N-N-l 

179  M0AIN>-R0aUI«C* 
ft0B(N|««0SIK»*C2 
aOCMlMWCIKIHI 
GO  TO  172 

176  ROA|N)-«Oa|N»*NuAUI«*2 
MOO  I N1  -4*06  (N»  -MB  I K I  «02 
ROC  INI -AOC  INI -IhlC  IK  I  ««Z 
GO  TO  173 

C  . 

iSO  GO  TO  416  , 

1  fit TORN 
fe.NO 


l 


\ 


mu  •  Fin  •• 


\ 
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8fiS8Si8«88S8IS88iSSSi!8li8S888SS 
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SUBROUTINE  INPOuT  iNPdO  1 

COMMON  KOUT.IEX.UEN.VR  ..  1NPIHJ  2 

COMMON  Ini  <381  .ILCXOBI  ,<R<30)  .TT2I30.10)  .TCP  1 30.  lot  .TKP13D.10)  .TMZINPOU  3 
1 130. 1 V) «1EP I3Q. ID) .TTHlIO). THE  130). TOR (30). ICHI301 .Til  130) .TNG  1301 INPOU  4 
2.0H1212) .HEC0R0I36) ,SOl?0)  INPOU  S 

COMMON  RHOUOI  'INPOU  6 

COMMON  HATEl  SO). OEM  S]  I  .TA(  SI  I  .H  (SOI  .RC  (SOI  .BA  (  SO).  INPOU  7 

1  AREA (  SOl.bMAi  SOI.HAV)  SO)  INPOU  B 

COMMON  ROA<  500).ROB<  SoO) .ROC (  S00)  INPOU  9 

COMMON  TPr ( 3) .NM6 (3) . THg t IS . 3) . NLU ( 15.3) . NH 1 ( IS. 3) .  InPOU  10 

IKHI I  IS* 3) .TTSEN12S.3) .ThSEN (25. 3) . TCPSENI25.3) .TLNCI2S.1S.3) •  .  INPOU  11 

2ISEN(3> >TP1 1301 .TTSI2S.15. 3) .TCh£m<25. 15.3) .VFZ. CMH  INPOU  12 

COMMON  TUPM2S.1S.3)  INPOU  13 

COMMON  NPR  INPOU  14 

COMMON  LCT .NPO. 1 I .NflM.NuMNiNL .DELHG.DELN.RFT .RHORA.RHORB.RhORC.TRAINPOU  15 
1CA, TRACB, TRACC. HMOOA.RHOOB.RMOOC.EA. EB.EC.BA.BB.8C.ps] a. PS1B.PSIC. INPOU  16 
2 TR ACM. PET .PETE.K5V.ETA.OTPR3.DTPR2.DTPRT . TPR3.TPR2 .THZRO.TMTIN.NT .  INPOU  IT 
3TNNT .oAMA .OMG.NU  >F JFH.F  JFS.  JF  •  JFHP.  JFNt INPUT •  OTMIN.BRP.NCONV. INPOU  18 

4EPSM.IRES  /  INPOU  19 

COMMON  1NCH.OTH0  /  INPOU  20 

COMMON  NN.NI.NOl  /  INPOU' 21 

COMMON  CHCHI.PYCRI  '  INPOU  22 

COMMON  TBRPI30)  INPOU  23 

COMMON  NR  INPOU  2* 

COMMON  TX<30.6).Fl(30.6l  .F2I30.6)  INPOU  25 

COMMON  NCOM  INPOU  26 

COMMON  N6PF.NFIS  INPOU  27 

COMMON  BREX.SVELL  INPOU  2« 

COMMON  BBBIS.3I .EEIS.3) ,FF (S.3) .PSI <5.31 .RMOO (5.3) .RHOR<5.3) •  INPOU  29 

1NOCOM<S0.3) ,DHC<S) .0MV(5>.RnOC<5).RH0V<S) .PIS) .PP <5> .TREF <S) ,8AIS) INPOU  3D 
2.0MGA<SI .NFI (S) .NLA (5) .TTSI 3D. 10) .TENT (39. 10) .TABU <30.101 ,TCBU <30. INPOU  31 
31«)«A<SO),NOBU.NBM2  INPOU  32 

COMMON  TRACIS. 31  INPOU  33 

COMMON  NBUFTIS)  INPOU  34 

COMMON/O AT A/ASTER. BLANK  INPOU  3S 

DIMENSION  TZSENI2S.3)  .TcZSEN  I2S.3)  .7 SURF  <2S) .  TSEN  <2S)  .  12 12S)  INPOU  36 

OINENSlON  I OPT <301  INPOU  37 

DIMENSION  KSVI5) .KMTLI6)  INPOU  38 

EQUIVALENCE  <DM1, DM1211) I . (DH2.0M12I2) )• ITS. TA)  INPOU  39 

*9?  FORMAT!  12. 3F10. SI  /  .  INPOU  AO 

A 99  FORMAT < 12A6)  >  INPOU  A1 

502  FORMAT  (6X12A6I  ,  INPOU  A2 

503  FORMAT  I//2AX31H - REACTION  KINETIC  EQUATION — ‘-/lH  )'  -  _  InPOU  A3 

SOA  FORMAT  UGaoTHORHO/OMMC  »  GAMMA''  !  BA«EAP  l-EA/T >  RHOOA  <  (RHOA-RHORA) / INPOU  AA 

1 RHOOA ) **PS I A  1/  21XS6H.  GAMMA  I  SB*EXP <-£b/T > RHOOB ( (RHOB-RMORB)/ INPOU  AS 

2RmOO0)««PSIB  )/  I9X58N* ( 1 -GAMMA ) <  BC*EXP l-EC/T ) RHOOC < IRHOC-RMORC )/ INPOU  A6 
3R«00C)*«PS1C  )l  I  INPOU  AT 

505  FORMAT  I/2AX32M - REACTION  KINETIC  CONSTANTS— —/1M  )  INPOU  A8 

506  FORMAT  <llX8MHEACTION2|*MRMO05XAMRM0H8XIMe7X3MPSI8XlME6X6MT  RE AC/2 INPOU  A9 

1 JX 1 OM ILB/CU  FT ) 6X7H 1 1/SFC) 1 2X7H  <D€G  R) 3A7H IDCG  R) I  INPOU  50 

507  FORMAT  I 1AXA1 ,2X2F9.2.2xE 10.A.F7.2.2XE10.A.F8.0)  INPOU  51 

510  FORMAT I12X31HMES1N  VOLUME  FRACTION.  GAMMA  •  FS. 3. 1 7h (MASS  FR ACT  I ON INPOU  52 

1  =  F5.3.1MI/1H  I  INPOU  S3 

511  FORMAT  (2AX32M - TIME  INCREMENT  INFORMATION - /1H  I  INPOU  5A 

512  FORMAT  IbXlSMlNlTIAL  TIME  (SEC) F7. 3.26X1 6MF INAL  TIME  ISCCIF7.2I  INPOU  55 

513  FORMAT  (1M  /6X 1 TnOUTPUT  INTERVAL  >F6.3.1X27nSEC  FROM  INITIAL  TIME  INPOU  56 

1UNTIL  FT. 3. AM  StC)  INPOU  57 

51A  FORMAT  I6X1TM0UTPUT  INTERVAL  -F6.3. 1 A8MSEC  FRONT 7. 3. 1 X9MSEC  UNT ILF INPOU  58 
17. 3. am  SEC!  "  INPOU  S9 

515  FORMAT  (6X17MOOTPUT  INTFRVAL  »F6. 3»  1X8HSEC  FR0^7.3.  IX2DMSEC  UNTILINPOU  60 

1  FINAL  TInE/Ih  )  INPOU  61 

516  FORMA)  ( 6X • 1 9MMAX 1 MUM  TIME  STEP  «.F A.2.8M  SECONDS)  INPOU  62. 

SIT  FORMAT  I/29X16M - NODAL  OATA - /1H  )  INPOU  63 

518  FORMAT (6X7AHNOOE  KATL  TEMPERATURE  RELATIVE  THICKNESS  NOOAL  OC INPOU  6A 

1PTM  CONT. RESISTANCE)  INPOU  65 

519  FORMAT  I7X73MNO.  NO.  (OEG.RANKINE)  AREA  (INCHES)  I  INCHES)  INPOU  66 

1  (SoFT-S— OEG/BTUI )  INPOU  67 

520  FORMAT (3X2I6»F12.2.£13«4<F9.5.F12.6.A1.£1S.A>  INPOU  68 

521  FORMA)  <  1AXA7HM1NIMUM  THICKNESS  OF  LAST  ABLATOR  NODE  < INCHES  I F INPOU  69 

17.A/lAX.10rlTMEHt  ARE  .I?.ADH  NOOELETS  ASSIGMEO  TO  EJkCM  ABLATING  NOINPOU  70 
2UE)  '  INPOU  71 

522  FORMA) (/18XA8M - MEAT  OF  FORMATION  OF  MATERIAL  CONSTITUENTS - /3TXINPOU  72 

18M IBTu/LB) /2 1 XThPLASTIC) IAAmCHAR I 7X3H6AS)  INPOU  73 

523  FORMAT  <20XF9.2,7XF9.2,)IXF9.2I  „  INPOU  7A 

524  FORMAT (/7X2SMENIHALPY  DATUM  TEMPERATURE  PF9.3.1X1 1MOEG  RANK INE I  INPOU  75 

525  FORHAI  I//22X36M - MATERIAL  THERMAL  PROPERTY  DATA - //6X14HMA TER INPOU  76 

1IAL  Nu.  1 1  OX 14HMATER I AL  NO.  210X26MMATERIAL  NOS.  3  THROUGH  10/6X1A1NPOU  77 
2HVIRGIN  PLASTIC15X4MCHAR23X7HBACK-UP)  /  INPOU  78 

526  FORMAT I/6X12HMATERIAL  NO. 12. 30X9HOENSITY  *F8 . 3. 1 X8MLB/CU  FT/  INPOU  79 

l  7x1 1HTEMPERATURE5X1 3HSPECIFIC  HEAT5X12HCONOUCTIV1TYSX8HSENINPOU  80 

2SIBLEAX ) OmEHISSI V I T Y/58X8HENTHALPY/9X7M (OEG  Rt 7X12H (BTU/L8-OEG) AX  1 INPOU  81 
36H (BTu/FT-sEC— OEg) 3X8M IflTU/LB) /  INPOU  82 

4  ( BxF B . 2*0 AF 7.4.9 xF 1 0.7 • 7XF9.2.4XF7 .4 ) )  /  INPOU  83 

5260  FORMA (I/6X12HMATERIAL  NO. 12.30X9HOENSITY  *FB.3» 1X8HL8/CU  FT/  INPOU  84 

1  7X11MTEMPEHATURESX13HSPECIFIC  HEATSX 12HCONOUCT I V I T Y5X8HSEN1 NPOU  8S 

2S1BEE  /S8X8MENTHAEPY/9X7MIDEG  R ) 7X 12H<8TU/L8-OEO> 4X I InPOU  86 

36MIBTU/FT-SEC-0EG) 3A8MI8TU/LB)/  INPOU  87 

4  C8XF6.2.BXF7.4.9XF I0.7.7XF9.2  ))  INPOU  68 

527  F OMHAT  </6A 12HMATER IAL  NO. I2.30X9HDENS1TY  *F8.3* 1X8HL8/CU  FT/  INPOU  89 

1  7X1 1MTEMPERATURF5X 1 3HSPEC IF  I C  HEXTSX 1 2MC0N0UCT I VITY/9X7H (DINPOU  90 


4-17 


«6  m  r<i2Hiaiu/va-ocki«ii*HiB?u/ri-$ec-oe6)/  iimi  *1 

3  iiiri.^un.t.wii.Tti  ■  imi  n 

531  FORMAT  c |H  /6A19HlEMPE0iTuat  IOCS  RISFII.2I  IIMU  *3 

5«  FORMAT  1//2WA6TH— OESIN  DECOMPOSITION  US  SENSIBLE  ENTHALPY - 1 IMMU  96 

S33  FOMNAI  IUIMMIUIM  liTIMJISril.21  INPOU  n 

S3*  FONNAI  I1H  //23A4#H— TINE  uepcnoent  MUWMI  COIN) I T IONS /IN  )  IMU  N 

535  FOMNAI  m.imi».U.*HMM«.M.«MtCDYt)H.l».«<MUIIOII.U.4IACHllflU  W 

ISx.0HPMESSuaE.JA.?HeL0M|H6/91.SMtsCCI.7A.6MPTN.3A.0HENTHALPV.3A,  HOW  90 
nmoi  MK.u.3icotFF.iwl««wcTia«/2UiMiini/iitialiiH(inwsiinu  ** 
30  ft-.ia.iohclo/»o  f i-.ya.smiaini .aa.t^araneter  /lox.tksecomd.  inpouka 
««x.mscconoi>  - —  INPOUIOI 

536  FOMMAT  IM.FI.2.U.I2.«l.2IFI.2.3»>ri.«i3I>n.5.UJI,)l  INPOUK2 

531  FOMNAI  I  IN  /U.IMCH/CHO  ■  PNI/CEaPCPhII-I.I  (ACRE  PNl  ■  2.MWM  INPOUI03 

IOOI/CmO.  8MP  IN  IA0LE)  IIT0U164 

530  FOMNAI  I//2TA30N - SURFACE  EOUILIOMIUN  DATA - 1  1NPOU1A5 

551  F0RNAfllNll#X6SHAER01lCAN  CHARM  I N6  MATERIAL  THERMAL  RESPONSE  AMO  AINP0U166 

10LAI10N  MMWAA/73U)OAGEIVlN  6TX2A6I  INPOUltT 

552  FOMmai  l«l)«HI|it;fU>M»NlilliMMFACEiUiniSMFACE/UiSN(KC)i  INPOU1M 

1TI.4HuPIN,SX'4MIEI».SI.9MR£CESSI0H/2#A.Th|0C6  Ml  .4A.6HRATE/3RX.  1ITOU109 
21#MIN1LS/SCC>>  X*  IIMIIK 

553  FOMNAI I0A24HMN1IIAL  {NTEHNAL  RADIUS. IX.F6. 3. 4X.21HAMCA  PMOP.TO  RAUWOUIU 

10IUS«*f 6.21  / '  INP0U1I2 

556  FORMAT  C0A?6M*INI  1 1AL  EXTERNAL  RADIUS.1X.F6.3.4A.21HAREA  PMOP.TO  RAHW0U113 
10IUS»*f4.2l  II#0MII6 

555  FOMNAI  CRXISMPPLAMAR  SURFACE!  INPOOI15 

556  FORMAT  I9>.4MTHN..BA.4HPHOB.5A.4M91EM.SX.9*»AOIATIOM/9X.SH1SCC1.  1NP0U116 

I7I.*M0P1N.4I,6MFACT0A.4I.9MHCAI  RATE/3BX.  1 1HCRTU/SO  FT-/60X,  INP0U1I7 

2TnSEC0NDI  I  HPWUi 

560  roMNAICl2.SE10.01  INPOUUO 

561  FoMNAT  CI2.F0.4.611)  IMP0UI20 

563  FOMNAT  112.13. 1 1.12. 12.TP  10.S/0F10.S1  1NP0U121 

566  FOMNAI  IA1.9X.2F IX.S.E  I0.3.F10.S.C10.3.FI0.S/A1  .OX.2F30.S.E10.3*'  IIW0U122 

1F1A. 5.  EI0.3.F10.S/AI.91.2F10.5.E10.3.F10.S.E10.3.F10.5.0X.I1I  II0>0U123 

566#  FOMNAI  IAl,9x2F|*.S.El#.3.Fl«.5.E10.3.ri0.51  INP0U126 

571  FOMNAT  I  I2.F  10.5.F  10.S.F I  ••0.F30.SI  1IW0U12S 

575  FOMNAT  1 1 1  .F9.S.  FF  10.S/OF  I0.SI  1IT0U126 

577  FOMNAT  C 1 1.F9.5.3F  10.SI  IMPOU127 

579  FOMNAI  II2X2EI2.5.2CI3.6.1II  IRP0U12R 

579#  FOMNAT  I6X.26MN0  RAOIUS  CORRECTION  ON  CM)  ICWQU129 

5791  FOMNAI  1 3F W.S.F9.6 tFS. 3.2F 9.3. 12.A01  HPW1M 

5792  F0MNAII//6X.3HP  ■•F9.4.6H  ATM//BX.3I25MTENPERATURE  EOSE  ENTM  3SW0U131 

16X.3I25M  COCO  Ml  AT  T-NAU.  II  IRP0UI12 

5793  FOMNAI  C//6A.3TI6XA0  SURFACE  EOUILIOMIUN  TAOEJE  OF  TYPE.  121  HMW333 

5796  FORMAT  C//6X.76MC0UAL  NASS  ANO  MEAT  TRANSFER  COEFFICIENTS  AND  E0UAUW0U136 

1L  DIFFUSION  CaCrriCIEMTs/6X.29NN0NUMt  SURFACE  VIE*  FACTOR  «.  INPOU135 
2FA.3I  ICW0UI36 

5795  FORMAT (//6X.16HN-OOT~OA5/CN'6.F7.6*OX.IOHPRCSSURC  -.F9.4.6M  ATM//  IITOU137 

1 7X.2C4HTEMP.SX.26HM-00T.-  CMEN.PHOO  SURFACE. 3XI/6X.2I36NI0CS  Ml  INP0U130 
2CHAR/CM  (tlU/UM  SPECIES. 2X1 /(5A.F0.2.2X.F7.4.2X.F0.2.6X.  AO.  IX.  INP0U1J9 
3Ffl.2.2A.F7.C.2X»F0.2.6X.A6l I  INP0U160 

5796  format i2fi«.o,f9.o.ii. jc9x.hi .fio.o)  inpouioi 

5797  FOMNAI  C//6A.45HMATI0  OF  NASS  TO  MEAT  TRANSFER  COEFFICIENTS  -.F6.VIIW0U142 

I  6X.2WHUTKUUAL  DIFFUSION  EXPONENT  •'.F6.3/6X.29MN0N1NAL  SURFACE  VIEIIT0U14J 
2*  FACTOR  ..F6.3.11H  I  OPTION  III  UW0U166 

579*  FORMAT  C0X.F9.2»6X.F9.2.3X.F9.2.4X.F9.INP0U165 

22.3A.F9.2.4X.F9.2I  INP0UI66 

5799  FOMMAI  C6x.66MMtAT  TRANSFER  COEFFICIENT  MULTIPLIED  ST  IR  INIT1AL/RIIW0U16T 
1  CURRENT  I  .•EX./6A.63NNHERE  EX  *  .I.O-.2NI/C1.-N)  AND  N  IS  THE  0URNIIWOU16* 
2 IMG  rate  EXPONENT. /OX. 6#H1N  THIS  PRORLCN  N  HAS  SEEN  SET  EQUAL  TO  .INP2JI4* 
3FA.SI  INPUU1S0 

so#  format  cofi#.st  ,  inpouisi 

5RI  FORMAT C/16A20ICDACJC  MALL  CONVtC TIONl (A9MOACK  HALLKX9IMESERV0IR/  1IW0U1S2 
113X23mC0EF  8TU/F  ISO-SFC-OlO  RBXi NCNISSI VI TV0X1 1HTOVE MATURE/  INPOU153 
21 7XF 10.6. 10XF 6.3. 1#XF1#.2I  INP0U1S6 

502  FORNATI5F10.5.9X.I1)  INP0U155 

DATA  OC.  ANN,  ASTER/ 1H  .IN*/  ITTOUISO 

InPUT-S  IIW0U1ST 

«OUt«»  IWOU1SO 

INCH >5  IMP0U159 

»!*«.«  INP0W100 

1)0  350#  1*1.5  INPOUIOI 

35##  NooFTlI)*#  INP0U102 


Input /output 

1  NPC-I 

WRITE  IA0UT.551INP6 

MEAD  C INPUT •  499.  E NO-999 I  IHCCORDIII.  Ial.341 
READ  (INPUT.  4991  (OCCOOOIII.  1*1.361 

00  TO  <9# 

IF  CEOF.  INPUTl  999.  990 
999  CONTINUE 
STOP 

990  CONTINUE 

WRITE  CKOUT.S#2I<R£CORO(II.I*1.36I 
•RITE  IROuT .503) 

WRITE  CXOUT (S#4I 
WRITE  CXOUT .SOS) 

WRITE  CKOUT.SM) 


INP0U163 

IIMU106 

IIW0U16S 

INP0U166 

INPOUIOI 

INP0UI40 

1NP0UI0* 

iirouiTo 

INP0U17I 

INP0UIT2 

1NP0U1T3 

EUH1TNINP0U1T6 

INP0U17S 

INP0U1I6 

INP0U1TT 

1NP0UIT0 

INPUU179 


read  iinput.saai  ar.rhooa.mmora.oa.psia.ea.traca.m.rhooo  .rnoms.oinpouio# 


IB.PSIb.U). TRACS  .C9.l)MOnC.(MORC.BC.  PSIC.EC.TRACC.NO (Ml  ItfNXIlBl 

mm  (Roui.st7)»9.«HOO».KMORA.iM.PSiA.t». r«ACA.09,RNooo.RHORB,B8.lNPoui82 
1PSIB.EB. TRACS  .C9.HH00C.hM0HC.8C.  PSIC.CC.THACC  imuil) 

«C*U  (INPUT. S03I  jr.MUWT.MT.NO.Nl.THZaO.TMFIN.UTPHT.OTPPZ.OTPRJ,  IMUIM 
10THO.bMP.TPR2.TPMJ.0ELH.DMl.0H2.0CLH6.6AMA.TZ  INP0UI85 

if t jf-n  in.iio.i7?  inpouim 

170  JF*2  IMMU107 

GO  10  172  INPOUIM 

171  JF«10  INPOUIM 

172  JFH» JF/2  INPOUIM 

Jf.jfH.JFH  IMPOU191 

JTMP*JfM.l  1NP0U192 

F  JFS*Jf  INP0U193 

FJFM-FJf  S/2.0  INP0U19* 

NOIsNU'Nl  1NP0U195 

IF  (NO)  181. 191. 182  INP0U196 

182  READ  (INPUT. 580)  (50(1)  .Ul.NO)  1NP0U197 

v  181  CONTINUE  INP0U198 

IF  (Nil  189.18*. 183  1NP0U199 

183  N0P*MU»1  1NP0U288 

M£U>  (INPUT. S80)  (SOU).  l«NOP.NOIt  INP0U291 

18*  1FI0ANAI  *08. *89. *09  1NP0U282 

*08  6AMA*RMOOC/(RHOOC*(RMOOA.NH008)-(RMOOA-RH006'I/6AI(AI  INP0U283 

*09  OnO*1.0-6ANA  1NP0U294 

NNO(1)*OANA*((4NOOA.RNOO(4)  .8NG*OHOOC  INPOU20S 

MmO(2)*6A*A*  (KMUfcA.RHORKI  .ONG-MMORC  1NP0U286 

OANAN*6ANA/(»NO(l)«lNNOOAmNOOO)  INP0U287 

■RITE  (KOUT.S10I6ANA.6AMAM  INP0U208 

IDNObU)  700.781.700  1M0U2M 

700  *RITC(K0UT. 79001  1NP0U210 

7900  FORMAT  (/23X.38H - OCCO^OSINO  BACK-UP  KINETICS - /I  INP0U211 

DO  78*  l-l.NOOU  INP0U212 

REAOIINPUT .S6*8IA9.RM00(I. II. BMOOII.I)  .800(1. 11.  PSIll.ll.EE4I.il.  1NP0U213 

1 TRAC  1 1. 1)  .89. 0*400(1 .2)  .8.400(1 .2)  .8M41.2).  PS1 II  .21  ,EE  II.2INPOU21* 

2I.TNACII.2I  .C9,KNOO(I,3|.8MOR(I,J>.BOO(I.3I.PSIII.JI.EEI1.3I.TRAC(1NPOU215 
31.31  1NP0U2I8 

DC  AO  (INPUT. 798 II  0HV(II,0HC(II.6A(I).rBCF(l>  INP0U217 

7901  FOKHAT  ( BOX. 2F 10. 5. 18X.2F 1 0.51  INP0U218 

MPITCIKOUT. 79821  I  1M0U219 

7982  FORMAT  (/28X.2«HUCC0MP0slM6  BACK-UP  NO.  .11)  INP0U228 

■KITE (ROUT. 50* >  INP0U221 

■RITE (KOUT. 5*7)  *9.RHOO|1.1I. RHOOII. ll.BaSll.D.PSIlI.lI.EEU.l).  INP0U222 
1  TRACl  1.1 1  .B9*RMuO(1.2)  .RM0RII.2)  .BMII.2).  PSII1.2I.EEII.2INP0U223 

21  .TRAC (I. 21  .C9.RH004l.3l  .RHUR(1.3)  .BM4 1.3)  .PSI  (1.31  .EC 1 1.3)  .TRACI IM0U22* 
31.31  UMU22S 

IF IGA (1)1  792.783.783  INP0U228 

702  GAIIl8RHOO(I.3l/INHOO(1.3l-  (RMOOll.l). RHOOII. 2)  l-IRHOOII.il*  1NP0U227 

1RH00I1.2)  l/GAIII  I  INP0U22S 

703  0HM(1I*1.-«A(1I  INP0U229 

RHOV  ( 1 l-GAI  D—IMrtOOlI.ll  .RHOOII  .2)  I  -ONBAI 1 )  *RNOO(1 .31  1NP0U238 

RHOC  llla(*A(I)*(KNOR(l.l|.RN0Rll  .2)  I  .ONBAII  )*RH0R(I.3I  INP0U231 

GAMAM*GA( I l/RHOV (l)*IRHnO(I .1 )  .RHOOII  .2)  )  IM0U232 

8R1TE  (ROUT  .5101  MIII.BaNAM  INP0U233 

P  1 1 )  *NH0F(I)/(RHUP(1  l-RHOCI  1 1 1  INP0U23* 

PP(I»-P(II*RHOC(II  1NP0U235 

DO  7038  J-1.3  1NP0U23* 

7030  FF41.J>*(  l.-PSI  (1.  Jl  I*0(IO(1.J1*IRMOOII.JI**I  l.-PSI  (1.01 1 1  INPUU237 

70*  CONTINUE  IM0U238 

701  CONTINUE  1NP0U239 

IFIDInU)  *1«.«19.*12  INPOU2*0 

*10  0TH8*5.  INP0U241 

*12  DTHlM.OI  INP0U292 

*13  1FITPM2I*1*«*I*.*15  INP0U2*3 

*1*  IPR2*IHFIN  1NP0U2** 

*15  IF(TPR3)*18.*18.*17  1*0(12*5 

*18  TPR3*THFIN  1M0U2** 

IF ( TPM2-THZRUI  *18.*18.*17  1NP0U2*7 

*18  DTPRI*0TPR2  1NP0U248 

OTPOZaOTPOJ  INP0U289 

TPR2>TPR3  INP0U2M 

80  TO  *16  1NP0U2S1 

*17  (MUTE  (ROUT. 511)  INP0U252 

TmF1N8ANAXIITHF1m.TP02.TPR3)  INP0U2SJ 

miTE  (KOUT.S12)TNZRO.T*1N  I*0U29* 

■RITE  IKOUT.5I3IOTPOT.TPR2  INP0U25S 

■RITE  (K0UI.5UI0TPO2.TpO2.TPR3  1M0U2M 

■RITE  (K0UT.S151UTPOJ.TpO3  1NPOU2S7 

■HITE  IR0UT.5191  OTNB  INP0U298 

IRACN8AM1N1 1  THACA.TRAC8.TRACC)  1NP0U299 

PCT£*HHOI|I/(HhO(1>-AMOi2)>  INP0U280 

PfT*PETE*RN0(2)  1*00261 

C  - — — —  NOOAL  PROPERTIES  INP0U282 

N.O  1NP0U28J 

RHS-10  INP0U28* 

NOBUCMbO  INPOUIM 

JI-JFh  INP0U288 

CALL  L COUNT  (-N4P4N-8.LCT  .NPO. RECORD  1 35) )  INP0U287 

382  DOITE  (ROUT  .5171  1*QUMS 

MITE  (KOUT.SIO)  INPOUIM 

■RITE  (KOUT.SIO)  1*00279 


B*ASItR  11*00271 

ME*DIINPOT.S«a>  IHATL II I  . TAI1I  .AREA  111  ,U£L  1 1 1  ,MA I l I  .RCI1I  , Irtl.NUMNI  11*00272 
A€««AI2I  INPOU273 

HSVMAtll  IMM027* 

“All I *0.0  11*00275 

00  *00  I«1.NUHN  INP0027* 

IF  1 1-21  *5*4.«S2.*53  INP0U277 

*51  RAIll«RAf  I-l>.|uu.ll-ll.i>CLIIII/2.«  11*00270 

60  TO  *01  11*00279 

*52  RAI2I-OCLI1I.6ELI2I/2.  1NP0U2M 

b.BLANA  HMUitl 

*61  OCL 1 1-1 1*0EL( 1*11/12.  INP0O202 

*5*1  RAVII >»RAIII  INP002S3 

IFIAll  *5*2. *5*3. *5*2  11*0029* 

*5*2  AREAlllxIANSIRSF.RAtll  M*vAE  1NP0U205 

60  TO  *5*  INP0O2S0 

*5*1  IF 1AREAI  In  *5**.*5**.*5*  INP0U2S7 

*5**  IFIMSv)  *5*5 .*5*0. *5*5  IMP0U2O* 

*5*5  AREA  I II sABSIRSV.HA I II I  INP0U299 

AE-1.0  UMIM 

60  TO  *5*  11*00291 

*5*0  AREA  111*1.  1NP0O292 

*5*  9RITE  (AOUT  .520)  I  .NATL 1 1 1  .TAI 1)  .AREA!  1}  .DEL 1 11  .RAIII  .W.ftCi  II  INP0U291 

IF (KAIL ( 11-21  *01.005.705  INP0029* 

*01  NL«1  1 1*00295 

00  *0*  J«1.J1  1NP0U290 

N»N»1  1  MR  1X129  7 

NOAINIsRMOOA  INP0O299 

rt0BlNI«RHO06  IKP00299 

*0*  ROCINIMIHOOC  UROUMO 

J1«JF  INP0U391 

60  TO  *00  11*00302 

*05  Ml. I  11*003*3 

no  *00  J*1.J1  IMII30* 

M.M*  1  11*001*5 

NOAINI-HHOMA  1NP0U3B0 

R0BIHI-RH0R8  INP0O3B7 

*00  M0CINI.RH0MC  11*003** 

J1«JF  INROUJ09 

60  TO  *00  INP0U11* 

705  IFIMATLIll-101  400.000.700  INPOOllI 

700  A.HAIitll-21  INP0U112 

V«M.flOAT(*)/2.«  INPOOllI 

ah.a/2  1NP0U11* 

IFIVKh-ahi  700.707.700  11*00115 

707  MOCbMll.llrtMHONIKH.il  1NP0U1I* 

HOC  OH  (1,21  *Rm0R  I  AH.  21  11*0031  T 

MOC0Mll.ll*RMOHIAN,3l  INPOOll* 

Atll-O.  11*00119 

60  TO  709  IMPOOJ2# 

700  am. AM. 1  11*00121 

MuCUMII.llrtHrtuOtAH.il  11*00122 

MOCUNI 1  .21  rtRrtOOf  Art.21  IIOOU3I3 

MutOHII.3lrtRrtOOIAM.3l  11*0012* 

Aill.i.o  11*00125 

70V  II  lArt-AHSl  7090.7092.7000  11*0022* 

7090  MFIIAHIrtl  11*00127 

AHS.An  \  11*00320 

MOOUtM.MOHOCH.  |  \  11*00129 

7092  MiAIAHTrtl  '■  11*011319 

*00  COMTI  VUC  INPOOllI 

IF  INOuo-MOoOCrtl  1 095 .709*. 7093  INP0O112 

7095  M41TE IKOOT .70901  INP00333 

7090  I 0MMA1 I 10A.55rtTOO  MONT  Dt COMPOS INK  BACK -OPS  IN  NOOAl  DATA  —  WIT  INP0O33* 

IJOHI  11*00335 

STOP  IMPW330 

7091  *01  Tt  IKOOT, 70971  11*00317 

7097  FOOHAI  llOI»5*rtTUO  FCl  OF  COMPOS  I  MO  BACK -OPS  IN  NOOAL  DATA  —  WIT  JINP0U31B 

I  OBI  INP00339 

SUP  11*00300 

709*  CONTINUE  IMP003AI 

*01  OClINuMNIrtOtllNUNNI/lf.  INP003A2 

CALL  5L0P0INUHN.MA.ANEA.EHAI  INP0UJ03 


IFIMOuUI  7090. 7099. T090 
7090  MOM2rtHLAIMOIVO|.l 
60  TO  T09| 

7099  MmM/.mOH 

70«1  CONTINUE 

IFIN59I  *031 ,*032.0011 
*031  RSVNrt-«5* 

■MITE IKOOT .55*INS0N.AC 
60  10  10* 

*012  901 II IKOOT .5551 
M  10  104 

•031  •OITllAOOT.S53IHSF.AE 
30*  MITE  IKOoT ,5211  OCLN.JF 
UELNrtOELN/U.O 

MCAUI IMPOT  .5021  HCON9.CP5N.INES.CNCNl.P7CNI.NeON 
II  ICrtCNIt  103,305.10* 


UNPOU30S 
11*001*0 
ll 1*003*7 
SI*OOJOD 
1NP001I9 
11*00350 
11MPOU3SII 


INPOUJSI 

1*0031* 

11*00150 

IMOM 

INPW3S7 

11*0035* 

11*00109 

11*003*0 


305  CHCB1.0.02  lIMUti 

30b  IF  IPlCNII  -HI.Ju7.JtB  1M0U302 

307  PrCF<I*0>9B  INP0UJ03 

304  T*INUl«t«ll>TueS  IIPWJM 

CALL  LCOUNT <9  .kCT.NP6.IK.COM>  1351 I  IM0U305 

■BUT  (ROUT .SOI >  rtCONV.EoSu.IBES  INPOUJO* 

CALL  L COUNI  IS  .LCT.NP0.NEC0MU3SI1  IW0UJO7 

MITT  |KGul.522l  1MUUJOO 

■BITE  IKOUT.S23IUNl.DN2.OU.nb  .  INP0UJO9 

CALL  LCOUNT  12  .LCT.NPb.BECONOtJSlI  INPOU370 

MITE  IKOUT. 529I1Z  ,  IMKM371 

IF(NUBU)  7vl0.7yl0.7S0 7  INP0U372 

7907  00  7900  3*1.NMU  INPOU373 

CALL  L COUNT  10.LCI  .NPb.BFCOBUI 351  >  INP0U379 

UNITE (KOUT. 7902!  J  IM0UJ75 

7909  F0NNAII/21A7NPLA>TlCllA*nCHA«17AJM6ASI  INP0UJ7A 

■KITE IKOUT. 79091  IMMU37T 

UNITE  IKOUT .523)  UNVIJt  .OHCIJI  .DELH6  IM0UJ70 

7900  MI1TEIK0UT.S29I  THEFIJI  INP0U379 

7910  CONTINUE  INPOUJOO 

- - MATERIAL  PA0PENT1ES  1M0U301 

CALL  LCOUNT 16  .LCT .NP6.NEC0M1I3 51!  IIPIH302 

310  «NITE  IKOUT. S2S>  IIWUN3 

CALL  LCOUNT  I3.LCT.NP6.HFCUMH35I  )  INPOUJO* 

UNITE  IKOUT. 7905)  INP0U3OS 

7905  FONMAI(/3A.42HueCONPOSlN6  BACK-UP  91061N  HAIER1ALS  22.29.2*.20.30.IM0UJO« 
1  CHAN  HATENIALS  23.2S.2T.29.31I  INPOUJO 7 

IT«0  IIPMUIO 

IL0I3I*!  1IP0U1« 

IL0I9) *1  INP0UJ90 

KT*1  iKPOUJOl 

350  IT*IT.l  1NP0UJ92 

BEAD  I1NPUT.S71I  NC.TT2llI.KT! .TCPllT .KT) .TMPI1T.KTI.TEPIIT.KTI  IIPM393 
IFINCI3S1.350.3S1  INP0UJ99 

351  Ini (KT*2)«1L0(KT«2! *1T-1  1NP0UJV5 

1BIkT*2!»IlOIKT*21  imoujo* 

TN2I1.KTI>0.  IMII397 

DO  357  1«2.1T  1ITOU390 

357  T HZ  1 1  «KT|  *TMZ  1 1*1  «KI  1. 1 TCPI I  »KT>. TCP  1 1*1  »KT)  1/2.*  I  TT2t  I  »KTI*fT2ll*IMQU399 


11. KID 

INPOUQOO 

CALL  LOOK  (KT*2.T2.TT2<T.KT!  .TH2I1.KT)  .O.O.O.HSM.OUN.l  1 

IMOUMI 

04  359  1*1.17 

IMMUtOE 

359  IhZI1.kI>«Ihzi I.kTI-hsh 

INPDW003 

CALL  LCOUNT  l*«ll  .LCI.NP6. BE  COM)  135)1 

lWWhM 

312  UNITE  IKOUT .5201  Kl.NHOIKTI  •  ITT2I1  .KT!  .TCPII.KT)  .TKPll.KT)  .THZll  .KTI 1MUU995 

1.TEPII.K)!.I>1.III 

UPOMO 

KI<KI«| 

iMOUbOT 

1 1*0 

1M9MM00 

IFINCI  3S0.3S3.JS3 

1M0U909 

35*  1« IKI-2)  3S0. 359.710 

INPOUOIO 

710  IF  INUOU)  711.9110.711 

1M0U91I 

711  1*1 

1NP0U012 

00  720  K*1.NU0U 

IMQU913 

17*0 

1NP0U91* 

712  1T»I7«I 

IM0U915 

BE *01 INPUT .571!  NC.TT5III.il. TCOUIlT.ll.TKOullT.il 

INPOUOIO 

IFINCI  713.712.113 

IM0U9IT 

713  1n1I1.2II*IT 

IM0U913 

lL0tl>21l*l 

IM0W019 

|NII*2II*1 

IMOUbM 

IfcNT 11.11*0. 

IM0U921 

04  7|.  3*2.11 

1M0U922 

71*  ICNTIj.llaTCNT I3*|.ll« ITC0UI3.II •TC0WI3-l.il )/2.9«TTSI3.1l*m<3-l|IMUkU 

1.ID 

IM0W929 

CALL  LOO*  121*1.  f*(F  Ml  .Tltll. II  .TENT  11.11. O.O.O.H5H.OIBI.  I) 

■■POMS 

DO  7IS  3*1.11 

IIPOU020 

715  ICNT<J.liaHNT<J.ll*«IN 

1NP0UO27 

|ri2*K*ll  717.710.717 

IMOUbM 

710  OnN«NhMCIKI 

IMOUbM 

00  TO  Tib 

IMOUbM 

71 7  000*<M09iKI 

IMM1 

7|0  CALL  LL0UNTI0*II.lCT.I#«.KCCON0I30I> 

1M0U.J2 

L*1«2I 

IMOUbJJ 

Ml  IE  IKOUT .52001  L.NBb.  ITI5I3.II  .TC0UI3.il  .T00UI3.il .TtNTI3.ll.  IM0M9J9 

13*1.111 

IIPOUOJO 

Ul«l 

IMOU.Jb 

11*0 

IIPWM7 

IFINCI  7I9.3S3.3S3 

,  1M0U.M 

719  IF  M*K*II  I20.7I2.7M 

IMOUbM 

7?0  CONTINUE 

IMOUbbO 

0||0  lu*0 

1  MOU.D 

L*0 

IM0UO42 

«|l  NCA0IIMVT.5O1!  KTMNOIkTI  .  IKNTLI II  .1*1*01 

1MOWOOJ 

IFIOMulKlII  35M.J5bO.JNK 

1MOU040 

35*0  LL»I  . 

IOPOU045 

L*L«I 

1WPUU000 

30UM 

1M0U.E7 

IF IKNILI II 1  JS»J.  JbSJi  3559 

1NP0WOO0 

Ml  *1*FL0ATILI 

1M0U0O9 

In*  l 

1M0UO50 

*-h 


tl*2  MT0U4SI 

JSS4  DO  I-LL.4  IIMKH 

iFIKMTUli)  JTCv.m4.MM  INPOUtSJ 

JTC5  jeO-JWJ.l  1M0U454 

4SVIJMUIXMTLIII  UT0U4SS 

J-WUIII  UNMM 

•NXJFTlJXl 
JTC4  CONTINUE 

1F4J0UI  35-4. 35.4, 3TCI  IMMH 

JSAI  CALL  OHOENI  ( JNU.KST)  INP0U4** 

»•  i*»*  imMi 

JTCI  U«I*.|  INP0U442 

HCM)  I  INPUT .4*7)  MC>  Mil  I. LI  .FI  1 11. LI  .72(11.1)  INP0U443 

IT  INC)  JTC2.JTCl.JTC2  INP0U444 

JTC2  IlO(J1«L>«1  INP0U445 

IhII3I<L)*M  INP0U4** 

IPUl-OM  INP0UV47 

CALL  LCOUNT  (5.lN.JM.4*il-(L«2)/4).LCT.M0«ll£C000t  J5)l  IM0U44* 

IriL-l)  J544.JS44.JTC7  INP0U444 

JS4*  MIIUUVI.4HI  DMHTf 

JS47  nolle  IK  OUT  .444)  L  IIPW4TI 

IF  (IN)  JS4J.J54J.JS42  MP0U472 

JS42  MITE  IKOUT .4*71)  .  INPOU47J 

JTCJ  IFIJNU)  JS4S.J545.JS44  *  INP0U47* 

JS44  MlIC  IKOUT  .4401)  IKSVIII.I-I.JOUI  IIMN7S 

V3S*5  MITt.IK0UT.444)  IT1|U)^||IA)^>II.I).M.MI  1M0U47* 

lM*i  1NP0U4T7 

_  F0BN4T(/25*«lM*.121.S(*t  I*I.1M.9NF2(M//IIM.J(5*.F10.4>))  IIT0U474 

445  F0PN4II//71.47N140LES  OF  dPTIOML  MSOKTIM  FUNCTIONS  FOK  TMCNlMOIMT* 


INP0W4M 

INP0U401 


1N4L  CONDUCT  IT  I  TT/25J.2JMK  •  FIIUNF  •  FtlUNKI 
44*  FORMAT  I/2JK.2DU -FUNCTION  TAME  NO.  .I1.12M  ASSItNED  TO) 

4471  F04MAI  IJ4I.IJHNA1N  MTEPIALI 
4401  FOKMTI2M.24MCCOMOSIN4  UCSHF  NO.  .11) 

IF  INC)  4ll.4lI.JS3 
355  IT*IT<1 

Of  AD  (INPUT. 571)  NC.TT2lll.KTI.TCPIIt.KTI.IKPIIt.NT> 

IF  INC)  JS4.3TC.JS4 
JS4  IlOIKWI.J 

InIIK|.2)*IlOIKW).|T-|  . 

■RlkT-2>*|L0<KT<2> 

CUL  LCOUNTITCII  .LCT.NP*.NCCONDf  35)  ) 

JI4  Will  IKOuT .527)  Kl .ONOlKI)  •  ITI2I  I.KT)  .TCPII.KT)  .TKPI1.KT)  .|4|,II)  1NP0U442 
1 7*4 

IF  INC)  4II.4II.JSJ 

c  - -  PTNIH.TSIS  Mi  CNTNKCPT 

J5J  NTI-K 
3*1  lN>l*Nll 
NT ■•••Nil 

KU  tlMUI.STs>NC.ITT|ll).l*IN.NTl>.<TN»(l).I*lN.NTl> 

1FINC)J*|.J*I.J42 
3*4  NT  INI  i-1 

342  IFlTIIINflll  JN4.J44.J4S 
J4S  IC0I2I.I 
1*12)41 
IhM2)4NT| 

Call  lCOuni(J4<i«ii.*i/si.lct.np*.occon»ijsi) 

JI4  Milt  IXOUl.SJJI 
IFN*# 

344  1N.II..I 

ifnminointi.ifnasi 
Mill  ikout.sjii«iiiiii.i*in.ifn  I 

Milt  IkOuT.SJJ)  (IN4I1I  .|4|N.IFNI 
IFINll-IFN)  J47.J4T.J40 

e  - -  F UNCI  IONS  OF  T INC 

3*7  NTM*4 
IS4# 

N0PI4* 

JTI  N|MM1M*| 

*E  *01 IMMT  .51 IINC.  7TNINTNI .  INC  INTMI  .TOP  ININ)  .TCMINIH)  .  IP  I  ININ)  l  WQH 
IIOPPINfNI  IMOU424 

IFiINmPINIhIIJTa.jTS.JT*  INPQU02I 

375  IUPPImIMIMPP  IM0U402 

J7*  l|«l  WMI) 

IFIICNINlNIIJ42.442.J4J  IN 

J42  11-2 

IrtfNt  ININI-2.IJ44.J44.J4J 
J44  !|«J 
34 J  IQPf  (NlNIall 

IFIII-ISIJ4S.J44.J45  IN 

345  noPImOPT.I  IN 

IS*  1 1  IN  I 

344  IFIHClJ7l.j7l.Jf2 
172  U0(ll*t 

InIIIIMTn  I 

ltl|l*IL0ll> 

CALL  lCOWNT  (4.J«NOPIMTNM.CIiMli*tCONOIJSI)  IN 

Mill  ISOul.SJAl  1M4U0J? 

15*4  IN 

00  J4I4  1*1. NlN 
11*100111) 


V 


*22 


IF  1 1 1-15)  3*7. 3*->. 3*7  imuMI 

1*7  IS- 1 1  IMUM 

60  10  13*71.  1*72.3*711.(1  IM0US*3 

3*71  *ftlt£(Ki)UT.S3-..l  IM0US4* 

GO  TO  3*7*  IM0US45 

3*72  MITE  (KCiUT. 5521  IMIMt 

60  to  3*75  IM0U547 

3*73  MItE(K0Ut.5s6)  1MUM 

60  10  3*75  INP0U94R 

3*9  60  TO  1 3*74.3*79. 3*751  *11  HMD3M 

3*7*  *OITE(KOUT.S36I  t  TNll)  .  I  l.T«(I)  .T0R(I>  .TCM(l)  .TPI  ( I)  .TRRP(l)  IM0U551 

60  TO  3*76  1M0USS2 

3*75  mlT-IROUT.5361  71(1(11  .11. Tl* III  .700(11  IMUM3 

3*76  ClMIlMbt  1NPUU55* 

IFINC-il  373.373.7*3  IMQU5S5 

373  00  3731  I-l.MTH  IM0U55* 

3731  TP1(1)-AL06(AMAA1(TP1(I|.. 666661)1  IM0U557 

CALL  L COUNT  12  .LCT.MC.HECOROOS))  1MU3M 

320  MITE  I  ROUT  .5371  IM 0US59 

C  - - - SUBfACE  ESUILISMlUN  TARLE  1NP0U5** 

PE  AO  (INPUT  .57961  OMS.KF2.Mea.UN.IR5T.MPF  .MIS. SMELL  INPOU561 

U1SNFIJ.I  IM0US62 

If  (NS Tl  29M.2966.29NI  IM0U563 

29*6  CMH-OMS  IM0U56* 

60  TO  2962  1NP0U565 

2961  1FIKMST-TT7)  2969.2963.2969  1M0U566 

2963  IFICMOMI  2967.2965.2967  1NPOU567 

2965  MITE  (ROUT .2966)  1 MOUSES 

CALL  (.COUNT 1 3.LCT.NP6.NF COHO < ]S> >  INP0U566 

2966  format  (//inasrhsurface  tamles  arc  the  sue  as  in  pnekious  pronlemi  imoustn 

60  TO  2912  IM0U5T1 

2967  MITE  (ROUT. 2966)  IM0U5T2 

29M  FORMAT I//I6A.72MREV10U5  SURFACE  T ARLES  CALLED  FOR  RUT  CM/CM  RAT  101NF0US73 

1  MAS  CMAMREO.  Null  JOS//)  IM0U57* 

STOP  IMOUSTS 

2669  MITE  (ROUT. 2916)  IMOUSTA 

2916  FORMAT  (//]**.  76MOCKIOUS  SURFACE  T ARLES  CALLED  FM  RUT  THIS  IS  FIRIMOUSTT 

1ST  PMORLEN.  SUIT  JOR//I  IMOUSTR 

STOP  IMOUSTR 

2962  RNST—7TT  I  MOUSE* 

2912  IF  (Ml  3266.3206.3203  IMOUSRI 

3261  If  l»S»)  2*2*. 2024.320*  IM0USR2 

3266  IF IRREAI  32 R*. 3264.3203  1M0USR3 

3263  M-l  I  MOUSE* 

60  TO  3201  1M0UBO9 

320*  IF  (MSI  I  126S.326S.3204  IMOUSRA 

320*  IF  IMSEMI  2613.2*61.2*13  IMOUSRI 

12*5  CONTINUE  IMOUSRO 

KLS—1  IMOUSR* 

NSEM—  I  IMOUSR* 

IP-1  IMOUSRI 

1PN-1  IM0USR2 

I«l  IMOUSRI 

Iwl  IMOUSR* 

J>«  IMOUSR* 

2**6  3-3*1  IMOUSR* 

IFINRPF)  2666 1. ROOM. 20M 1  IM0USR7 

20**6  RCAOIINCM.S761)  NSK«DRS.fLMC(3.1.1P>  .fTS(3.1.IPI  .6LO.TCMMI3.I.IP1IMOUSR6 
l.TSEMI  J)  •  JMS.ISUMF  (3)  IMOUSTR 

TRPFIJ.l.IPIM.  IMRU666 

00  TO  20*02  IMOUSRI 

2*661  REA0IIMCM.S7R6IKSK.hMS.TLMC 1 3.1  .IP) .TTSI3.1. IP) .HER. TCMNI  J.I.IPI  IMOU662 
I.TSCMIJI.JM.  T1MF  (3)  .TrPF  (3.1.IP)  IMQU6R3 

2RM2  COMTlNNE  IMOUSR* 

STM  FORMAT (3FS.S^«.*^t.].|F6.3.ll.2ltA*MI.E IR.3)  IM0U6RS 

IF  I JM6)  2*17.2*11.2*21  IMOUSR* 

2017  TSURF  (31MLAMI  IM0U6R7 

2*21  COMtlNUC  1MRU6M 

IF ITTSI3.I.1P) I  2RRl.tR32.2RRI  IMOUSR* 

2*01  IIS(3.I.IP1«IISI3«I.1PI*I.R  IM0U6I6 

ICMEMI3.I  .IP|6TCMEMU.I.IPI*I.R  1M6U6II 

TlC3M31*TtEN13l*I.S  IM0U6I2 

SR  10  2*65  IM0UR13 

2RR3  lTS(3.1.|P|6-T1k»3.l.|P|  IMSU6I* 

2665  IFI6LSI  2SR0.2RRt.IRST  IM0U61S 

IRRT  16*6  IMRU6I6 

IFIASHR.S1  2*2*. 20 1 1.262*  IM0U6I1 

2*69  RINAS  IMQM6IR 

20 1 1  IF  IMSEMI  2RR2.2626.I62R  I  MRU*  1 9 

20*2  IFIJM6I  2RM.2664.2R6*  I  MOM  626 

206*  MSCM63-1  IM0U62I 

I1EM4IPMUREM  IM0M622 

If  IMSLM-ll  2*26.2*26.20**  IM0U623 

20*0  DO  2*6*  L-IMREM  I  MRU*  2 1 

ttSEM(L.IPt*mu..l.lP>  IM0U62S 

trSEMIL.IPI*tCMJML.I.lPI  IMOUSR* 

2*6*  IMSEUIC.  IPl*tttRlL)  IMRMS27 

CAU.  SLSFOIMEM.IIKOII.IPI.tMICMII.IPI.ICPStlMI.IPII  1 1  POM  2* 

CALL  ILRPO(siEN.ItSEN||.IP» .tlEEMII.IPI .tCRSCM(I.IPI)  1M0US26 

Lit*  IMSEMI)  23*1  IM0US36 


•M>«CiM)lKII 


iril>l>ii2>ij).;(ij>2i|ji 
2013  CM.L  lC0UNflll.2-NM.LCl 
Ml  It  IKOUI.SJo! 

Milt  IKOu7.b7v7l  ClM.BLti.llfZ 
oo  loiTiiH.iiiiiirfris 
/•111  atlltMOIlT.Zalkl' 

2*112  roM«*i<M.*«Hri>MME  wmi  jsto  raa  suoace  cococr  teoos  mo 

IN6  CONNECTION) 

GO  10  MIU 

»IM  BBI It  IKOUT .20612) 

20131  CONTINUE 

If  (Ml  ZOI0.2B10.20I* 

2*1*  Ml  It  IK0UT.S79y>  OOEt 
60  10  201$ 
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29  iFiNA)  229.129.129 
229  IFIXIC-AKJI)  34.S3.33 
129  IFIX1<J9-X1C>  34.S3.33 
S3  3-1 

GO  TO  3 
4  J-J.l 
3  4*4-1 

44*4 

00  2  A-J.MM 
A)  (4.|)  —A)  (41 

1||N«1I«I|(K 
2  4*4-1 

I 1 (31-1 1C 
1  A|<J)«XIC 
R£TuRw 
140 


0R0CRD35 

OROCRO 34 


OR0CRO47 


C  UTILITY  9UOT14L  04©C*1  —  ROUTINE  TO  REORDER  INTEGER  ARRAYS 
SUt**0uU4r  OROCUKNX.IPmI 
0I4LNS104  1PHI74I 
H*HA 

IF14XI  9.4.1 
9  4«— 4 

1  IF (4-1 >6.6.2 

2  *  ■  § 

00  4  3-2.4 
1F(4A>».7.7 

7  IF  <U»4tJ-i>-IPH|J>)  4.4.3 
0  IF  «IPH(J-U-|R»1(J>|  3.3.4  x 

3  IT  •  |RH(>*I> 

IMH(J-I)  -  IPHIJ1 
|PH(J>  •  IT 

«  ■  1 

4  continue 

IF  (XI  4.4.2 
4  RETURN 
£40 


OROCRI  1 
OROCRI  2 
OROCRI  3 
0R0CR1  4 
OROCRI  S 
OROCRI  6 
OROCRI  7 
OROCRI  0 
OROCRI  9 
OROCRI 19 
OROCRI 11 
OROCRI 12 
0R0CR113 
OROCRI 14 
OROfJUlS 
OROCRI 14 
0R0CR117 
OROCRI 10 
OROCRI 19 
0R0CR129 
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C  UTILITY  BOUT  It*  StOU A  —  AEoROENS  UP  TO  TO UP  DEPENDENT  VARIABLE 
C  ARRAYS  VASCO  ON  IMC  RESULTS  OF  OROCRO  CAU. 

SUBROUTINE  SCOUaIN.L.A.R.C.OI 
01  At  Ml.  I  ON  AUl.alH.CUl.Dtll.Llll 
IS*# 

00  3#  I l*| »N 
1*11 

21  J*L 111 

Lin*i 

IriJ-II  22.30.22 

22  IF 4 1ST  2S.23.2S 

23  SAsAHI 

se*em  #r 

SC*C(1) 

SO*OUI 
IS*  I 

2*  A)II*AIJI 
Btll*MIJI 
C(I)*C(J) 

D(It«OIJI 

WJ 

00  TO  21 

2S  IFIIS-JI  26.20.2b 
20  IS*0 
A|II«SA 
B(1>«S0 
C|I)*SC 
0tlt«$0 
30  CONTINUE 
RETURN 
ENO 


StOU  A  | 
StOU  A  2 
St  QUA  3 
SC  QUA  4 
SEOUA  S 
SCOUA  6 
SEOUA  T 
SCOUA  0 
SEOUA  9 
SEOUA  10 
SEOUA  II 
SEOUA  12 
SEOUA  13 
SEOUA  lb 
SEOUA  IS 
SEOUA  16 
SEOUA  IT 
SEOUA  It 
SEOUA  19 
SEOUA  29 
SEOUA  21 
SEOUA  22 
SEOUA  23 
SEOUA  2b 
SEOUA  25 
SCOUA  26 
SCOUA  27 
SCOUA  20 
SCOUA  29 
SCOUA  3# 
SCOUA  31 


utility  routine  scouas  —  Reorders  up  to  five  dependent  variable 

SCOUAS  1 

ARRAYS  MASCO  ON  THE  results  of  OROEOO  CALL 

SCOUAS  2 

SOB* OUT I**  SCOUASIN.L.A.O.C.O.E) 

SCOUAS  3 

DIMENSION  All I. Ml 11 «C(lt. Dill. LI 11 

SCOUAS  - 

dimension  cut 

SCOUAS  s 

IS*# 

SCOUAS  6 

00  39  11*1. N 

SCOUAS  7 

1*11 

SCOUAS  0 

21 

J*LU» 

SCOUAS  9 

Llll-I 

SCOUAS  1 A 

IFIJ-II  22. 3#. 22 

scouas i; 

22 

1FIISI  2S.23.2S 

SC  QUASI  2 

23 

SA*AI 1 1 

SE QUASI J 

SB*M( 1 1 

scauAsib 

SC*Cill 

SCOUASIS 

SONHII 

SCOUAS lb 

SC*Cltt 

SCOUAS 17 

!S*1 

SC QUASI# 

2* 

AlIIMIJI 

SC0UAS19 

Milt«aiJt 

SC0UAS2# 

C|It*CIJI 

SC0UAS21 

0111*0131 

SCOUAS22 

EIII*EU> 

K0UAS23 

l*J 

SEOUAStb 

00  TO  21 

ICOUAStS 

25 

iriIS>Jt  26.20.26 

SC0UAS26 

?• 

is*# 

SCQUAS27 

A|II*bA 

SC0UAS20 

out«so 

SC0UAS29 

CllfSC 

SC0UAS3# 

uin*so 

SC0UAS3I 

Elll*SC 

SC0UAS32 

3* 

CONTINUE 

RETURN 

END 

UQUU39 
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c  utility  buutihc  slow  —  bchekates  slopes  or  dependent  vaoiam.es 
C  HaSEO  ON  OUAiMAIlC  firs  TmKjMN  three  SUCCESSIVE  iamle  points 
SUBROUTINE  SLOPUINUHA, A.P.EA) 

c  slope  evaluation  MuriNr 
c 

01 HENS 10*  Alll,  Pill,  FNIII,  nil 
SO  En|2I  ■  I  PI2I  -  Pill  I  /  |  *121  -  XIII  I 
ENI1I  ■  ENI2I 
/HIM.* 

OC  «  ENIII 
00  SO  I  •  I  ,  M MX 
IPO  •  I  •  1 
IPT  «  I  •  2 
IT  *  IPO  -  NUN* 
ir  till  S3  ,  J1  .  32 
si  os«oc 

GO  TO  «1 
32  BO  TO  AO 
S3  EOT  •  XIII  -  XI1P0I 
ATT  •  A I IPO)  -  AI1PIJ 
ATO  •  AIIPTI  •  Alll 
AA  *  Pill  /  I  AOl  *  ATO  I 
XOTT*XOT«XTT 
ST  AS*PI1*II/A0TT 

AC  •  PI IPT I  /  I  XTT  •  XTO  > 

AAA  •  AA  •  ATI 
AM  «  AO  •  ATO 
ACC  •  AC  •  XOT 
0*  •  OC 
OH  *  ENI|| 

OC  •  EN  IIPOI 

CniIPoi  •  AO  •  IAOT  -  ATTI  •  ACC  -  AAA 
ENIIPII  ■  AC  •  IATT  -  XTO  I  •  AAA  •  AM 
ENf II*  AA  •  I A  TO  -  AOTI  •  AOS  •  ACC 
3*  0E  »  ENIII 

iril-21  3*,*1,3S 

35  EnIII  *  I  OC  •  uA  I  ft. 

*1  EniII  •  lENIII  •  gi  I  /l, 

*0  A0*AIII-A|I-1I 
3*  CONTINUE 
RETURN 
CNO 


SLOP* 

SLOP* 

SLOPS 

SLOPS 

SLOPS 

SLOP* 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 


SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 

SLOPS 


/ 


*-J2 


,-ah. 


m-V 

,  *  it &«*«>  •  :*  .  • 
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